t>-  •  -\r 


AD-A256  766 

IlllBlilllHIlllil 


1  IX  VJ  i.1  JL J  ±—J  W  *  A  **  A  **  *  -*» ' 

MIXING  LAYERS 

bv 

L.-K.  Tsene.  P.-K.  Wu  and  G.  M.  Faeth 


15  August  1992 


College  of  Engineering 

DEPARTMENT  OF  AEROSPACE  ENGINEERING 

The  University  of  Michigan 
Ann  Arbor,  Michigan  48109-2140 


DISCLAIM!  NOTICE 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTIC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
COLOR  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY  ON  BLACK 
AND  WHITE  MICROFICHE. 


GDL/GMF-92-02 


HIGH-SPEED  MULTIPHASE 
MIXING  LAYERS 


by 

L.-K.  Tseng,  P.-K.  Wu  and  G.  M.  Faeth 


15  August  1992 


Prepared  for: 


department  of  the  Navy 
Office  of  Naval  Research 
Mechanics  Division,  Code  1132P 
800  N.  Quincy  Street 
Arlington,  Virginia  22217-5000 


Grant  No.  N00014-89-J-1199 
Gabriel  D.  Roy 
Scientific  Program  Officer 


92-27688 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  This  FACE  rWJien  Dote  Enlmrmd) 


|  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  REFORT  NUMBER 

2.  GOVT  ACCESSION  NO. 

3  RECIPIENT’S  catalog  number 

4.  TITLE  (end  Submit) 

High-Speed  Multiphase  Mixing  Layers 

s  type  of  report  a  period  covered 

Final  Report 

11/15/88  -  6/30/92 

«•  performing  org.  report  number 

GDL/GMF-92-02 

7.  AUTHOR^ 

L.-K.  Tseng,  P.-K.  Wu  and  G.  M.  Faeth 

«.  contract  or  grant  number^ 

N00014-89-J -1199 

>■  FCRFORMINO  ORGANIZATION  NAME  ANO  AOd'rESS 

Department  of  Aerospace  Engineering 

The  University  of  Michigan 

Ann  Arbor,  MI  48109-2140 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

It.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Naval  Research 

800  North  Quincy  Street 

Arlington,  VA  22217-5000 

12.  AEPOHT  DATE 

August  1991 

13.  number  of  pages 

i«.  MONITORING  AGENCY  NAME  A  ADDRESS  (II  dlllerent  from  Controlling  Olllce) 

is.  security  class,  f oi  o> i*  report.) 

Unclassified 

1So.  DECLASSIFICATION'  DOWNGRADING 
SCHEDULE  . 

IB.  DISTRIBUTION  STATEMENT  (ol  thl <  Kopwi)  - - - 

Approved  for  public  release;  distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  (ol  Ihm  ebetrect  entered  In  Block  JO.  II  dlllerent  Iron,  Roport) 

IB.  SUFFLEMENTARY  NOTES 

Multiphase  flow,  primary  breakup,  atomization 

*0.  ABSTRACT  ( Continue  on  rororot  elde  II  neceeearr  mtd  Identity  by  block  number) 

This  is  the  final  report  of  a  research  program  whose  objective  was  to 
study  the  structure  and  liquid  breakup  properties  of  high-speed  multiphase 
mixing  layers.  The  research  is  relevant  to  both  gas  injection  into  liquids, 
e.g.,  liquid  metal  combustion  processes,  and  liquid  injection  into  gases, 
e.g.,  pressure-atomized  sprays.  The  investigation  was  divided  into  two 
phases,  considering  the  structure  of  multiphase  mixing  layers  and  primary 
breakup  at  liquid  surfaces,  respectively. 

DO  (  j ah*7j  1473  COITION  OF  I  NOV  «S  IS  OBSOLETE 

s/n  0102-  lf-  ou-  460 1  i ;  — Un c l.as s  i f  i  ed _ 

11  SECURITY  CLASSIFICATION  OF  This  RASE  fBBort  0«o  Mntered) 


Unclassified 

SKCURITY  CLASSIFICATION  OF  THIS  PACK  (Whit  Dm • 

20.  Abstract  (continued) 


The  structure  of  multiphase  mixing  layers  was  studied  using  large  round  liquid 
jets  (9.5  mm  diameter)  in  still  air  at  pressures  of  1-8  atm;  with  both  slug  flow  and 
fully-developed  turbulent  pipe  flow  at  the  jet  exit.  Measurements  included  mean 
liquid  volume  fractions  using  gamma-ray  absorption  and  dispersed-phase  properties 
using  single-and  double-pulsed  holography.  It  was  found  that  mixing  was  affected 
strongly  by  the  liquid/gas  density  ratio  and  the  degree  of  flow  development  at  the  jet 
exit,  with  the  smallest  gas/liquid  density  ratios  and  fully-developed  turbulent  pipe 
flow  yielding  the  fastest  mixing  rates.  Predictions  of  flow  structure  based  on  the 
locally-homogeneous  flow  approximation  (where  velocity  differences  between  the 
phases  are  neglected)  were  reasonably  good  at  large  liquid  volume  fractions  but 
progressively  overestimated  the  rate  of  flow  development  as  liquid  volume  fractions 
decreased  and  liquid/gas  density  ratios  increased.  Liquid  volume  fractions  were  low 
in  the  dispersed  flow  region  so  that  primary  and  secondary  breakup  dominated  flow 
properties  rather  than  collisions.  Separated  flow  effects  were  important  in  the 
dispersed  flow  region,  but  tended  to  decrease  as  liquid/gas  density  ratios  decreased 
due  to  reduced  drop  sizes  from  secondary  breakup. 

Primary  breakup  was  studied  along  the  surface  of  large  round  liquid  jets  (3.6-9.5 
mm  diameter,  with  water,  n-heptane  and  glycerol  mixtures)  in  still  gases  (helium,  air. 
Freon  12)  at  pressures  of  1  and  2  atm,  with  both  slug  and  fully-developed  turbulent 
pipe  flow  at  the  jet  exit.  Measurements  included  flash  photography,  pulsed 
shadowgraphy  and  single-  and  double-pulsed  holography.  Drop  sizes  after  primary 
breakup  satisfied  Simmons'  universal  root  normal  distribution  and  could  be 
characterized  by  the  Sauter  mean  diameter  alone.  The  behavior  of  nonturbulent 
primary  breakup  could  not  be  explained  using  existing  unstable  wave  growth 
theories  but  could  be  correlated  by  a  mechanism  based  on  the  vorticity  within 
boundary  layers  along  the  passage  walls;  the  drops  formed  were  intrinsically 
unstable  to  secondary  breakup  which  appears  to  be  the  main  aerodynamic  effect  for 
these  conditions.  When  liquid/gas  density  ratios  were  greater  than  500,  primary 
breakup  of  turbulent  liquids  was  controlled  by  liquid  turbulence  properties  alone, 
yielding  successful  correlations  for  both  onset  conditions  and  the  subsequent 
evolution  of  drop  sizes  after  primary  breakup  with  distance  along  the  mixing  layer. 
Aerodynamic  effects  become  important  at  lower  liquid/gas  density  ratios  and 
reduced  the  distance  required  for  onset  and  drop  sizes  at  onset  by  enhancing  the 
growth  of  small  scale  surface  distortions.  Drops  sizes  after  primary  breakup 
downstream  of  the  onset  condition  were  reduced  by  aerodynamic  effects  as  well,  due 
to  merging  of  primary  and  secondary  breakup.  Both  these  effects  on  turbulent 
primary  breakup  were  correlated  successfully  using  phenomenological  theories. 


S/N  0102-  LF.0M.6601 

iii  Unclassified 


SCCuniTY  CLASSIFICATION  OF  THIS  PAO*fW*M  DM  CfltWMQ 


ACKNOWLEDGMENTS 


The  authors  wish  to  acknowledge  the  assistance  of  the  staff  of  the  Gas  Dynamics 
Laboratories  of  the  Department  of  Aerospace  Engineering:  D.  R.  Glass,  W.  C.  Eaton,  G. 
L.  Gould,  T.  L.  Griffin,  D.  Huddleson  and  T.  M.  Larrow  for  help  with  the  experiments; 
Debbie  Laird  and  S.  C.  Bauerle  for  help  with  publications  and  financial  management. 


T""'  "”  ™CTED  l 


iv 


|  Aooesslon  For  | 

HTIS 

GRA&I 

Sf 

DTIC 

TAB 

□ 

Uiianr; 

ounce  d 

□ 

Ju~ti 

r  t  not.  1  on 

Bv. 

Dl  st  r 

J  rut.  ion/ 

Av'>  1 

Codes 

,A  i’ll  and/or 

Dist 

> 

|  Special 

i  i 

TABLE  OF  CONTENTS 


Page 

ACKNOWLEDGMENT .  iv 

LIST  OF  TABLES .  vi 

LIST  OF  FIGURES .  vii 

NOMENCLATURE .  ix 

1.  INTRODUCTION .  1 

2.  STRUCTURE  OF  MULTIPHASE  MIXING  LAYERS .  6 

2.1  Introduction .  6 

2.2  Experimental  Methods .  7 

2.3  Theoretical  Methods .  7 

2.4  Results  and  Discussion .  10 

2.5  Conclusions .  21 

3.  PRIMARY  BREAKUP  AT  LIQUID  SURFACES .  21 

3.1  Introduction .  21 

3.2  Experimental  Methods .  22 

3.3  Nonturbulent  Liquid  Breakup .  24 

3.4  Turbulent  Liquid  Breakup .  32 

3.5  Conclusions .  43 

REFERENCES .  46 

APPENDIX  A:  FAETH  (1990) .  49 

APPENDIX  B:  RUFF  ET  AL.  (1992) .  54 

APPENDIX  C:  TSENG  ET  AL.  (1992a) .  64 

APPENDIX  D:  TSENG  ET  AL.  (1992b) .  73 

APPENDIX  E:  WU  ET  AL.  ( 199 1 ) .  93 

APPENDIX  F:  WU  ET  AL.  (1992) .  104 

APPENDIX  G:  WU  AND  FAETH  (1992) .  126 


v 


LIST  OF  FIGURES 

Figure  Caption  Page 

1  Sketch  of  gas  injection  into  liquids .  2 

2  Sketch  of  test  apparatus .  8 

3  Sketch  of  holocamera .  9 

4  Liquid  volume  fractions  along  axis  for  fully-developed  flow .  11 

5  Liquid  volume  fractions  along  axis  for  slug  flow .  12 

6  Radial  profiles  of  liquid  volume  fractions  for  fully-developed  flow .  13 

7  Radial  profiles  of  liquid  volume  fractions  for  slug  flow .  14 

8  Drop  size  distributions .  15 

9  Drop  properties  for  fully-developed  flow  at  x/d  =  25  and  1  atm .  17 

10  Drop  properties  for  fully-developed  flow  at  x/d  =  25  and  4  atm .  18 

1 1  Drop  properties  for  slug  flow  at  x/d  =  25  and  1  atm .  19 

1 2  Drop  properties  for  slug  flow  at  x/d  =  25  and  4  atm .  20 

13  Sketch  of  test  apparatus .  23 

14  Sketch  of  holocamera . \ .  25 

15  Streamwise  variation  of  SMD  for  nonturbulent  breakup .  26 

16  Correlation  of  SMD  for  fully-developed  nonturbulent  primary  breakup 

using  existing  theories .  28 

17  Sketch  of  injector  passage  boundary  layer  vorticity  mechanism .  29 

1 8  Correlation  of  SMD  for  fully-developed  nonturbulent  primary  breakup 

based  on  the  boundary  layer  vorticity  mechanism .  30 

19  Mass  averaged  drop  velocities  after  turbulent  primary  breakup .  33 

20  Sketch  of  aerodynamically-enhanced  turbulent  primary  breakup  at  the 

liquid  surface .  34 

21  SMD  at  the  onset  of  turbulent  primary  breakup .  36 

22  Length  to  initiate  turbulent  primary  breakup .  37 

23  Sketch  of  turbulent  breakup  processes .  39 


Vll 


LIST  OF  FIGURES  (continued) 

Eigms  Caption  Page 

24  SMD  after  turbulent  primary  breakup  when  aerodynamic  effects  are 

small .  40 

25  SMD  after  aerodynamically-enhanced  turbulent  primary  breakup .  42 

26  Turbulent  primary  breakup  regime  map .  44 


NOMENCLATURE 


Cb  constant  in  nonturbulent  primary  breakup  correlation 

Cs  empirical  constant  for  SMD  after  secondary  breakup 

Csa  empirical  constant  for  enhanced  breakup 
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Up  drop  streamwise  velocity 
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vp  drop  crosstream  velocity 

Wejj  Weber  number  based  on  density  of  phase  i  and  length  j,  p^j  u2/o 


x 


streamwise  distance 


Greek 


a 

P 

A 

P 

a 

Subscripts 


volume  fraction 
molecular  viscosity 

turbulent  crosstream  integral  length  scale,  d/8  for  fully-developed 
turbulent  pipe  flow 

density 

surface  tension 

aerodynamic  secondary  breakup  time 
Rayleigh  breakup  time 


f 

FD 

g 

i 

o 

Superscripts 


liquid  phase  property 

fully-developed  nonturbulent  primary  breakup  regime 
gas-phase  property 
property  at  onset  of  breakup 
jet  exit  condition 


(  )  time-average  quantity 

(~)  Favre-averaged  quantity 

O’  root-mean-squared  time-averaged  fluctuation 


1.  INTRODUCTION 


This  investigation  considered  the  properties  of  high-speed  multiphase  mixing  layers 
involving  gases  and  liquids.  Such  flows  are  encountered  in  a  number  of  applications  of 
interest  to  the  Navy  and  can  involve  injection  of  either  gases  into  liquids  or  liquids  into 
gases.  Typical  applications  of  gas  injection  into  liquids  include  liquid  metal  combustion  for 
underwater  propulsion  systems,  direct  contact  condensers,  metal  processing  systems,  gas 
dissolution  systems,  reservoir  destratificaticn  systems  and  nuclear  reactor  pressure 
suppression  systems,  among  others.  Typical  applications  of  liquid  injection  into  gases 
primarily  involve  numerous  liquid  atomization  processes. 

The  present  investigation  seeks  a  better  understanding  of  both  gas  and  liquid 
injection  processes,  however,  it  is  a  direct  outgrowth  of  an  earlier  ONR-sponsored  study  of 
gas  injection  into  liquids  (Chuech  et  al.  1989,  Loth  and  Faeth  1988,  1989,  1990).  This 
earlier  investigation  was  motivated  by  interest  in  developing  the  stored  Chemical  Energy 
Propulsion  System  (SCEPS)  for  underwater  propulsion  (Bush  and  Merrill  1981;  Hughes  et 
al.  1981).  This  process  involves  injection  of  gaseous  sulfur  hexafluoride  into  a  high- 
temperature  molten  bath  containing  molten  lithium  and  the  liquid  combustion  products  of 
the  reaction  between  sulfur  hexafluoride  and  lithium.  Unsteady  injector  flow  at  low  gas 
flow  rates  causes  clogging  of  injector  passages  because  the  liquid  combustion  products 
solidify  at  relatively  high  temperatures,  ca.  1050-1 100K.  Thus,  the  injectors  usually 
operate  as  underexpanded  jets  which  yields  a  steadier  flow  that  eliminates  the  clogging 
problem.  Pressures  within  the  reaction  vessel  primarily  are  controlled  by  the  vapor 
pressure  of  the  liquid  reactant/product  mixture  within  the  molten  bath  and  usually  are 
subatmospheric,  however,  they  can  vary  over  a  relatively  wide  range  depending  on  the 
power  level  and  the  extent  of  reaction. 

The  actual  SCEPS  combustion  process  is  a  difficult  environment  for  detailed 
measurements;  therefore,  most  studies  of  the  process  have  involved  analogous  processes  of 
noncondensing  and  condensing  gas  jets  into  liquids  (Chan  et  al.  1987;  Chen  and  Faeth 
1982;  Kemey  et  al.  1972;  Loth  and  Faeth  1988,  1989,  1990;  Weimer  et  al.  1973).  The 
general  picture  that  emerges  from  these  studies  is  illustrated  in  Fig.  1,  which  specifically  is 
based  on  air  injection  into  water  from  Loth  and  Faefh  (1988,  1989,  1992).  Basic  to  the 
SCEPS  process,  the  flow  leaves  the  passage  as  a  soni^  gas  jet.  Expansion  to  the  ambient 
pressure  is  initiated  by  an  expansion  fan  at  the  periphery  of  the  jet  exit,  which  accelerates 
the  bulk  of  the  gas  flow  to  supersonic  velocities.  The  expansion  waves  reflect  as 
compression  waves  from  the  sonic  line  region  of  the  multiphase  mixing  layer  at  the  flow 
periphery;  these  compression  waves  subsequently  steepen  to  form  shock  waves.  The  shock 
waves  cross  the  flow  and  then  are  reflected  as  expansion  waves  from  the  sonic  line  region 
of  the  mixing  layer,  tending  to  repeat  this  pattern  in  a  succession  of  shock  cells.  The  shock 
cells  eventually  disappear,  ending  the  external  expansion  region,  when  the  sonic  line  within 
the  mixing  layer  reaches  the  axis.  Clearly,  the  process  is  unusually  complicated,  involving 
patches  of  supersonic  flow,  multiphase  flow,  liquid  breakup  into  drops,  gas  breakup  into 
bubbles,  and  reactive  phenomena  that  yield  a  number  of  intermediate  solid  and  immiscible 
liquid  phases  (Avery  and  Faeth  1975;  Chan  et  al.  1988;  Chen  and  Faeth  1983;  Loth  and 
Faeth  1988,  1989,  1990). 

Even  though  the  SCEPS  process,  and  other  noncondensing  and  condensing  gas 
injection  processes  like  it,  is  very  complex,  some  progress  has  been  made  toward  modeling 
the  flow  for  engineering  purposes  (Avery  and  Faeth  1975;  Chan  et  al.  1987,  1988;  Chen 
and  Faeth  1982,  1983;  Chuech  et  al.  1989;  Faeth  1987,  1990;  Loth  and  Faeth  1988,  1989, 
1990).  These  methods  generally  simplify  the  treatment  of  multiphase  flow  and  reaction 
effects  by  adopting  the  locally-hornogeneous-flow  (LHF)  approximation,  which  implies 
local  kinematic  equilibrium  (velocities  of  all  the  phases  locally  the  same)  and 
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Figure  1  Sketch  of  gas  injection  into  liquids 


thermodynamic  equilibrium  (temperatures  of  all  phases  the  same  and  the  chemical  potential 
of  each  species  equal  in  all  phases).  Turbulent  mixing  generally  is  treated  using  a 
turbulence  model,  with  recent  approaches  employing  Favre-averaged  k-e-g  models. 
Finally,  the  complexities  of  the  external  expansion  region  have  either  been  treated  by  full 
gas-dynamic  analyses  or  approximated  as  an  equivalent  adapted  jet  (if  the  gas  or  vapor 
containing  region  is  large  in  comparison  to  the  external  expansion  region).  Models  of  this 
type  have  been  reasonably  successful  for  treating  gross  features  like  the  length  of  the  gas  or 
vapor  containing  region,  when  evaluated  using  measurements  (Chen  and  Faeth  1982, 
1983;  Loth  and  Faeth  1988,  1989,  1990).  However,  more  detailed  features  of  the  flow, 
like  the  structure  and  liquid  breakup  properties  of  the  multiphase  mixing  layer,  that  are 
crucial  to  the  development  of  the  flow,  have  not  been  predicted  very  satisfactorily  using  the 
LFtF  approach.  This  deficiency  that  has  been  attributed  to  effects  of  separated  flow,  e.g., 
lack  of  kinematic  and  thermodynamic  equilibrium  between  the  phases  (Loth  and  Faeth 
1988,  1989,  1990). 

Even  though  the  gas  dynamic  features  of  underexpanded  gas  jets  in  liquids  are 
complicated,  it  is  generally  within  the  capabilities  of  current  technology  to  treat  these 
effects.  In  contrast,  the  properties  of  the  multiphase  mixing  layer  along  the  periphery  of  the 
flow  are  not  well  understood  (Faeth  1987,  1990).  In  particular,  the  rate  of  generation  of 
drops  at  the  surface,  which  controls  mixing  rates  between  the  phases,  and  the  size  and 
velocities  of  the  drops  after  breakup,  which  control  the  importance  of  separated  flow 
effects,  are  unknown.  Thus,  further  progress  to  develop  more  effective  models  of  gas 
injection  processes  primarily  depends  on  a  better  understanding  of  the  multiphase  mixing 
layer;  therefore,  studying  this  portion  of  the  flow  was  the  primary  objective  of  the  present 
investigation. 

Several  features  of  multiphase  mixing  layers  must  be  addressed  for  applications 
reverent  to  gas  injection  processes.  First  of  all  the  liquid  phase  generally  is  turbulent  so  that 
effects  of  liquid  turbulence  on  the  structure  of  the  mixing  layer  and  primary  breakup  of 
liquid  at  the  liquid  surface  must  be  known.  Secondly,  gas  injection  systems  operate  over  a 
relatively  wide  range  of  liquid/gas  densities  ratios  and  effects  of  this  important  variable 
must  be  considered.  Additionally,  properties  like  liquid  viscosity  and  surface  tension  vary 
widely  as  well,  and  their  effects  must  be  addressed.  Finally,  the  velocity  difference  across 
the  mixing  layer  is  known  to  affect  the  structure  of  the  mixing  layer  and  this  behavior  must 
be  understood.  These  features  all  were  specifically  addressed  during  the  present 
investigation. 

In  view  of  the  present  limited  understanding  of  multiphase  mixing  layers,  however, 
it  also  was  necessary  to  neglect  some  features  of  gas  injection  processes  in  order  to  keep 
the  study  tractable.  First  of  all,  the  presence  of  gas  dynamic  features  like  shock  and 
expansion  waves  in  the  mixing  layer  are  substantial  complications;  therefore,  such  effects 
were  not  studied  pending  the  development  of  background  information  on  the  simpler 
constant  pressure  multiphase  mixing  layer.  One  feature  of  the  flow  that  helps  to  mitigate 
lack  of  consideration  of  strong  gas  dynamic  features,  like  shock  waves,  is  the  fact  that  the 
sonic  line  is  located  at  some  distance  of  the  surface,  so  that  the  most  important  part  of  the 
mixing  layer  near  the  liquid  surface  only  involves  subsonic  flow  with  no  direct  interactions 
with  shock  waves.  The  next  simplification  was  to  reverse  the  process  to  consider  a  high¬ 
speed  liquid  flow  in  still  gas,  rather  that  a  high-speed  gas  flow  over  a  relatively  still  liquid. 
This  was  done  to  illuminate  problems  of  overall  flow  instability  (like  those  encountered  at 
low  flow  rates  for  the  SCEPS  system)  that  would  complicate  study  of  mixing  layer 
properties,  and  to  provide  well-defined  liquid  turbulence  properties.  This  allowed  provision 
of  various  relative  velocities  across  the  mixing  layer  but  is  only  appropriate  to  the  extent 
that  flow  properties  are  independent  of  the  source  of  the  momentum.  Fortunately,  for 
multiphase  mixing  layers,  separated-flow  effects  are  the  dominant  feature  and  the  direction 


of  the  velocity  vector  with  respect  to  the  liquid  has  little  effect  on  primary  breakup  at  the 
liquid  surface.  Due  to  the  density  difference,  however,  the  overall  structure  of  the 
multiphase  mixing  layer,  is  affected  by  the  direction  of  the  velocity  vector;  nevertheless, 
gaining  understanding  of  conditions  where  the  momentum  is  carried  by  the  liquid  still  is  a 
reasonable  first  step  toward  understanding  gas  injection  processes.  Naturally,  the  present 
experimental  configuration  yields  information  that  is  directly  applicable  to  pressure 
atomization,  particularly  the  near-injector  region. 

This  is  a  final  report  of  accomplishments  during  the  investigation.  The  present 
report  is  brief,  additional  details  can  be  found  in  other  articles,  papers  and  theses  generated 
during  the  course  of  the  study,  see  Table  1  for  a  summary  of  these  publications.  In 
addition,  the  articles  by  Faeth  (1990),  Ruff  et  al.  (1992),  Tseng  et  al.  (1992a,b),  Wu  et  al. 
(1991,  1992)  and  Wu  and  Faeth  (1992)  are  appended  to  provide  a  convenient  reference. 
Finally,  Tseng  (1991)  and  Wu  (1992)  provide  detailed  descriptions  of  the  research  as  well 
as  complete  tabulations  of  data. 


Table  1.  Summary  of  Publications 


Faeth,  G.  M.  (1990)  Structure  and  atomization  properties  of  dense  turbulent  sprays. 
Twenty-Third  Symposium  ('International-)  on  Combustion,,  The  Combustion 
Institute,  Pittsburgh,  1345-1352. 

Chuech,  S.  G.,  Lai,  M.-C.  and  Faeth,  G.  M.  (1989)  Structure  of  turbulent  sonic 
underexpanded  free  jets.  AIAA  J.  27,  549-559. 

Loth,  E.  and  Faeth,  G.  M.  (1989)  Structure  of  undeTexpanded  round  air  jets  submerged  in 
water.  Int.  J.  Multiphase  Flow  15,  589-603. 

Ruff,  G.  A.,  Wu,  P.-K.,  Bernal,  L.  P.  and  Faeth,  G.  M.  (1992)  Continuous-and 
dNpersed-phase  structure  of  dense  nonevaporating  pressure-atomized  sprays.  L. 
Prop.  Power  8,  280-289. 

Tseng,  L.-K.,  Ruff,  G.  A.  and  Faeth,  G.  M.  (1992a)  Effects  of  gas  density  on  the 
structure  of  liquid  jets  in  still  gases.  AIAA  J.  30,  1537-1544. 

Tseng,  L.-K.,  Wu,  P.-K.  and  Faeth,  G.  M.  (1992b)  Dispersed-phase  structure  of  pressure 
atomized  sprays  at  various  gas  densities,  J.  Prop.  Power,  in  press. 

Wu,  P.-K.  and  Faeth,  G.  M.  (1992)  Aerodynamic  effects  on  primary  breakup  of  turbulent 
liquids.  Atomization  and  Spravs.  submitted. 

Wu,  P.-K.,  Ruff,  G.  A.  and  Faeth,  G.  M.  (1991)  Primary  breakup  in  liquid/gas  mixing 
layers.  Atomization  and  Spravs  1,  421-440. 

Wu,  P.-K.,  Tseng,  L.-K.  and  Faeth,  G.  M.  (1992)  Primary  breakup  in  gas/liquid  mixing 
layers  for  turbulent  liquids.  Atomization  and  Spravs.  in  press. 

Papers: 

Faeth,  G.  M.,  Tseng,  L.-K.  and  Wu,  P.-K.  (1990)  Structure  and  mixing  properties  of 
high-speed  gas/liquid  mixing  layers.  27th  JANNAF  Combustion  Meeting, 
Cheyenne,  WY. 
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Faeth,  G.  M.,  Tseng,  L.-K.  and  Wu,  P.-K.  (1991)  Primary  breakup  of  round  pressure 
atomized  sprays.  28th  JANNAF  Combustion  Meeting,  San  Antonio,  TX. 

Ruff,  G.  A.  and  Faeth,  G.  M.  (1990)  continuous-and  dispersed-phase  structure  of  dense 
nonevaporating  pressure-atomized  sprays.  AIAA  28th  Aerospace  Sciences  Meeting, 
Reno,  NV,  AIAA  Paper.  No.  90-0464. 

Tseng,  L.-K.  and  Faeth,  G.  M.  (1991)  Dispersed-phase  structure  of  pressure-atomized 
sprays  at  various  gas  densities.  ICLASS-91,  NIST,  Gaithersburg,  MD. 

Tseng,  L.-K.,  Ruff,  G.  A.  and  Faeth,  G.  M.  (1991)  Effects  of  ambient  gas  density  on  the 
structure  of  pressure-atomized  sprays.  AIAA  29th  Aerospace  Sciences  Meeting, 
Reno,  NV,  AIAA  Paper  No.  91-0690. 

Tseng,  L.-K.,  Wu,  P.-KJ.  and  Faeth,  G.  M.  (1992)  Dispersed-phase  structure  of  pressure 
atomized  sprays  at  various  gas  densities.  AIAA  30th  Aerospace  Sciences  Meeting, 
Reno,  NV,  AIAA  Paper  No.  92-0230. 

Wu,  P.-K.  and  Faeth,  G.  M.  (1991)  Primary  breakup  of  round  turbulent  liquid  jets  in  still 
air.  ICLASS-91,  NIST,  Gaithersburg,  MD. 

Wu,  P.-K.  and  Faeth,  G.  M.  (1993)  Aerodynamic  effects  on  primary  breakup  of  turbulent 
liquids.  31st  AIAA  Aerospace  Sciences  Meeting,  Reno,  NV,  AIAA  Paper  No.  93- 
0903. 

Wu,  P.-K.,  Ruff,  G.  A.  and  Faeth,  G.  M.  (1991a)  Primary  breakup  in  liquid/gas  mixing 
layers.  AIAA  29th  Aerospace  Sciences  Meeting,  Reno,  NV,  AIAA  Paper  No.  91- 
0285. 

Wu,  P.-K.,  Tseng,  L.-K.  and  Faeth,  G.  M.  (1991  Breakup  of  round  pressure- atomized 
sprays.  28th  JANNAF  Combustion  Meeting,  San  Antonio,  TX. 

Wu,  P.-K.,  Tseng,  L.-K.  and  Faeth,  G.  M.  (1992)  Primary  breakup  in  gas/liquid  mixing 
layers  for  turbulent  liquids.  AIAA  30th  Aerospace  Sciences  Meeting,  Reno,  NV, 
AIAA  Paper  No.  92-0462. 

Wu,  P.-K.,  Hsiang,  L.-P.  and  Faeth,  G.  M.  (1993)  Aerodynamic  effects  on  primary  and 
secondary  spray  breakup.  First  International  Symposium  on  Liquid  Rocket 
Combustion  Instability,  The  Pennsylvania  State  University,  University  Park,  PA. 

Theses: 

Tseng,  L.-K.  (1991)  Near-injector  structure  of  nonevaporating  pressure-atomized  sprays  at 
various  ambient  densities.  Ph.D..  thesis.  The  University  of  Michigan,  Ann  Arbor, 
MI. 

Wu,  P.-K.  (1992)  Liquid  surface  breakup  for  nonturbulent  and  turbulent  liquids.  Ph.D. 
thesis.  The  University  of  Michigan,  Ann  Arbor,  MI. 


The  investigation  was  divided  into  two  phases,  one  considering  the  overall  structure 
of  multiphase  mixing  layers,  and  the  other  specifically  treating  primary  breakup  processes 
at  the  liquid  surface.  Highlights  of  each  of  these  phases  of  the  research  are  discussed  in  the 
following,  in  turn. 


2.  STRUCTURE  OF  MULTIPHASE  MIXING  LAYERS 

2.1  Introduction 

This  phase  of  the  investigation  involved  study  of  multiphase  mixing  layers  formed 
along  the  sides  of  nonevaporating  liquid  jets  injected  into  still  gases.  Experiments  were 
limited  to  the  atomization  breakup  regime,  where  liquid  breakup  into  drops  and  the 
development  of  a  multiphase  mixing  layer,  begins  right  at  the  jet  exit  because  this  regime  is 
most  important  for  practical  applications  such  as  gas  injection  into  liquids  and  pressure- 
atomized  sprays.  Jet  exit  conditions  included  both  fully-developed  turbulent  pipe  flow  and 
low  turbulence  intensity  slug  flow  in  order  to  examine  the  effect  of  liquid  turbulence  on 
flow  structure.  The  measurements  were  made  using  large  diameter  (9.5  mm)  round  liquid 
jets  in  order  to  extend  the  region  where  the  structure  of  the  multiphase  mixing  layer  could 
be  observed;  however,  the  complete  breakup  properties  of  these  flows,  as  the  entire  liquid 
column  disintegrates,  was  not  considered.  Effects  of  various  liquid/gas  density  ratio  were 
studied  by  injecting  water  into  still  air  at  pressures  in  the  range  1  -8  atm.  The  measurements 
were  used  to  evaluate  predictions  of  flow  structure  based  on  the  LHF  approximation. 

Discussion  of  past  work  will  be  limited  because  reviews  treating  multiphase  jets  and 
dense  sprays  have  appeared  recently  (Faeth  1987,1990).  Past  measurements  have 
established  the  main  features  of  the  near  injector  region  of  liquid  jets  in  still  gases  within  the 
atomization  breakup  regime.  The  flow  consists  of  a  liquid  core,  much  like  the  potential  core 
of  a  single-phase  jet,  surrounded  by  a  dispersed  drop/gas  mixing  layer.  The  liquid  core 
becomes  shorter  as  liquid/gas  density  ratios  decrease  due  to  increased  entrainment  rates 
(Chehroudi  et  al.  1985).  Other  features  of  these  flows,  however,  have  only  be  studied  for 
injection  of  water  into  air  at  atmospheric  pressure  using  relatively  large  diameter  (19.1  and 
9.5  mm)  round  liquid  jets  (Ruff  et  al.  1988,  1991).  These  measurements  included  liquid 
volume  fraction  distributions  using  gamma-ray  absorption,  drop  sizes  and  velocities  using 
double-pulse  holography  and  air  entrainment  rates  using  laser  velocimetry.  The 
measurements  of  mean  liquid  volume  fraction  distributions  showed  that  the  rate  of 
development  of  the  mixing  layer  increased  with  increased  liquid  turbulence  levels. 
Predictions  based  on  the  LHF  approximation  were  effective,  including  effects  of  liquid 
turbulence,  for  liquid  volume  fractions  greater  than  0.2,  however,  they  progressively  failed 
as  the  flow  became  dilute.  This  difficulty  with  the  LHF  approach  was  found  to  be  due  to 
significant  separated  flow  effects  in  regions  where  liquid  volume  fractions  were  low. 
Additionally,  because  entrainment  rates  mainly  are  governed  by  processes  within  dilute 
portions  of  the  flow,  the  LHF  predictions  generally  overestimated  entrainment  rates  due  to 
separated  flow  phenomena.  However,  the  predictions  suggested  improved  performance  of 
LHF  methods  as  the  Wedf  increased,  with  this  limit  approached  more  rapidly  as  the 
liquid/gas  density  ratio  decreases. 

The  objective  of  the  present  investigation  was  to  explore  effects  of  liquid/gas 
density  ratio  on  the  mixing  properties  of  liquid  jets  in  still  gases  and  the  adequacy  of 
predictions  of  flow  properties  based  on  the  LHF  approximation.  Experiments  involved 
measurements  of  liquid  volume  fraction  distributions  within  large  round  water  jets  (9.5  mm 
initial  diameter)  in  still  air  at  pressures  of  1-8  atm.  The  measurements  were  compared  with 
predictions  based  on  the  LHF  approach  to  determine  whether  this  methodology  could  treat 
effects  of  varying  liquid/gas  density  ratios. 
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The  discussion  of  the  study  begins  with  a  brief  description  of  experimental  and 
theoretical  methods.  Results  then  are  considered,  treating  liquid  volume  fractions  and 
dispersed-phase  properties  in  turn.  The  present  description  is  brief,  see  Tseng  (1991),  Ruff 
et  al.  (1992)  and  Tseng  et  al.  ( 1992a, b)  for  more  details.  For  convenience,  the  latter  three 
papers  can  be  found  in  Appendices  B,  C  and  D. 


2.2  Experimental  Methods 

Apparatus.  Experimental  methods  were  similar  to  Ruff  et  al.  (1988,  1991)  except 
that  the  flow  was  contained  in  a  large  pressure  vessel  so  that  the  pressure  of  the  ambient  air 
could  be  changed.  A  sketch  of  the  apparatus  appears  in  Fig.  2.  The  arrangement  involves  a 
steady  water  jet  injected  vertically  downward  within  a  1.5m  diameter  x  4.5  m  high  pressure 
vessel.  Instrumentation  was  mounted  rigidly  so  that  flow  structure  was  measured  by 
traversing  the  injector.  Water  collected  at  the  bottom  of  the  chamber  and  was  drained 
continuously.  Two  injectors  were  used,  identical  to  those  of  Ruff  et  al.  (1988,  1991);  one 
yielding  nonturbulent  slug  flow,  the  other  yielding  fully-developed  turbulent  pipe  flow. 

Instrumentation.  Distributions  of  liquid  volume  fractions  were  measured  using 
gamma-ray  absorption  similar  to  Ruff  et  al.  (1988).  This  was  based  on  an  iodine- 125 
isotope  source  within  a  lead  casket  having  a  small  (1.6  mm  diameter)  outlet  aperture.  The 
gamma  rays  passing  through  the  fiow  were  detected  through  a  second  lead  aperture  using  a 
Bicron  x-ray  probe  and  an  EG&G  Ortec  single-channel  analyzer  and  counter  timer. 
Absorption  measurements  (20,000-25,000  counts)  were  made  for  30-60  parallel  paths 
through  the  flow;  these  results  then  were  deconvolved  to  find  the  radial  variation  of  time- 
averaged  liquid  volume  fraction. 

The  appearance  of  the  liquid  surface  and  the  properties  of  drops  in  the  dispe.sed- 
flow  region  were  observed  using  single-and  double-pulsed  holography.  A  sketch  of  the 
holocamera  appears  in  Fig.  3.  It  was  based  on  an  off-axis  arrangement  using  a  ruby  laser 
with  a  pulse  duration  of  20  ns  (and  pulse  separations  as  small  as  1  p.s)  to  stop  the  flow. 
This  arrangement  allowed  objects  as  small  as  2  pm  to  be  observed  and  drops  having 
diameters  as  small  as  5  pm  to  be  measured,  while  double  pulsing  allows  drop  velocities  to 
be  measured.  Measurements  were  obtained  over  6x6x4  mm3  volume,  using  at  least 
three  holograms  per  position.  The  holograms  were  reconstructed  and  analyzed  using  a  high 
magnification  video  system  in  conjunction  with  a  Gould  FD  5000  Image  Processor. 

Test  Conditions.  Mean  flow  conditions  were  the  same  for  slug  and  fully-developed 
flows.  This  involved  a  water  flow  rate  of  3.47  kg/s,  an  average  jet  exit  velocity  of  49.1 
m/s,  Refd  =  462,000,  Ohd  =  0.00121,  Wefd  =  312,000  and  Wegd  =  380,760,  1520  and 
3040  at  pressures  of  1 ,  2,  4  and  8  atm. 


2.3  Theoretical  Methods. 

Flow  properties  were  predicted  using  the  LHF  model  of  Ruff  et  al.  (1988,  1991). 
In  addition  to  the  LHF  approximation,  the  major  assumptions  of  the  model  were  as 
follows:  steady  (in  the  mean)  axisymmetric  flow  with  no  swirl,  boundary-layer 
approximations  apply,  negligible  kinetic  energy  and  viscous  dissipation  of  the  mean  flow, 
buoyancy  only  affects  the  mean  flow,  equal  exchange  coefficients  of  all  species  and 
phases,  and  negligible  mass  transport  between  the  phases  (no  evaporation).  These 
assumptions  are  either  conditions  of  the  experiments  or  are  justified  by  past  practice  except 
for  the  LHF  approximation  which  will  be  evaluated. 


WATER  INLET 


Figure  2  Sketch  of  test  apparatus 
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1800  mm 


Figure  3  Sketch  of  holocamera 


Under  these  approximations  the  flow  field  can  be  found  from  a  simplified  version 
of  the  conserved-scalar  formulation.  This  involved  solving  Favre-averaged  governing 
equations  for  conservation  of  mass,  stream  wise  mean  momentum,  mean  mixture  fraction, 
turbulence  kinetic  energy,  and  the  rate  of  dissipation  of  turbulence  kinetic  energy.  The 
specific  formulation,  all  empirical  constants,  calibration  of  the  approach  for  variable-density 
single-phase  jets,  the  specification  of  initial  and  boundary  conditions,  and  the  details  of  the 
numerical  computations,  can  be  found  in  Ruff  et  al.  (1988,  1991). 

2.4  Results  and  Discussion 

Mean  Liquid  Volume  Fractions.  Measured  and  predicted  time-averaged  mean  liquid 
volume  fractions  along  the  jet  axis  are  plotted  in  Figs.  4  and  5  for  fully-developed  and  slug 
flow  jet  exit  conditions,  respectively.  Predictions  for  slug  flow  are  shown  for  passage 
length  to  diameter  ratios,  L/d  =  0  and  5,  which  bounds  the  potential  degree  of  flow 
development  within  the  injector  passage.  Measurements  of  Ruff  et  al.  (1988)  for  similar 
test  conditions  at  1  atm  also  are  shown  on  the  plot:  they  agree  very  well  with  present 
measurements.  The  substantial  effect  of  jet  exit  turbulence  levels  on  mixing  rates  is  clearly 
evident,  with  fully-developed  flow  yielding  the  fastest  mixing  rates.  The  faster  mixing  rates 
as  liquid/gas  density  ratios  decrease  (or  pressure  increases)  also  is  indicated  by  the 
measurements. 

The  LFIF  predictions  are  in  reasonably  good  agreement  with  the  measurements  for 
fully-developed  flow  (Fig.  4)  and  provide  the  correct  trend  of  effects  of  ambient  pressure. 
However,  the  conditions  of  Fig.  4  represent  relatively  low  levels  of  mixing.  For  example, 
although  mean  liquid  volume  fractions  as  small  as  0.03  are  reached,  the  entire  range 
considered  still  represents  Favre-averaged  mixture  fractions  greater  than  0.85  due  to  the 
relatively  large  liquid/gas  density  ratios  of  the  flow  (Tseng  et  al.  1992a).  There  are 
significant  differences  between  predictions  for  pure  slug  flow  (L/d  =  0)  and  allowance  for 
boundary  layer  growth  within  the  injector  passage  (L/d  =  5),  see  Fig.  5.  These  conditions 
bound  the  range  of  possibilities  for  present  tests,  and  it  is  encouraging  that  the  two 
predictions  tend  to  bound  the  measurements.  Nevertheless,  there  are  significant 
discrepancies  between  predictions  and  measurements  at  x/d  =  100  and  4  and  8  atm.  This 
region  involves  rapid  streamwise  variation  of  flow  properties  and  dilute  flows  where  both 
predictions  overestimate  the  rate  of  development  of  the  flow  —  behavior  that  can  be 
attributed  to  effects  of  separated  flow  based  on  subsequent  results. 

Predicted  and  measured  radial  profiles  of  mean  liquid  volume  fractions  at  1  atm  are 
illustrated  for  fully  developed  and  slug  flow  in  Figs.  6  and  7  (results  at  other  pressures  are 
qualitatively  similar).  The  results  of  Ruff  et  al.  (1988)  also  a  plotted  on  the  figures:  they 
generally  agree  with  present  measurements  within  experimental  uncertainties.  The 
comparison  between  predicted  and  measured  mean  liquid  volume  fractions  is  reasonably 
good,  however,  the  measurements  generally  involve  large  mixture  fractions  as  noted 
earlier.  In  particular,  as  mixture  fractions  become  small  at  large  x/d,  predictions  clearly 
overestimate  rate  of  flow  development,  indicating  significant  effects  of  separated  flow  as 
the  flow  becomes  dilute. 

Drop  Properties.  The  failure  of  the  LHF  approximation  within  the  dispersed  flow 
region  can  be  proven  directly  based  on  measurements  of  drop  properties.  An  important 
simplification  was  observed  with  respect  to  drop  size  properties  that  is  illustrated  in  Fig.  8. 
This  is  a  plot  of  drop  size  distributions  according  to  Simmons’  (1977)  universal  root 
normal  distribution  for  various  positions  along  the  liquid  surface  and  across  the  multiphase 
mixing  layer  considering  both  fully-developed  and  slug  slow  conditions.  Remarkably,  the 
data  generally  satisfies  this  distribution  function  with  the  ratio  MMD/SMD  =  1.2:  a  few 
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Figure  4  Liquid  volume  fractions  along  axis  for  fully-developed 
flow 
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Figure  5  Liquid  volume  fractions  along  axis  for  slug  flow 
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Figure  6  Radial  profiles  of  liquid  volume  fractions  for 
fully-developed  flow 
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Figure  8  Drop  size  distributions 
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exceptions  at  large  cumulative  volume  percentages  are  felt  to  be  due  to  sampling  limitations 
for  large  drops.  This  is  a  two-parameter  distribution;  therefore,  taking  the  best  estimate 
value,  MMD/SMD  =  1.2,  implies  that  the  entire  distribution  can  be  represented  by  the  SMD 
alone.  Thus,  drop  size  results  will  be  reported  in  terms  of  SMD  in  the  following. 

Typical  distributions  of  drop  properties  across  the  mixing  layer  for  fully-developed 
flow,  illustrating  flow  properties  for  turbulent  primary  breakup,  appear  in  Figs.  9  and  10 
(see  Tseng  (1991)  and  Tseng  et  al.  (1992b)  for  results  at  other  test  conditions).  These 
measurements  are  for  x/d  =  25  at  pressures  of  1  and  4  atm.  Properties  illustrated  include  the 
variation  of  ellipticity  (mass  averaged  ratio  of  maximum  to  minimum  drop  dimensions), 
SMD  and  drop  velocities  across  the  dispersed  flow  region.  Predictions  of  flow  velocities 
using  the  LHF  approach  also  are  shown  on  the  plots. 

Many  of  the  features  observed  in  Figs.  9  and  10  are  similar  to  the  earlier 
observations  of  Ruff  et  al.  (1991)  at  atmospheric  pressure,  however,  significant  effects  of 
reducing  the  liquid/gas  density  ratio  are  observed  as  well.  Thus,  at  atmospheric  pressure, 
the  largest  values  of  ep  and  SMD  are  in  the  region  of  the  liquid  surface,  supporting  the 
presence  of  primary  breakup  yielding  large  drops  which  subsequently  undergo  secondary 
breakup.  At  the  higher  pressure,  however,  ep  approaches  unity  and  drop  sizes  are  smaller 
near  the  surface,  with  smaller  variations  of  these  properties  across  the  width  of  the  flow: 
this  suggests  merging  of  primary  and  secondary  breakup  processes  near  the  liquid  surface 
as  the  liquid/gas  density  ratio  decreases.  In  general,  there  are  substantial  differences  among 
the  velocities  of  drops  of  different  size,  which  is  direct  proof  of  the  deficiency  of  the  LHF 
approximation  in  this  region  (the  LHF  predictions  also  clearly  are  not  very  good  for  any 
drop  size).  Increasing  pressures  tend  to  reduce  drop  velocities  but  the  velocities  of  small 
drops  are  surprising  low  and  suggest  small  gas  velocities.  This  implies  that  momentum 
exchange  between  the  liquid  and  gas  is  not  very  effective.  This  occurs  due  to  the  large 
relaxation  time  of  the  large  drops  that  carry  most  of  the  liquid  momentum,  aided  by  the  fact 
that  liquid  volume  fractions  in  the  dispersed  flow  region  are  very  low  (less  than  0.1%). 
Thus,  this  region  intrinsically  is  subject  to  separated  flow  effects,  with  relatively  large 
drops  moving  a  large  velocities  in  a  dilute  multiphase  flow. 

The  structure  of  the  dispersed  phase  and  the  effect  of  varying  liquid/gas  density 
ratio  is  quite  different  for  slug  flow.  Some  typical  results  for  x/d  =  25  and  pressures  of  1 
and  4  atm  are  illustrated  in  Figs.  1 1  and  12.  First  of  all,  drop  sizes  are  clearly  much  smaller 
for  slug  then  for  fully  developed  flow  (cf.  Figs.  9  and  10,  and  11  and  12)  even  though 
fully-developed  flow  mixes  more  rapidly  as  a  whole.  Values  of  en  and  SMD  are  largest 
near  the  surface,  suggesting  effects  of  secondary  breakup  as  before,  however,  drop 
ellipticity  and  SMD  across  the  mixing  layer  are  far  more  uniform  for  slug  flow  than  for 
fully-developed  flow.  Decreasing  the  liquid/gas  density  ratio  also  has  much  less  effect  on 
drop  properties  for  slug  than  for  fully-developed  flow.  However,  drop  velocities  vary 
significantly  with  drop  size,  suggesting  significant  effects  of  separated  flow  as  before. 

Taken  together,  these  results  demonstrated  rather  different  primary  breakup 
properties  for  turbulent  and  nonturbulent  liquids  but  important  effects  of  separated  flow  in  a 
dilute  dispersed  flow  region  for  both.  Reduced  liquid/gas  density  ratios  for  water/air  flows 
at  pressures  of  1-8  atm  also  yield  evidence  of  aerodynamic  breakup  effects  with  drop  sizes 
after  primary  breakup  particularly  reduced  for  turbulent  liquids.  The  key  to  treating 
separated-flow  phenomena  is  an  understanding  of  the  outcome  of  primary  breakup  (drop 
size  and  velocity  distributions);  therefore,  this  issue  was  addressed  in  some  detail  during 
the  second  phase  of  the  investigation,  which  will  be  considered  next. 
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Figure  10  Drop  properties  for  fully-developed  flow 
at  x/d  =  25  and  4  atm 
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Figure  12  Drop  properties  for  slug  flow  at  x/d  =  25 
and  4  atm 
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2.5  Conclusions 


The  structure  and  mixing  properties  of  nonevaporating  multiphase  mixing  layers 
along  the  sides  of  liquid  jets  in  still  gases  were  studied  for  various  liquid/gas  density  ratios. 
The  major  conclusions  of  the  study  are  as  follows: 

1 .  Decreasing  liquid/gas  density  ratios  increases  the  rate  of  mixing  and  the  flow  width, 
analogous  to  the  behavior  of  single-phase  mixing  layers. 

2.  Increasing  turbulence  levels  in  the  liquid  increase  the  mixing  rates  and  flow  width 
as  well. 

3.  Use  of  the  LHF  approximation  yielded  good  results  for  mean  liquid  volume 
fractions  greater  than  0.1  including  predictions  of  effects  of  liquid  turbulence. 
However,  these  conditions  correspond  to  relatively  low  levels  of  mixing  (mixture 
fractions  greater  than  0.85)  due  to  the  large  liquid/gas  density  ratios  of  the  present 
flows. 

4 .  Effects  of  separated  flow  were  observed  for  mean  liquid  mass  fractions  less  than 
0.85,  causing  poor  performance  of  predictions  based  on  the  LHF  approximation  in 
this  region. 

5.  Local  drop  size  distributions  throughout  the  dispersed  flow  region  satisfied 
Simmons'  (1977)  universal  root  normal  distribution  function  with  MMD/SMD  = 
1 .2,  and  can  be  characterized  completely  by  the  SMD  alone. 

6 .  Drop  sizes  after  primary  breakup  were  larger  for  turbulent  than  laminar  liquids  even 
though  the  rate  of  mixing  is  more  rapid  for  turbulent  liquids. 

7.  Effects  of  decreasing  the  liquid/gas  density  ratio  on  drop  sizes  after  primary 
breakup  were  relatively  small  for  nonturbulent  liquids  but  caused  drop  sizes  to 
progressively  decrease  for  turbulent  liquids. 

8 .  Drops  sizes  were  largest  at  the  liquid  surface  while  mean  liquid  volume  fractions 
within  the  dispersed  phase  were  less  than  0.1%;  this  implies  that  secondary  breakup 
is  a  dominant  feature  of  the  mixing  layer,  rather  than  collisions. 

Thus,  pending  resolution  of  existing  uncertainties  about  primary  and  secondary 
breakup,  the  LHF  approach  offers  a  useful  treatment  of  high-speed  multiphase  mixing 
layers.  It  should  be  recognized,  however,  that  the  LHF  approach  will  over-estimate  the  rate 
of  development  of  the  flow  by  a  degree  that  can  not  be  quantified  until  uncertainties  about 
breakup  are  quantified.  This  process  will  be  treated  next. 


3.  PRIMARY  BREAKUP  AT  LIQUID  SURFACES 

3.1  Introduction 

The  objective  of  this  phase  of  the  investigation  was  to  study  primary  breakup  along 
the  liquid  surface  in  multiphase  mixing  layers.  Based  on  the  results  discussed  in  Section  2, 
the  investigation  was  divided  into  two  parts,  considering  primary  breakup  of  nonturbulent 
and  turbulent  liquids.  However,  aerodynamic  effects  were  considered  in  both  cases. 
Experimental  methods  were  similar  to  the  study  of  multiphase  mixing  layers,  mainly 
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involving  observations  of  primary  breakup  along  the  surfaces  of  liquid  jets  in  still  gases 
using  single-and  double-pulsed  holography. 

A  number  of  reviews  of  liquid  atomization  processes  have  appeared  (Clift  et  al. 
1978;  Faeth  1987,  1990;  Harrje  and  Reardon  1972;  Lefebvre  1980;  Levich  1962; 
McCarthy  and  Malloy  1974;  Reitz  and  Bracco  1982);  therefore,  only  the  main  features  of 
past  work  will  be  considered.  Beginning  with  nonturbulent  primary  breakup,  past  work 
has  sought  to  explain  mean  drop  sizes  after  primary  breakup.  Two  major  approaches  have 
been  proposed:  (1)  aerodynamic  breakup  based  on  linear  wave  growth  theory,  which  was 
proposed  initially  by  Taylor  (1940),  with  subsequent  development  by  Ranz  (1959)  and 
Levich  (1962)  and  later  application  by  Reitz  and  Bracco  (1982),  among  others;  and  (2) 
sheltered  wave  growth  theory  due  to  Mayer  (1961)  and  Adelberg  (1968)  which  is  based  on 
the  surface  wave  growth  model  of  Jeffreys  (1926).  Intrinsic  to  both  approaches  is  that 
primary  breakup  is  induced  by  instabilities  growing  on  the  liquid  surface  due  to 
aerodynamic  effects  in  spite  of  extensive  evidence  of  strong  effects  of  vorticity  in  the  liquid 
at  the  passage  exit  on  primary  breakup  properties  (Faeth  1990,  Wu  et  al.  1991).  Thus, 
these  methods  have  not  been  very  successful  and  have  only  been  evaluated  using  data  that 
was  probably  affected  by  secondary  breakup,  yielding  a  long  history  of  continued  use  due 
to  the  absence  of  alternatives  and  limited  acceptable  ranges  of  application  (Harrje  and 
Reardon  1972).  Strong  effects  of  liquid  vorticity  on  spray  breakup  have  been  reported  for 
some  time.  First  of  all,  the  early  studies  of  pressure  atomization  of  DeJuhasz  et  al.  (1932) 
and  Lee  and  Spencer  (1933)  showed  that  atomization  quality  differed  for  nonturbulent  and 
turbulent  liquids.  Next,  Grant  and  Middleman  (1966)  and  Phinney  (1966)  observed  that  jet 
stability  and  the  onset  of  breakup  were  affected  by  the  presence  of  turbulence.  Additionally, 
Chehroudi  et  al.  (1985)  observed  effects  of  jet  exit  turbulence  on  the  length  of  the  liquid 
core  for  pressure  atomization.  Finally,  Hoyt  and  Taylor  (185)  found  virtually  no  effect  of 
breakup  properties  when  gas  velocities  were  varied  at  atmospheric  pressure  for  fixed  liquid 
flow  conditions:  clear  proof  that  most  of  the  existing  data  base  at  atmospheric  pressure  was 
not  affected  by  aerodynamic  phenomena  in  spite  of  continued  use  of  aerodynamic  theories 
to  correlate  measurements  over  limited  data  ranges. 

The  present  investigation  seeks  to  clarify  some  of  the  difficulties  of  current 
understanding  of  primary  breakup.  The  approach  was  to  make  measurements  of  primary 
breakup  properties,  exploiting  recent  progress  in  the  use  of  holography  for  measurements 
of  this  type  (Ruff  et  al.  1991).  This  was  accompanied  by  the  use  of  phenomenological 
theories  in  an  attempt  to  explain  and  correlate  the  measurements,  in  order  to  capture  the 
physics  of  the  process  with  a  tractable  investment  in  analysis  and  computations. 

The  discussion  begins  with  a  brief  description  of  experimental  methods.  Results  are 
then  considered,  treating  the  findings  of  measurements  and  phenomenological  analysis  for 
the  primary  breakup  of  nonturbulent  and  turbulent  liquids  in  turn.  The  present  description 
is  brief,  see  Wu  (1992),  Wu  and  Faeth  (1992)  and  Wu  et  al.  (1991, 1992)  for  more  details. 
For  convenience,  the  latter  three  papers  may  be  found  in  Appendices  E,  F  and  G. 


3.2  Experimental  Methods 

Apparatus.  The  apparatus  consisted  of  a  pneumatically-driven  piston/cylinder 
arrangement  as  illustrated  in  Fig.  13.  For  nonturbulent  breakup  experiments  the  cylinder 
outflow  was  through  a  smooth  convergent  passage  with  injection  into  room  air.  For 
turbulent  breakup  experiments,  the  convergent  passage  was  followed  by  a  constant  area 
passage  (41  diameters  long)  to  provide  fully-developed  turbulent  pipe  flow  at  the  exit;  in 
this  case,  injection  was  either  into  room  air  or  into  a  windowed  chamber  filled  with  helium, 
air  or  Freon  12  at  pressures  of  1  and  2  atm.  Liquid  sample  sizes  within  the  cylinder  were 
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Figure  13  Sketch  of  test  apparatus 
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roughly  500  ml  with  passage  diameters  in  the  range  3.6-6.4  mm.  Instrumentation  was 
mounted  rigidly  so  that  various  positions  were  studied  by  traversing  the  entire 
injector/chamber  assembly. 

The  piston/cylinder  arrangement  was  filled  while  venting  gas  trapped  below  the 
piston  face.  Operation  was  initiated  by  admitting  high-pressure  air  to  the  upper  side  of  the 
piston  through  a  solenoid  value  which  forced  the  liquid  through  the  jet  passage.  Total  tunes 
of  injection  were  180-1400  ms  due  to  the  finite  liquid  sample  size;  however,  this  was  not 
problematical  because  flow  development  times  were  generally  less  than  1  %  of  the  flow  time 
for  conditions  of  interest.  Jet  exit  velocities  were  calibrated  using  an  impact  plate. 

Instrumentation.  The  measurements  consisted  of  pulsed  shadowgraphy  and  single- 
and  double-pulse  holography.  The  arrangement  of  the  holocamera  was  the  same  as  Ruff  et 
al.  (1991)  as  illustrated  in  Fig.  14.  An  off-axis  arrangement  was  used  based  on  a  ruby  laser 
with  a  20  ns  pulse  duration  that  could  be  double-pulsed  with  separation  times  as  small  as  1 
jis.  The  optics  allowed  measurements  of  drop  diameters  as  small  as  5  p.m,  while  double 
pulsing  allowed  drop  velocities  to  be  measured.  The  holograms  were  reconstructed  and 
analyzed  using  a  Gould  FD  5000  Image  Processor.  The  shadowgraphs  were  made  by 
blocking  the  reference  beam,  with  other  aspects  of  these  measurements  the  same  as 
holography. 

Test  Conditions.  Test  conditions  for  nonturbulent  primary  breakup  involved  water, 
n- heptane  and  glycerol  mixtures  (21,  42  and  63%  by  mass)  in  still  air  at  atmospheric 
pressure.  Present  work  used  a  5  mm  injector  diameter  with  liquid  velocities  in  the  range 
16.5-133.8  m/s.  This  was  supplemented  by  measurements  using  a  9.5  mm  diameter 
injector,  with  velocities  of  49-60  m/s  and  ambient  air  at  pressures  of  1-4  atm  from  Ruff  et 
al.  (1991,  1992)  and  Tseng  et  al.  (1992a,b). 

Test  conditions  for  turbulent  primary  breakup  involved  water,  n-heptane  and 
glycerol  (42%  dry  mass)  in  still  helium,  air  and  Freon  12  at  pressures  of  1  and  2  atm. 
Injector  diameters  were  3.6  and  6.4  mm  with  liquid  velocities  in  the  range  16-109  m/s. 
These  results  also  were  supplemented  by  measurements  from  Ruff  et  al.  (1991,  1992)  and 
Tseng  et  al.  (1992a,b)  for  a  9.5  mm  diameter  injector. 


3-3  Nonturbulent  Liquid  Breakup 

Similar  to  the  results  discussed  in  Section  2,  drop  sizes  after  primary  breakup 
satisfied  Simmons’  (1977)  universal  root  normal  distribution  with  MMD/SMD  =  1.2. 
Thus,  drop  size  information  will  be  presented  solely  in  terms  of  the  SMD  in  the  following. 

It  was  found  that  drop  characteristics  varied  with  distance  along  the  mixing  layer  at 
first  but  eventually  approached  a  regime  where  the  rate  of  variation  was  relatively  slow. 
This  latter  regime  will  be  termed  the  fully-developed  primary  breakup  regime  in  the 
following.  Approach  to  fully  developed  conditions  was  slower  for  more  viscous  liquids 
and  tended  to  correlate  as  a  function  of  Refx  as  illustrated  for  SMD  in  Fig.  15.  The  data  are 
relatively  limited,  however,  the  fully-developed  primary  breakup  regime  seems  to  be 
reached  for  Refx  on  the  order  of  10^.  The  values  of  ep  for  nonturbulent  primary  breakup 
generally  were  in  the  range  1.2- 1.5  for  the  fully-developed  regime  after  reaching  maximum 
values  in  the  range  1. 5-2.0  for  Refx  near  106  (Wu  et  al.  1991). 

A  correlation  of  the  SMD  after  nonturbulent  primary  breakup  in  the  fully-developed 
regime  was  attempted  based  on  the  existing  theoretical  correlations  of  Taylor  (1940), 
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Levich  (1962),  Mayer  (1961)  and  Adelberg  (1968).  All  these  correlations  can  be  expressed 
in  terms  of  the  same  variables  as  shown  in  Fig.  16.  The  Taylor  (1940)  and  Levich  (1962) 
correlations  are  identical  at  small  values  of  the  independent  variable  (the  Taylor  number) 
and  have  been  Fitted  to  the  measurements  of  Ruff  et  al.  (1991)  for  illustrative  purposes. 
The  limit  of  primary  breakup  proposed  by  Levich  (1962)  at  a  Taylor  number  on  the  order 
of  unity  also  is  shown  on  the  plot.  The  correlations  of  Mayer  (1961)  and  Adelberg  (1968) 
are  plotted  using  recommended  empirical  coefficients  although  these  are  known  to  vary 
widely  for  particular  injector  designs  (Hanje  and  Reardon  1972). 

Clearly,  none  of  the  theoretical  expressions  illustrated  in  Fig.  16  provide  a 
successful  correlation  of  present  measurements.  In  particular,  results  for  a  given  liquid  do 
not  yield  the  correct  velocity  dependence,  while  changes  of  liquids  yield  shifts  of  behavior 
that  are  not  anticipated  by  any  of  the  theories.  The  correlations  of  Mayer  (1961)  and 
Adelberg  (1968)  provide  the  best  approximation  of  velocity  effects,  which  probably 
accounts  for  their  continued  use  for  liquid  rocket  engine  applications  (Harrje  and  Reardon 
1972). 


Initially,  an  alternative  approach  was  developed  based  on  the  development  of 
surface  boundary  layers  on  the  windward  side  of  waves  along  the  liquid  surface  (Wu  et  al. 
1991).  This  was  successful  for  treating  data  from  a  given  injector  but  was  problematical 
due  to  unexplained  shifts  in  measurements  resulting  from  changes  in  the  length  of  the 
injector  passage.  Coupled  with  the  observations  of  small  aerodynamic  effects  from  Hoyt 
and  Taylor  (1985),  this  led  to  a  modified  approach  based  on  vorticity  developed  by  the 
injector  passage  boundary  layers  (Wu  et  al.  1992).  Thus,  present  considerations  will  be 
limited  to  this  mechanism. 

A  sketch  of  the  injector  passage  vorticity  mechanism  appears  in  Fig.  17.  In  this 
case,  it  is  assumed  that  drop  generation  is  due  to  vorticity  produced  by  the  wall  boundary 
layer  formed  when  liquid  passes  through  the  injector  passage.  This  is  plausible  because  the 
surface  only  can  be  distorted  and  ligaments  ejected  when  vorticity  exists  near  the  liquid 
surface.  Thus,  once  the  liquid  leaves  the  injector  passage,  boundary  layer  vorticity  starts  to 
distort  the  liquid  surface  and  causes  ligaments  to  be  shed.  After  a  certain  distance,  liquid 
surface  disturbances  have  fully  developed  and  the  resulting  drop  sizes  remain  almost 
constant  in  the  fully-developed  regime.  It  then  is  assumed  that  drop  sizes  in  the  fully- 
developed  primary  breakup  regime  are  proportional  to  the  wall  boundary  layer  thickness  at 
the  jet  exit.  Due  to  the  large  contraction  ratio,  it  is  further  assumed  that  the  wall  boundary 
layer  remains  laminar.  Then,  applying  predictions  for  a  laminar  flat  plate  boundary  having 
a  nozzle  contraction  curve  length,  Lp,  there  results: 

WegSMD  =  CfiCLp/d^Wegd/Re^  (1) 

Measurements  for  the  fully-developed  nonturbulent  primary  breakup  regime  are 
plotted  in  terms  of  the  variables  of  Eq.  (1)  in  Fig.  18.  Experimental  results  include  present 
measurements  and  those  of  Ruff  et  al.  (1991)  and  Tseng  et  al.  (1992b).  The  best  fit  of  the 
measurements  also  is  shown  on  the  plot;  it  has  a  slightly  different  power  than  Eq.  (1)  as 
follows: 


WegSMD  =  6.8 1  [(Lp/dl^Wegd/Re^]0-85  (2) 

The  reduction  of  the  power  on  the  right  hand  side  of  Eq.  (2)  from  unity  to  0.85  is 
statistically  significant  but  is  not  large  in  view  of  the  qualitative  nature  of  the  development 
of  the  correlation.  The  correlation  coefficient  of  the  fit  is  0.98,  which  is  excellent. 
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Figure  16  Correlation  of  SMD  for  fully-developed  nonturbulent  primary  breakup 
using  existing  theories 


gure  17  Sketch  of  injector  passage  boundary 
layer  vorticity  mechanism 
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Figure  18  Correlation  of  SMD  for  fully-developed  nonturbulent  primary  breakup  based 
on  the  boundary  layer  vorticity  mechanism 


The  excellent  performance  of  Eq.  (2)  suggests  that  boundary  layer  growth  inside 
the  injector  passage  determines  drop  dimensions  in  the  fully-developed  breakup  regime. 
This  implies  that  the  surrounding  gas  plays  no  role  in  the  breakup  process  through 
aerodynamic  effects  even  though  relative  velocities  are  large  for  present  test  conditions, 
with  the  possible  exception  being  a  contribution  to  the  reduction  of  the  power.  This  agrees 
with  the  observations  of  Hoyt  and  Taylor  (1985)  that  there  is  no  significant  change  of 
liquid  surface  structures  in  either  co-flowing  or  counter-flowing  air  streams  flowing  at  a 
variety  of  relative  velocities  at  atmospheric  pressure.  Additionally,  the  results  of  Tseng  et 
al.  (1992)  for  gas  density  variations  over  the  range  1-4  atm  showed  modest  effects  of  gas 
density  on  drop  sizes  after  primary  breakup.  These  results  also  are  seen  to  be  effectively 
correlated  by  Eq.  (2)  in  Fig.  18.  Thus,  it  appears  that  liquid  vorticity  dominates  low  levels 
of  vorticity  due  to  aerodynamic  drag  for  present  test  conditions,  which  accounts  for  the 
failure  of  the  aerodynamic  drag  theories  seen  in  Fig.  16.  This  also  highlights  a  limitation  of 
present  work;  namely,  that  only  two  injectors  of  similar  general  passage  design  were 
considered  (Wu  1991).  Based  on  past  experience  (Harije  and  Reardon  1972)  variations  of 
Cb  with  different  passage  designs  can  be  anticipated.  Finally,  the  connection  between 
these  results  and  gas  injection  processes  for  nonturbulent  liquids  is  not  clear.  However, 
this  is  not  a  major  limitation  because  most  practical  gas  injection  processes  involve 
turbulent  liquids. 

Another  interesting  feature  of  the  results  of  Fig.  18  is  the  potential  for  secondary 
breakup.  First,  of  all,  it  generally  is  agreed  that  drops  are  unstable  to  secondary  breakup 
when  their  Weber  number  exceeds  a  critical  value,  as  follows  (Lefebvre  1980,  Ruff  et  al. 
1992): 


Wegp  —  pg  dp  u02/o  >  6.5  (3) 

Equation  (3)  has  been  written  assuming  that  the  relative  velocity  of  large  drops  subject  to 
secondary  breakup  is  equal  to  the  relative  velocity  of  the  mixing  layer,  which  is  reasonable 
based  on  existing  measurements  (Ruff  et  al.  1991, 1992).  Now  the  universal  root  normal 
drop  size  distribution  implies  that  99.7%  of  the  spray  mass  is  in  the  diameter  range  0.098  < 
dp/SMD  <  3.5.  This  information  and  Eq.  (3)  implies  that  WegSMD  =1-8  and  66  define 
limiting  conditions  where  virtually  no  drops  and  all  drops  undergo  secondary  breakup. 

The  limits  of  secondary  breakup  are  shown  in  Fig.  18  along  with  the  primary 
breakup  correlation.  Noting  that  it  is  difficult  to  achieve  primary  breakup  for  values  of  the 
independent  variable  much  lower  than  the  present  test  range,  it  is  evident  that  drops  formed 
by  nonturbulent  primary  breakup  generally  are  unstable  to  secondary  breakup.  This  helps 
explain  the  durability  of  primary  breakup  correlations  like  those  of  Taylor  (1940)  and 
Levich  (1962):  any  contamination  of  measurements  by  secondary  breakup  causes 
reductions  in  the  range  of  WegsMD»  tending  to  support  the  trends  of  these  relationships, 
since  past  measurements  of  drop  sizes  were  carried  out  after  the  sprays  were  dilute, 
providing  time  for  completion  of  secondary  breakup  processes.  Additionally,  the  results  of 
secondary  breakup  do  not  depend  on  jet  diameter,  which  helps  to  explain  the  weak,  or 
absent,  jet  diameter  effects  seen  in  early  correlations  of  primary  breakup  (Faeth  1990). 

Finally,  it  was  argued  that  effects  of  ambient  gas  density  on  present  observations  of 
nonturbulent  primary  breakup  were  small.  However,  parameters  controlling  secondary 
breakup  (see  Eq.  (3))  respond  rapidly  to  gas  density  changes.  Thus,  higher  gas  densities 
than  those  considered  by  Tseng  et  al.  (1992b)  could  involve  merging  of  primary  and 
secondary  breakup  similar  to  observations  of  turbulent  primary  breakup  to  be  discussed 
next:  this  should  be  explored. 
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3.4 


Turbulent  Liquid  Breakup 

Drop  Size  Distribution.  Similar  to  the  other  conditions,  drop  sizes  after  turbulent 
primary  breakup  satisfied  Simmons'  (1977)  universal  root  normal  distribution  with 
MMD/SMD  =  1.2.  Thus,  drop  sizes  will  be  represented  by  SMD  alone  in  the  following. 

Drop  Velocities.  Drop  velocities  after  turbulent  primary  breakup  were  measured 
over  the  full  range  of  test  conditions.  The  measurements  of  volume  (mass)  averaged 
streamwise  and  crosstream  drop  velocities  are  plotted  as  a  function  of  distance  along  the 
mixing  layer  in  Fig.  19.  Results  for  pf/pg  <  500  are  shown  as  darkened  and  half-darkened 
symbols,  while  results  for  pf/pg  >  500  are  shown  as  open  symbols,  in  order  to  highlight 
properties  at  the  low  liquid/gas  density  ratio  regime  where  aerodynamic  effects  are 
important.  This  practice  will  be  followed  on  subsequent  plots  of  data  as  well. 

The  normalized  drop  velocities  shown  in  Fig.  19  provide  crude  correlations  in 
terms  of  x/d  but  other  factors  are  involved.  Near  the  jet  exit,  up/u0  is  roughly  0.6,  which 
increases  to  roughly  0.9  for  x/d  >  20.  This  behavior  is  consistent  with  streamwise  drop 
velocities  after  turbulent  primary  breakup  being  roughly  the  same  as  streamwise  liquid 
velocities,  with  lower  velocities  near  the  jet  exit  due  to  the  retarding  effect  of  the  passage 
wall.  Superimposed  on  this  behavior,  however,  is  a  clear  trend  that  measurements  for 
pf/pg  <  500  have  significantly  lower  streamwise  velocities  than  the  rest.  This  behavior  can 
be  attributed  to  effects  of  aerodynamic  drag. 

Similar  trends  with  respect  to  density  variations  are  much  less  clear  for  Vp/uQ  in  Fig. 
19.  In  general,  crosstream  drop  velocities  are  comparable  to  crosstream  velocity 
fluctuations  in  the  liquid,  e.g.,  v0'/u0  =  0.058  for  fully-developed  turbulent  pipe  flow 
(Hinze  1975;  Tennekes  and  Lumley  1972).  The  absence  of  an  effect  of  aerodynamic  drag 
in  this  case  can  be  explained  due  to  a  compensating  aerodynamic  acceleration  effect  over  the 
tips  of  liquid  protuberances  that  will  be  discussed  next  in  connection  with  the  onset  of 
breakup. 

Onset.  Phenomenological  analysis  to  find  drop  properties  at  the  onset  of  breakup  is 
based  on  the  concept  illustrated  in  Fig.  20.  It  is  assumed  that  onset  involves  the  formation 
of  a  drop  from  a  turbulent  eddy  of  characteristic  size  £  and  a  characteristic  crosstream 
velocity  relative  to  the  surrounding  liquid,  vjj.  The  smallest  drops  that  can  be  formed  either 
are  comparable  to  the  Kolmogorov  scale  of  the  turbulence  or  the  smallest  eddy  that  has 
sufficient  mechanical  energy  needed  to  form  a  drop,  whichever  is  larger.  For  present 
conditions,  drop  energy  was  the  limiting  factor  and  only  this  limit  was  pursued.  The 
mechanical  energy  available  consists  of  the  kinetic  energy  of  the  eddy  relative  to  its 
surroundings,  plus  the  added  mechanical  energy  due  to  the  pressure  drop  caused  by 
acceleration  of  the  surrounding  gas  over  the  tip  of  the  protuberance.  Then  equating  these 
energies,  assuming  that  the  eddy  is  in  the  inertial  range  of  the  turbulence  spectrum,  to  the 
surface  tension  energy  of  the  resulting  drop,  there  results: 
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Figure  20  Sketch  of  aerodynamical ly-enhanced  turbulent  primary 
breakup  at  the  liquid  surface 


(SMDj/A)(l  +  Csa(pg/pf)(u0^o,)2(A/SMDi)2^)3/5 


=  Csid^oO^We'f  (4) 

where  turbulence  properties  in  the  liquid  are  assumed  to  be  the  same  as  at  jet  exit 
conditions,  C^  and  CSi  are  empirical  constants  and  A  is  the  radial  integral  scale  (which  is 
taken  to  be  equal  to  d/8  for  fully-developed  turbulent  pipe  flow,  see  Hinze  (1975)). 

Present  measurements  of  SMDj  are  plotted  in  terms  of  Eq.  (4)  in  Fig.  21.  In  doing 
this,  was  optimized  to  a  value  of  0.04,  based  on  taking  v07u0  =  0.058  for  fully- 
developed  turbulent  pipe  flow  (Hinze  1975).  The  correlation  of  the  data  all  conditions 
generally  is  within  anticipated  experimental  uncertainties.  The  power  of  WefA  is  not  -3/5  as 
suggested  by  Eq.  (4),  however,  and  the  results  can  be  correlated  best  by  the  following 
empirical  fit  that  is  shown  on  the  plot: 

(SMDj/A)(l  +  0.04(pg/pf)(uoA?o,)2(A/SMDi)2/3)3/5 

=  2.5(Uo£o')6/5We;°-69  (5) 

The  reduction  in  the  power  of  WefA  from  Eq.  (4)  to  Eq.  (5)  is  not  statistically  significant 
while  the  correlation  coefficient  of  the  fit  0.97,  which  is  excellent.  Aerodynamic  effects  are 
carried  by  the  term  multiplied  by  Csa  =  0.04  which  largely  depends  on  the  liquid/gas 
density  ratio  and  is  small  for  pf/pg  >  500. 

Given  SMDj,  it  was  assumed  that  the  drop-forming  eddy  convects  along  the  surface 
with  a  velocity  uQ  for  the  distance  xj  required  to  form  a  drop.  The  time  of  convection  was 
taken  to  be  the  Rayleigh  breakup  time  of  the  protruding  ligament  found  from  Weber 
(1931),  after  ignoring  a  liquid  viscosity  term  that  is  small  for  present  test  conditions.  Then 
using  Eq.  (4)  to  specify  SMDj,  the  following  expression  for  xj  is  obtained: 

(Xi/A)(l  +  Csa(pg/pf)(^o’)2(A/SMDi)2^)9/10 

=  Cxi(^o')9/5We^/10  (6) 


where  Cxi  is  a  constant  of  proportionality.  Similar  to  Eqs.  (4)  and  (5),  aerodynamic  effects 
in  Eq.  (6)  largely  are  controlled  by  the  liquid/gas  density  ratio. 

Adopting  Csa  -  0.04  as  before,  the  measurements  of  xj  are  plotted  in  terms  of  the 
variables  of  Eq.  (6)  in  Fig.  22.  The  correlation  of  the  measurements  is  reasonably  good  in 
view  of  the  relatively  large  experimental  uncertainties  of  xj.  The  best  fit  of  the 
measurements  can  be  expressed  as  follows: 
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Figure  21  SMD  at  the  onset  of  turbulent  primary  breakup 


Figure  22  Length  to  initiate  turbulent  primary  breakup 
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(xj/A)(l  +  0.04(pg/pf)(u<^o,)2(A/SMDi)2^)9/10 


=  12(uo/^o')9/5  We£63  (7) 

with  the  correlation  coefficient  of  the  fit  being  0.89.  Thus,  the  difference  between  the 
powers  of  WefA  in  Eqs.  (6)  and  (7)  is  statistically  significant  but  is  not  large  in  view  of  the 
approximations  used  to  develop  the  correlating  expression. 

Drop  Sizes  after  Onset.  Aerodynamic  secondary  breakup  times  scale  like 
JUpf/Pg)1/2u0  for  an  object  of  size  fi  and  relative  velocity  u0  (Ranger  and  Nicholls  1969; 
Hsiang  and  Faeth  1992,  1993).  Thus,  for  conditions  where  aerodynamic  effects  are 
important,  the  Rayleigh  breakup  times  of  ligaments  increase  more  rapidly  than  secondary 
breakup  times  as  A  increases.  This  yields  two  possible  turbulent  primary  breakup 
mechanisms  after  the  onset  conditions:  (1)  nonaerodynamic  primary  breakup,  and  (2) 
merged  primary  and  secondary  breakup  when  aerodynamic  effects  are  important.  The 
configuration  of  the  two  mechanisms  is  illustrated  in  Fig.  23. 

When  aerodynamic  effects  are  not  important,  it  is  assumed  that  the  SMD  is 
proportional  to  the  largest  drop  that  can  be  formed  at  a  particular  position  because  the  SMD 
is  dominated  by  the  largest  drops.  Then  adopting  the  Rayleigh  breakup  mechanism  to 
relate  drop  size  and  position,  similar  to  the  treatment  of  onset  conditions,  the  following 
expression  is  obtained: 


SMD/A  =  Csx(x/(A  Wejf))2^  (8) 

where  Csx  is  an  empirical  constant. 

Under  the  present  assumptions,  principles  relating  SMD  to  x  are  the  same  for  both 
the  onset  of  breakup  and  at  larger  distances,  which  suggests  that  SMD/A  should  be 
correlated  as  a  function  of  x/(A  We^'54)  from  Eqs.  (5)  and  (7).  Present  measurements  and 
those  of  Ruff  et  al,(1991)  for  pf/pg  >  500,  where  aerodynamic  effects  are  not  important, 
are  correlated  in  terms  of  these  variables  in  Fig.  24.  Two  coiTelations  are  shown  on  the 
plot  as  well:  (1)  the  best  fit  of  all  the  data,  and  (2)  the  correlation  based  on  onset  conditions 
alone.  The  best  fit  correlation  of  all  the  data  is: 

SMD/A  =  0.69[  x/(A  Wt^'54  ]0-57  (9) 

which  yields  a  correlation  coefficient  of  0.98,  which  is  excellent.  Notably,  results  for 
injection  into  still  air  and  helium  at  atmospheric  pressure  are  essentially  the  same,  in  spite  of 
the  substantial  difference  of  pf/pg,  proving  that  aerodynamic  effects  are  small  for  these 
conditions. 

Differences  between  the  two  correlations  illustrated  in  Fig.  24  are  small  so  that  drop 
formation  both  at  the  onset  of  breakup  and  at  distances  farther  along  the  multiphase  mixing 
layer  appear  to  be  similar  when  aerodynamic  effects  are  weak.  The  SMD  of  the  smallest 
drops,  formed  at  the  onset  of  breakup,  are  in  the  range  10-20  pm  which  is  roughly  ten 
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Figure  23  Sketch  of  turbulent  breakup  processes 
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Figure  24  SMD  after  turbulent  primary  breakup  when  aerodynamic  effec 


times  larger  than  Kolmogorov  length  scales  for  these  test  conditions.  The  characteristics  of 
turbulent  breakup  should  change  when  Kolmogorov  scales  are  involved  which  is  an 
interesting  issue  for  future  study.  The  largest  SMD  seen  in  Fig.  24  approach  the  radial 
integral  scales  of  the  turbulence;  much  larger  values  are  unlikely  for  round  liquid  jets  due  to 
the  finite  length  of  the  liquid  core.  This  is  shown  on  the  figure  by  indicating  the  core  length 
from  Grant  and  Middleman  (1966): 

Lc/d  =  8.5  We°'32  (10) 

where  the  range  of  Lc/d  shown  in  the  figure  is  due  to  the  variation  of  Wefd  over  the  present 
test  range.  A  more  recent  correlation  of  Chehroudi  et  al.  (1985)  yields  similar  values  of  Lc 
but  a  somewhat  broader  range.  That  the  SMD  approaches  A  when  x  approaches  Lc,  and 
thus  is  comparable  to  the  diameter  of  the  liquid  core  itself  in  this  region,  is  consistent  with 
the  liquid  column  breaking  up  as  a  whole. 

When  aerodynamic  effects  are  important,  it  is  assumed  that  turbulent  primary  and 
secondary  breakup  merge  as  illustrated  in  Fig.  23.  In  this  regime,  Eq.  (8)  is  used  to 
estimate  the  drop  size  after  primary  breakup  with  properties  after  secondary  breakup  found 
from  the  recent  results  of  Hsiang  and  Faeth  (1992,  1993).  Assuming  that  the  appropriate 
relative  velocity  is  uQ,  the  expression  for  the  SMD  after  secondary  breakup  of  a  drop 
having  a  diameter  dp  is  as  follows  (Hsiang  and  Faeth  1992): 

SMD/dp  =  Cs(pf/pg)l/2(W/(pfdpIi0))l/2  (11) 

Then  assuming  that  dp  is  proportional  to  the  SMD  after  primary  breakup,  Eqs.  (8)  and  (11) 
can  be  combined  to  yield  the  SMD  after  merged  primary  and  secondary  breakup  as  follows: 

pgSMDu02/a  =  CsCs1f(x/A)l/3(pg/pf)3/4We5/6Re-ii/2  (12) 

where  Cs  in  Eq.  (12)  is  modified  from  Eq.  (11)  to  account  for  a  range  of  initial  drop  sizes 
rather  than  a  single  drop. 

Present  measurements  of  aerodynamic  turbulent  primary  breakup,  along  with  some 
results  from  Tseng  et  al.  (1992b),  are  plotted  according  to  the  variables  of  Eq.  (12)  in  Fig. 
25.  The  following  best  fit  correlation  of  the  measurements  also  is  shown  on  the  plot: 

pgSMDu02/c  =  12[(x/A)l/3(pg/pf)3/4we^6Re^/211-17  (13) 

which  yields  an  overall  correlation  coefficient  of  the  fit  of  0.92.  Additionally,  adopting  the 
best  fit  values  of  Cs  and  Csx  in  Eq.  (12)  yields  the  following  theoretical  prediction  for  the 
SMD: 


pgSMDu02/G  =  12.9(x/A)l/3(pg/pf)3/2We^6Re^/2  (14) 

which  also  is  shown  on  the  plot.  The  standard  deviation  of  the  constant  in  Eq.  (14)  is 
21%;  thus,  the  differences  between  Eqs.  (13)  and  (14)  are  not  statistically  significant.  All 
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the  measurements  are  in  reasonably  good  agreement  with  Eqs.  (13)  and  (14)  except  for  the 
highest  pressure  results  of  Tseng  et  al.  (1992b)  which  yield  significantly  smaller  drops  near 
the  liquid  surface.  An  explanation  of  this  behavior  is  that  these  flows  involved  the  largest 
drop  concentrations  of  all  the  measurements  so  that  drop  properties  near  the  surface  could 
be  influenced  by  small  drops  migrating  from  other  parts  of  the  multiphase  mixing  layer  due 
to  turbulent  dispersion. 

Breakup  Regimes.  The  present  measurements  suggest  three  regimes  of  turbulent 
primary  breakup:  (1)  non-aerodynamic  turbulent  primary  breakup;  (2)  aerodynamically- 
enhanced  turbulent  primary  breakup,  observed  at  onset  conditions;  and  (3)  aerodynamic 
turbulent  primary  breakup,  involving  merged  turbulent  primary  and  secondary  breakup. 
The  results  also  indicated  that  the  liquid/gas  density  ratio  and  the  relative  characteristic 
Rayleigh  breakup  times  of  ligaments  and  secondary  breakup  times  of  liquid  objects  fixed 
the  boundaries  of  these  regimes.  The  ratio  of  the  characteristic  times  can  be  expressed  as 
follows  (Wu  and  Faeth  1992): 

tR/tb  =  (pg/pf)1/2(xWefA/A)l/3  (15) 

The  resulting  primary  breakup  regimes,  based  on  present  measurements  and  those 
of  Ruff  et  al.  (1991)  and  Tseng  et  al.  (1992b),  are  illustrated  in  terms  of  pf/pg  and  tR/tb  in 
Fig.  26.  The  results  yield  the  following  criteria  for  transitions  between  regimes 

pf/pg  =  500,  non-aerodynamic/aerodynamic  transition  (16) 

tR/tb  =  4,  enhanced/merged-aerodynamic  transition  (17) 

Thus,  Eqs.  (4)  and  (5)  always  provide  onset  conditions  —  whether  aerodynamic  effects  are 
important  or  not.  Equation  (9)  provides  the  SMD  subsequently,  when  aerodynamic  effects 
are  not  important.  Finally,  Eq.  (12)  provides  the  SMD  after  onset,  when  aerodynamic 
effects  are  important  in  the  merged  primary  and  secondary  breakup  regime. 

3.5  Conclusions 

The  mechanisms  of  primary  breakup  of  nonturbulent  and  turbulent  liquids  were 
studied  within  multiphase  mixing  layers.  Various  liquids  and  ambient  gas  conditions  were 
studied  to  highlight  effects  of  both  liquid  properties  and  liquid/gas  density  ratio  —  the  latter 
controlling  the  importance  of  aerodynamic  phenomena.  Confirming  earlier  observations, 
all  breakup  processes  yielded  drop  size  distributions  that  satisfied  Simmons'  (1977) 
universal  root  normal  distribution  with  MMD/SMD  =  1.2.  Additionally,  drops  formed  by 
primary  breakup  are  intrinsically  unstable  to  secondary  breakup,  which  is  a  dominate 
process  in  the  multiphase  mixing  layer.  Other  conclusions  differ  for  nonturbulent  and 
turbulent  primary  breakup  as  discussed  in  the  following. 

Major  conclusions  concerning  nonturbulent  primary  breakup  are  as  follows: 

1 .  Drop  properties  after  primary  breakup  varied  with  distance  along  the  mixing  layer 
but  became  relatively  uniform  in  a  fully-developed  primary  breakup  regime  for  Refx 
greater  than  106,  see  Fig.  15. 
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Figure  26  Turbulent  primary  breakup  regime  map 
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2 .  Existing  correlations  of  SMD  after  primary  breakup  based  on  aerodynamic  wave 
growth,  due  to  Taylor  (1940),  Levich  (1962),  Mayer  (1961)  and  Adelberg  (1968), 
were  ineffective  for  present  measurements  because  aerodynamic  effects  were  small 
and  were  dominated  by  low  levels  of  liquid  vorticity  at  the  start  of  the  multiphase 
mixing  layer,  see  Fig.  16. 

3.  A  reasonably  good  correlation  of  SMD  after  fully-developed  primary  breakup  — 
for  various  liquids,  injector  diameters,  relative  velocities  and  liquid/gas  density 
ratios — was  achieved  using  Eq.  (2),  see  Fig.  18.  This  expression  was  based  on 
drop  sizes  being  proportional  to  the  thicknesses  of  boundary  layers  in  the  injector 
passage,  with  the  best  empirical  fit  fortuitously  accounting  for  aerodynamic  effects 
as  well  as  over  the  present  test  range. 


Major  conclusions  concerning  turbulent  primary  breakup  are  as  follows: 

1 .  The  presence  of  aerodynamic  phenomena  for  primary  breakup  is  largely  controlled 
by  the  liquid/gas  density  ratio,  and  appear  when  this  ratio  is  less  than  500,  see  Fig. 


2.  Drop  sizes  at  the  onset  of  turbulent  primary  breakup  could  be  explained  by  equating 
the  surface  tension  energy  required  to  form  a  drop  to  both  the  kinetic  energy  of  a 
corresponding  liquid  eddy  relative  to  its  surroundings,  assisted  by  the  pressure 
reduction  over  the  tips  of  liquid  protrusions  when  aerodynamic  effects  are 
important.  Phenomenological  analysis  based  on  these  ideas  yielded  a  good 
correlation  of  drop  sizes  at  the  onset  of  breakup,  see  Fig.  21,  as  well  as  the  distance 
along  the  mixing  layer  required  for  breakup  to  begin,  see  Fig.  22. 

3 .  Drop  sizes  from  primary  breakup  progressively  increased  with  distance  along  the 
mixing  layer,  which  could  be  explained  by  associating  the  SMD  with  the  largest 
drops  that  had  sufficient  residence  time  to  form  by  Rayleigh  breakup  of  protruding 
liquid  elements  when  aerodynamic  effects  were  small,  see  Fig.  24. 

4.  For  conditions  where  secondary  breakup  times  become  small  in  comparison  to  the 
Rayleigh  breakup  time  of  ligaments,  processes  of  primary  and  secondary  breakup 
merge  to  reduce  drop  sizes.  The  reduction  of  drop  sizes  could  be  explained  by 
combining  the  turbulent  primary  breakup  mechanism  with  secondary  breakup 
results  of  Hsiang  and  Faeth  (1992,  1993),  see  Fig.  25. 

5.  Mass-averaged  drop  velocities  after  primary  breakup  approximate  mean  and  rms 
velocity  fluctuations  in  the  liquid  in  the  streamwise  and  crosstream  directions, 
respectively.  However,  there  was  a  tendency  for  streamwise  velocities  to  be 
reduced  somewhat  when  aerodynamic  effects  were  important,  see  Fig.  19. 
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The  structure  and  breakup  properties  of  the  dense-spray  region  of  nonevaporating  pressure-atomized  sprays 
were  studied.  The  multiphase  mixing  layer  near  the  injector  exit  was  emphasized,  considering  large-scale  |9.J- 
mm  injector  diameter)  water  jets  injected  vertically  downward  in  still  room  air.  Phase-discriminating  laser 
vetocimetry  and  double- pulse  holography  were  used  to  measure  phase  velocities  and  drop-size  properties  for 
both  nonturbulent  and  turbulent  jet  exit  conditions.  Present  test  conditions  involved  two  types  of  primary 
breakup:  1)  aerodynamic  breakup  for  nonturbulent  jets  where  properties  could  be  correlated  using  earlier 
aerodynamic  breakup  theories;  and  2)  turbulent  breakup  for  turbulent  jets  where  drop  properties  could  be 
related  to  liquid  turbulence  properties.  Both  mechanisms  yielded  Weber  numbers  exceeding  secondary  drop 
breakup  limits  near  the  liquid  surface.  Significant  effects  of  separated  flow  were  observed  for  present  test 
conditions;  however,  scaling  analysis  suggests  reduced  effects  of  separated  flow  at  higher  Injector  velocities  and 
ambient  pressures — largely  due  to  finer  atomization. 


Nomenclature 

C„  =  aerodynamic  breakup  coefficient 

C„  =  drop  breakup  time  coefficient 

Cc  =  liquid  core  length  coefficient 

CD  -  drop  drag  coefficient 

C„  C;  =  turbulent  breakup  coefficients 

d  =  injector  exit  diameter 

d„  =  drop  diameter 

k  -  Favre-averaged  turbulence  kinetic  energy 

L  =  injector  passage  length 

Lc  =  liquid  core  length 

N  =  number  of  measurements 

Oh  =  jet  Ohnesorge  number,  ji/(p,dcr)1'2 

Re  =  jet  Reynolds  number,  p/u0d/pr 

r  =  radial  d;stance 

5  =•  separated-flow  factor,  Eq.  (9) 

r  =  time 

u  =  streamwise  velocity 

u,  =  relative  streamwise  velocity 

V  =  volume  of  sample 

v  =  radial  velocity 

We,  =  jet  Weber  number  based  on  phase  i,  p,u20dltr 
We„,  =  drop  Weber  number  number  based  on  phase  i, 
P,u}d,l<r 

We,  =  Favre-averaged  spray  Weber  number,  Eq.  (8) 
x  =  streamwise  distance 

A  =  radial  integral  scale 

M  =  molecular  viscosity 

p  =  density 

a  =  surface  tension 

=  drop  breakup  time,  Eq.  (18) 
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re  =  liquid  core  residence  time,  Eq.  (12) 

Tf  =  flow  residence  time,  Eq.  (10) 

Tp  =  drop  response  residence  time,  Eq.  (16) 


Subscripts 

f 

g 

P 

s 

0 


liquid-phase  property 
gas-phase  property 
drop  property 
liquid  surface  property 
injector  exit  condition 


Superscripts 

C),  (')’  =  time-averaged  mean  and  root-mean-squared 
fluctuating  quantities 

('),  (")'  =  Favre-averaged  mean  and  root-mean-squared 
fluctuating  quantities 
(  )*  =  critical  breakup  condition 


Introduction 

THERE  have  been  numerous  efforts  to  develop  methods 
to  analyze  spray  processes.  The  present  investigation  seeks 
to  contribute  to  this  methodology  by  studying  the  dense-spray 
region  near  the  exit  of  the  injector  passage  for  pressure  atom¬ 
ization.  Experiments  were  limited  to  relatively  large-scale  (9.5- 
mm  injector  diameter)  nonevaporating  round  water  jets  in¬ 
jected  into  still  air  at  normal  temperatures  and  pressures. 
Earlier  measurements  of  liquid  volume  fractions,  entrainment 
rates,  and  dispersed-phase  properties  in  the  dense-spray  re¬ 
gion  of  these  flows12  were  extended  to  provide  gas-phase  and 
additional  dispersed-phase  properties.  The  new  measure¬ 
ments,  along  with  additional  information  obtained  from  the 
data  of  Ruff  et  al.,1  were  used  to  continue  study  of  locally 
homogeneous-flow  (LHF)  analysis  of  the  process,  e  g.,  anal¬ 
ysis  based  on  the  assumption  of  infinitely  fast  interphase  trans¬ 
port  rates  so  that  both  phases  have  the  same  velocity  and  are 
in  thermodynamic  equilibrium  at  each  instant  and  point  within 
the  flow.  Present  measurements  were  limited  to  the  atomi¬ 
zation  breakup  regime  where  a  multiphase  mixing  layer  along 
the  edge  of  the  flow  begins  to  develop  right  at  the  injector 
exit.3  4  Jet  exit  conditions  involved  both  slug  and  fully  devel¬ 
oped  turbulent  pipe  flow  because  past  work  had  shown  sen¬ 
sitivity  of  dense-spray  properties  to  liquid-phase  turbulence 
levels.1 2 

For  atomization  breakup,  the  flow  near  the  injector  exit 
involves  a  liquid  core  (much  like  the  potential  core  of  a  single- 
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h  e  jet),  surrounded  by  a  multiphase  mixing  layer  that  be- 
onTto  develop  right  at  the  injector  exit.3"5  The  dense-spray 

enon  is  normally  considered  to  include  both  the  liquid  core 
and  the  multiphase  mixing  layer  up  to  the  point  where  the 
liquid  core  disappears.  There  have  been  several  studies  of  the 
length  of  the  liquid  core,  taken  as  the  length  of  unbroken 
liquid  extending  from  the  passage  exit.6-1*  Findings  indicate 
that  this  length  is  influenced  by  the  breakup  regime,  turbu¬ 
lence  properties  at  the  injector  exit,  and  the  gas/liquid  density 
ratio.  For  atomization  breakup  of  typical  liquids  in  gases  at 
atmospheric  pressure,  however,  the  dense-spray  region  ex¬ 
tends  quite  far  from  the  injector,  ca.  200-400  injector  di¬ 
ameters.  Thus,  dense  sprays  are  an  important  feature  of  spray 
injection  processes  due  both  to  their  extent  and  their  influence 
on  drop  properties  at  the  start  of  the  dilute-spray  region. 

The  multiphase  mixing  layer  is  also  affected  by  the  breakup 
regime,  jet  exit  turbulence  properties,  and  the  density  ratio 
of  the  flow.‘-2-,~11  Ruff  et  al.‘  provide  information  on  liquid- 
phase  properties  in  the  multiphase  mixing  layer  for  water  jets 
in  still  air  at  atmospheric  pressure,  using  double-pulse  hol¬ 
ography  to  find  liquid  element  (drop)  sizes  and  velocities.  It 
was  found  that  the  multiphase  mixing  layer  had  liquid  volume 
fractions  less  than  1%  for  their  test  conditions.  The  inner 
portion  of  the  mixing  layer  contained  large  irregularly  shaped 
liquid  elements  and  drops  while  the  proportion  of  spherical 
drops  increased  and  drop  sizes  decreased  with  increasing  ra¬ 
dial  distance,  suggesting  significant  effects  of  secondary  breakup 
in  the  flow.  The  velocities  of  large  drops  were  generally  much 
larger  than  small  drops,  implying  that  separated  flow  effects 
were  important  as  well.  Increased  turbulence  levels  at  the  jet 
exit  had  a  substantial  effect  on  the  structure  of  the  flow, 
increasing  the  number  and  size  of  irregular  liquid  elements 
by  promoting  ejection  of  liquid  from  the  liquid  surface  and 
increasing  the  width  of  the  liquid-containing  region  of  the 
multiphase  mixing  layer.  Wu  et  al.*  also  observed  significant 
effects  of  jet  exit  conditions  on  the  rate  of  spread  of  the  outer 
edge  of  the  multiphase  mixing  layer. 

The  complexities  of  dense  sprays — involving  stripping  of 
liquid  from  the  all-liquid  core,  the  presence  of  irregular  liquid 
elements,  secondary  breakup,  and  turbulent  dispersion  of  drops 
toward  the  edge  of  the  flow — make  the  LHF  approximation 
attractive  as  a  means  of  circumventing  detailed  descriptions 
of  these  phenomena.  Ruff  et  al.2  found  that  LHF  predictions 
were  reasonably  effective  for  estimating  distributions  of  liquid 
volume  fractions  for  atomization  breakup  in  the  region  where 
mean  liquid  volume  fractions  were  greater  than  0.2.  However, 
the  multiphase  mixing  layer  dominates  flow  properties  at  lower 
liquid  volume  fractions  and  exhibits  significant  effects  of  sep¬ 
arated  flow,  which  limits  the  effectiveness  of  the  LHF  ap¬ 
proximation.1  Other  evaluations  of  the  LHF  approximation 
find  varying  degrees  of  success,  at  times  yielding  encouraging 
results  and  at  other  times  overestimating  the  rate  of  devel¬ 
opment  of  the  flow  as  the  spray  becomes  dilute.”- 12  Unfor¬ 
tunately,  methods  for  determining  conditions  when  the  LHF 
approximation  is  appropriate  for  dense  sprays  have  not  been 
developed.  The  main  difficulty  is  that  success  of  the  LHF 
approximation  depends  on  drops  being  small  enough  to  re¬ 
spond  quickly  to  changes  in  gas  properties,  while  information 
on  both  drop  sizes  and  gas  properties  in  dense  sprays  is  very 
limited. 

The  objective  of  the  present  investigation  was  to  seek  a 
better  understanding  of  the  properties  of  dense  sprays  by 
studying  both  gas-  and  dispersed-phase  velocities  in  the  mul¬ 
tiphase  mixing  layer  for  the  same  test  conditions  as  Ruff  et 
al.1 2  Coupled  with  the  existing  measurements  of  liquid-phase 
properties  (drop  sizes  and  velocities),  this  provides  a  relatively 
complete  picture  of  the  structure  of  these  flows.  Furthermore, 
knowledge  of  both  drop  size  distributions  and  phase  velocities 
allows  direct  evaluation  of  primary  breakup  properties,  the 
propensity  for  secondary  drop  breakup,  and  quantitative  es¬ 
timates  of  effects  of  separated  flow.  Mean  and  fluctuating  gas 
velocities  near  the  edge  of  the  multiphase  mixing  layer  were 


measured  using  phase-discriminating  laser  veiocimetry  (LV). 
whereas  double-pulse  holography  was  used  near  the  liquid 
surface  where  LV  was  no  longer  feasible.  The  measurements 
were  compared  with  predictions  based  on  the  LHF  approxi¬ 
mation  to  help  provide  a  measure  of  separated-flow  effects. 
Scaling  considerations  of  drop  breakup  and  response  were 
also  used  to  interpret  dense-spray  properties  and  gam  insight 
concerning  use  of  the  LHF  approximation. 

The  paper  begins  with  brief  descriptions  of  experimental 
methods  and  the  approach  used  for  the  LHF  computations. 
Experimental  results  are  then  described,  considering  phase 
velocities,  liquid  breakup  properties,  and  separated-flow  pa¬ 
rameters.  The  paper  concludes  by  considenng  scaling  of  var¬ 
ious  processes  in  the  flow.  The  present  discussion  is  brief: 
additional  details  and  a  complete  tabulation  of  data  can  be 
found  in  Ruff.13 

Experimental  Methods 

Apparatus 

The  experimental  apparatus  was  identical  to  past  work.1 2 
and  will  be  described  only  briefly.  The  arrangement  involved 
large-scale  (9.5-mm  injector  diameter)  water  jets  injected  ver¬ 
tically  downward  in  still  room  air.  City  water  was  supplied  to 
the  injector  by  a  centrifugal  pump  with  the  water  flow  rate 
measured  using  a  paddle-wheel  flow  meter  that  was  calibrated 
by  collecting  water  for  timed  intervals. 

The  slug  and  fully  developed  flow  injectors  were  also  iden¬ 
tical  to  Ruff  et  al..1-2  and  are  described  by  Ruff.13  Measure¬ 
ments  of  mean  and  fluctuating  velocities  across  the  exit  of  the 
injectors13  showed  ihat  the  slug  flow  injector  provided  rela¬ 
tively  uniform  velocities,  with  a  streamwise  turbulence  inten¬ 
sity  of  roughly  1%.  at  the  jet  exit,  while  the  fully  developed 
flow  injector  yielded  jet  exit  properties  that  approximated 
fully  developed  turbulent  pipe  flow.1*  15  The  tests  involved  a 
mean  jet  velocity  of  56.3  m/s,  yielding  other  jet  exit  param¬ 
eters  as  follows:  water  flow  rate.  3.99  kg/s;  Reynolds  number, 
534,000;  Weber  numbers  based  on  gas  and  liquid  densities. 
500  and  412 .000,  respectively;  and  Ohnesorge  number,  0.0012 1 . 
These  conditions  are  well  within  the  atomization  breakup 
regime  and  properly  yielded  a  multiphase  mixing  layer  that 
began  at  the  exit  of  the  injector.3* 

Instrumentation 

Laser  Veiocimetry 

The  ambient  air  was  seeded  with  condensed  oil  particles 
having  diameters  less  than  1  /zm  for  LV  measurements.  The 
seeding  particles  had  a  flat  frequency  response  to  the  gas 
motion  up  to  about  30  kHz,  which  was  adequate  for  present 
LV  measurements.  A  large  enclosure  (3  x  3  x  4  m  high) 
was  constructed  around  the  spray  facility  to  avoid  contami¬ 
nating  laboratory  equipment  with  seeding  particles.  For  an 
enclosure  of  this  size,  the  induced  counterflow  velocity  is  less 
than  2%  of  the  axial  velocities  in  the  region  measured  near 
the  edge  of  the  jets. 

A  phase-discriminating  LV  system,  similar  to  that  of  Mo- 
darress  et  al.,1&  was  used  to  avoid  biasing  gas-phase  velocity 
signals  with  liquid  velocity  signals.  The  LV  signal  was  obtained 
using  the  green  line  (514.5  nm)  of  an  argon-ion  laser  (4W, 
Coherent.  INNOVA  90-4)  in  the  dual-beam  forward-scatter 
mode.  A  3.75:1  beam  expander  was  used  to  minimize  the 
dimensions  of  the  measuring  volume  and  to  improve  signal- 
to-noise  ratios.  The  plane  of  the  LV  beams  was  rotated  to 
measure  both  streamwise  and  crosstream  velocities.  The  LV 
was  frequency  shifted  (40  MHz  Bragg  cell,  TSI  model  9180- 
12)  to  eliminate  effects  of  directional  bias  and  ambiguity.  The 
receiving  optics  observed  the  probe  volume  of  the  LV  at  an 
angle  of  30  deg  from  the  forward-scattering  direction,  yielding 
a  measuring  volume  having  a  diameter  of  60  /im  and  a  length 
of  110  #xm. 
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The  phase-discrimination  system  involved  surrounding  the 
LV  measuring  volume  with  the  beam  from  a  5-mW  HeNe 
laser  directed  at  an  angle  of  15  deg  from  the  LV  axis,  with 
collection  optics  in  the  forward-scattering  direction  also  at  an 
angle  of  15  deg  from  the  LV  axis.  This  yielded  a  region  viewed 
by  the  discriminator  that  had  a  diameter  of  0.6  mm  and  a 
length  of  1.3  mm.  surrounding  the  LV  measuring  volume. 
Thus,  droplets  that  could  graze  or  cross  the  LV  measuring 
volume  yielded  a  scattering  signal  on  the  discriminator  output 
so  their  velocities  could  be  eliminated  from  the  velocity 
record. 

The  LV  scattering  signal  was  recorded  using  a  photomul¬ 
tiplier  (TSI  rr  del  9160)  and  processed  using  a  burst  counter 
(TSI  model  1990  C).  The  measurements  involved  low  burst 
densities  (one  seeding  particle  in  the  measuring  volume)  and 
high  data  densities  (time  between  validated  velocity  signals 
small  in  comparison  to  integral  time  scales);  therefore,  the 
analog  output  of  the  processor  was  time  averaged  to  yield 
unbiased  time  averages.  This  involved  low-pass  filtering  (Ith- 
aco  model  4213)  of  the  output  signal  before  it  was  digitized 
(LeCroy  models  8212A/8  and  8800A)  and  transferred  to  a 
microcomputer  for  processing  and  storage. 

The  performance  of  the  phase-discriminator  system  was 
checked  by  measuring  LV  data  rates  with  and  without  seeding 
particles  present.  It  was  required  that  the  LV  data  rate  at 
least  double  for  the  measurements  to  be  considered  valid. 
Additionally,  no  measurement  was  accepted  if  more  than  40 % 
of  the  LV  signal  was  rejected  due  to  signals  from  the  phase 
discriminator  system.  Evaluation  of  this  approach,  using  ve¬ 
locities  of  small  particles  measured  by  double-pulse  holog¬ 
raphy  as  an  upper  bound  for  gas  velocities,  indicated  signif¬ 
icant  upward  bias  of  LV  gas  velocity  measurements  by  particles 
if  these  limits  were  exceeded. 

LV  signals  were  averaged  for  2  min  (128.000  samples)  to 
provide  repeatable  values  of  mean  and  fluctuating  gas  veloc¬ 
ities.  Experimental  uncertainties  (95%  confidence)  are  esti¬ 
mated  to  be  less  than  8%  for  mean  streamwise  velocities,  and 
less  than  15%  for  streamwise  and  cross-stream  velocity  fluc¬ 
tuations,  largely  dominated  by  finite  sampling  times.  Mea¬ 
surements  were  repeatable  well  within  these  limits. 

Double-Pulse  Holography 

LV  data  rates  from  drops  alone  were  too  high  near  the 
liquid  surface  so  that  double-pulse  holography  was  used  in 
this  region.  This  approach  is  tedious  and  less  accurate,  due 
to  difficulties  in  accumulating  sufficient  velocity  samples  and 
reduced  particle  response  because  only  particles  having  di¬ 
ameters  of  roughly  5  pm  and  greater  could  be  measured. 
Nevertheless,  holography  provided  access  to  portions  of  the 
multiphase  mixing  layer  that  were  no  longer  feasible  using 
LV. 

Ruff  et  al. 1  describe  the  holocamera  and  reconstruction 
systems.  Velocity  data  were  obtained  over  6x6x4  mm 
volumes,  using  at  least  three  holograms  per  position.  The  data 
were  spatially  averaged  over  the  width  of  the  measuring  vol¬ 
umes,  or  *  one-half  the  distance  between  adjacent  radial  po¬ 
sitions,  whichever  was  smaller.  Velocity  measurements  were 
based  on  the  motion  of  the  centroid  of  the  image  and  were 
correlated  as  a  function  of  diameter  using  a  linear  least  squares 
fit,  considering  drops  having  diameters  less  than  30  /*m  for 
present  estimates  of  gas  velocities.  This  typically  involved 
correlation  of  50- 150  individual  velocity  determinations.  The 
value  given  by  the  fit  at  a  diameter  of  5  which  is  the 
lower  end  of  the  range  that  could  be  resolved,  was  used  as 
the  estimate  of  mean  gas  velocities.  Unfortunately,  sample 
sizes  were  too  small  to  obtain  reliable  estimates  of  velocity 
fluctuations. 

Gradient  broadening  and  bias  errors  of  mean  gas  velocities 
were  not  significant  for  present  conditions,  yielding  experi¬ 
mental  uncertainties  (95%  confidence)  estimated  to  be  less 
than  30%,  largely  governed  by  sampling  limitations  (except 
as  noted  later).  Experimental  uncertainties  (95%  confidence) 


for  Sauter  mean  diameter  (SMD)  and  mass  median  diameter 
(MMD)  were  both  less  than  10%,  largely  governed  by  sam¬ 
pling  limitations,  similar  to  past  work.2 


Theoretical  Methods 

Predictions  of  flow  properties  were  based  on  the  LHF  ap¬ 
proximation.  neglecting  evaporation  of  the  liquid,  similar  to 
Ruff  et  al. 1  Other  major  assumptions  of  the  analysis  are  as 
follows:  steady  (in  the  mean)  axisymmetric  flow  with  no  swirl; 
negligible  kinetic  energy  and  viscous  dissipation  of  the  mean 
flow;  buoyancy  only  affects  the  mean  flow;  and  equal  ex¬ 
change  coefficients  of  all  species  and  phases.  Under  these 
assumptions,  the  flowfield  can  be  found  using  a  simplified 
version  of  the  conserved-scalar  formalism  of  Lockwood  and 
Naguib17  but  based  on  mass-weighted  (Favre)  averages,  fol¬ 
lowing  Bilger.18  Initial  conditions  were  based  on  the  LV  ve¬ 
locity  measurements  of  jet  pioperties  at  the  injector  exit,1-' 
supplemented  by  available  information  on  the  properties  of 
fully  developed  turbulent  pipe  flows  when  appropriate.14 15 

Justification  of  the  assumptions  and  other  details  of  the 
LHF  predictions  can  be  found  in  Ruff  et  al.1  Notably,  the 
formulation  has  been  successfully  calibrated  for  a  variety  of 
constant-  and  variable-density  single-phase  round  jets. 14  Ad¬ 
ditionally,  the  same  formation  has  been  used  successfully  to 
estimate  the  structure  of  turbulent  round  air  jets  injected  into 
water,  which  involves  the  same  density  ratio  as  the  present 
flow.20  Thus,  the  main  issue  of  the  present  evaluation  is  the 
adequacy  of  the  LHF  approximation  in  the  near-injector  re¬ 
gion  for  liquid  injection  into  gases. 

Results  and  Discussion 

Phase  Velocities 
Mean  Velocities 

Measured  mean  velocities  in  the  multiphase  mixing  layer 
for  slug  and  fully  developed  flow  jet  exit  conditions  are  il¬ 
lustrated  in  Figs.  1  and  2.  Velocities  are  plotted  as  a  function 
of  r/x,  the  radial  similarity  variable  for  turbulent  jets,  to  in¬ 
dicate  the  width  of  the  flow;  however,  flow  properties  do  not 
exhibit  similarity  in  the  r/x  coordinate  system.  Three  sets  of 
velocity  measurements  are  shown:  continuous  (gas)-phase  mean 
velocities  found  using  the  phase -discriminating  LV;  dispersed 
(drop)-phase  mean  velocities  for  5-jum-diameter  particles, 
measured  using  double-flash  holography  (which  are  taken  to 
be  representative  of  mean  gas-phase  velocities  in  the  region 
where  LV  was  no  longer  feasible);  and  Favre-averaged  ve¬ 
locities  of  the  mixture  as  a  whole  for  comparison  with  pre¬ 
dictions.  The  Favre-averaged  velocities  were  found  by  sum¬ 
ming  over  the  sample  volume  V  containing  N  drops,  as  follows: 

£  (Trl6)prdl,Up,  +  (V  ~  2  (tr/6)dp, )p,ut 

u  =  ^ - # -  (l) 

£  (7r/6)p^,  +  (V  -  2  {rrm\,)p, 

»  -  I  i-l 

where  drop  size  and  velocity  determinations  included  the  full 
range  of  the  data  reported  by  Ruff  et  al.1  Mean  gas  velocities 
in  Eq.  (1)  were  obtained  from  the  LV  measurements,  if  avail¬ 
able,  or  from  holography  measurements  (based  on  5-^m-di- 
ameter  drops),  otherwise.  The  range  of  positions  of  the  liquid 
surface,  observed  using  holography,  is  indicated  by  the  cross- 
hatched  regions  in  the  figures.  Finally,  LHF  predictions  of 
Favre-averaged  velocities  are  also  illustrated  on  the  plots. 
Two  predictions  are  shown  for  slug  flow,  ignoring  and  allow¬ 
ing  for  boundary-layer  development  along  the  walls  of  the 
injector  passage,  the  latter  using  a  flow  development  length 
of  Lid  =  5.  These  limits  should  bound  the  properties  of  the 
slug  flow  injector;  however,  differences  between  the  two  pre¬ 
dictions  are  not  very  significant. 
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Fig.  I  Mean  phase  velocities  for  slug  flow  and  atomization  breakup. 

The  results  for  slug  flow  jet  exit  conditions  (Fig.  1)  exhibit 
reasonably  good  agreement  between  mean  velocity  measure¬ 
ments  found  by  LV  and  particle  tracking  in  the  region  where 
they  overlap.  Gas  velocities  remain  quite  low  in  the  mixing 
layer,  and  only  increase  slightly  near  the  liquid  surface  as 
distance  from  the  injector  exit  increases.  Near  the  injector 
exit,  Favre-averaged  velocities  are  significantly  greater  than 
gas  velocities  except  near  the  outer  edge  of  the  mixing  layer. 
Farther  from  the  injector,  however,  differences  between  phase 
velocities  are  only  significant  near  the  liquid  surface.  Never¬ 
theless,  separated-flow  effects  are  important  throughout  the 
mixing  layer  and  LHF  predictions  of  mean  velocities  are  not 
satisfactory. 

The  results  for  the  fully  developed  jet  exit  conditions  (Fig. 
2)  exhibit  greater  differences  between  mean  velocities  found 
by  LV  and  particle  tracking  than  for  slug  flow  jet  exit  con¬ 
ditions.  In  general,  velocities  measured  by  particle  tracking 
are  biased  upward  from  the  LV  results.  This  is  felt  to  be  the 
result  of  the  breakup  of  larger  liquid  elements  that  are  more 
common  for  fully  developed  flow  than  for  slug  flow,  gener¬ 
ating  small  drops  with  relatively  high  initial  velocities.  Mea¬ 
surement  of  the  velocity  of  these  drops  before  they  relax  to 
the  gas  velocity  would  tend  to  bias  present  velocity  measure¬ 
ments  toward  higher  velocities.  The  extent  of  this  effect  is 
difficult  to  quantify;  therefore,  the  particle  velocities  of  Fig. 
2  are  at  best  representative  of  an  upper  bound  on  gas  veloc¬ 
ities.  Similar  to  slug  flow,  gas  velocities  remain  relatively  low 
throughout  the  mixing  layer,  Favre-averaged  velocities  are 
generally  significantly  greater  than  gas  velocities  indicating 
significant  effects  of  separated  flow,  and  LHF  predictions  are 
not  very  satisfactory  as  a  result. 

Cos  Velocity  Fluctuations 

Measured  and  predicted  gas  velocity  fluctuations  are  illus¬ 
trated  in  Figs.  3  and  4  for  slug  and  fully  developed  flow  jet 


Fig.  2  Mean  phase  velocities  for  slug' Row  and  atomization  breakup. 

exit  conditions.  The  measurements  are  limited  to  streamwise 
and  radial  velocity  fluctuations  found  by  LV,  which  generally 
only  covers  the  outer  half  of  the  mixing  layer.  Predictions  are 
based  on  the  isotropic  approximation,  ii"2  =  v"2  =  2kl3,  and 
are  Favre  averages  under  the  LHF  assumption.  Measured 
velocity  fluctuations  are  nearly  isotropic  near  the  outer  edge 
of  the  mixing  layer,  which  is  similar  to  the  behavior  of  single¬ 
phase  jets.  Streamwise  velocity  fluctuations  are  substantially 
greater  than  radial  velocity  fluctuations,  however,  as  the  spray 
becomes  dense  (this  is  particularly  evident  for  xld  =  12.5). 
Such  behavior  is  typical  of  dense  sprays  in  regions  where 
relative  velocities  are  significant.5  21  22  Other  than  indicating 
that  velocity  fluctuations  should  increase  in  the  dense  portions 
of  the  mixing  layer,  in  a  qualitative  way,  the  LHF  predictions 
of  velocity  fluctuations  are  not  very  useful. 

Liquid  Breakup 

Drop-Size  Distributions 

Table  1  is  a  summary  of  drop-size  properties  for  nontur- 
bulent  and  turbulent  jet  exit  conditions.  These  properties  were 
relatively  independent  of  position  along  the  liquid  surface, 
and  over  the  flow  cross  section  if  the  region  near  the  liquid 
surface  was  excluded.  Thus,  SMD  and  MMD/SMD  ratio  have 
been  grouped  accordingly  in  Table  1  with  properties  over  the 
cross  section  represented  by  results  at  xld  =  100.  Gas  veloc¬ 
ities  along  the  liquid  surface  varied  irregularly  due  to  exper¬ 
imental  uncertainties  so  that  only  the  range  of  values  is  sum¬ 
marized  in  Table  1.  Velocities  in  the  liquid  core  remain  close 
to  jet  exit  velocities,  ca.  56  m/s.1 

For  each  jet  exit  turbulence  level,  SMD  along  the  liquid 
surface  is  substantially  greater  than  over  the  cross  section, 
providing  evidence  for  significant  effects  of  secondary  breakup. 
Additionally,  liquid  turbulence  substantially  increased  drop 
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Fig.  3  Gas-phase  velocity  fluctuations  for  slug  flow  and  atomization 
breakup. 


sizes  everywhere  in  comparison  to  nonturbulent  conditions: 
roughly  7:1  after  primary  breakup  at  the  surface,  and  roughly 
2:1  over  the  cross  section.  In  view  of  the  larger  drop  sizes  for 
a  turbulent  liquid,  it  is  somewhat  surprising  that  this  flow' 
mixes  fastest  (cf.  Figs.  1  and  2).  This  behavior  is  due  to  faster 
primary  breakup  for  the  turbulent  liquid,1  perhaps  with  en¬ 
hanced  radial  drop  velocities  due  to  radial  velocity  fluctua¬ 
tions  in  the  liquid. 

In  spite  of  effects  of  turbulence  on  primary  breakup,  as 
well  as  secondary  breakup,  the  MMD/SMD  ratio  is  roughly 
1.2  for  all  conditions  in  Table  1.  This  value  agrees  with  ob¬ 
servations  for  a  wide  variety  of  nonevaporating  dilute  sprays.23 
This  motivated  examination  of  whether  the  drop-size  distri¬ 
butions  in  the  dense  sprays  satisfy  Simmons’  universal  root- 
normal  drop-size  distribution  function.23  These  results  are  il¬ 
lustrated  in  Fig.  5,  considering  every  position  measured  for 
both  nonturbulent  and  turbulent  jet  exit  conditions.  It  is  ev¬ 
ident  that  the  root-normal  distribution  with  MMD/SMD  = 
1.2  provides  a  reasonably  good  fit  of  all  the  measurements. 
Notably,  except  for  x/d  =  12.5,  where  effects  of  flow  devel¬ 
opment  are  strongest,  individual  evaluations  of  the  MMD/ 
SMD  ratio  yield  a  value  of  roughly  1.2  as  well.  The  correlation 
of  Fig.  5  is  very  helpful  because  the  entire  drop-size  distri¬ 
bution  can  be  related  in  a  simple  manner  to  SMD.23  However, 
an  explanation  of  this  behavior  and  more  evidence  that  it  is 
generally  suitable  for  dense  sprays  is  needed. 

Nonturbulent  Primary  Breakup 

Present  results  indicate  different  breakup  mechanisms  for 
nonturbulent  and  turbulent  liquids  at  otherwise  identical  con¬ 
ditions  in  the  atomization  breakup  regime  with  low  Oh.  It  is 
generally  agreed  that  atomization  breakup  at  low  Oh  involves 
a  stripping  mechanism  for  nonturbulent  liquids.2  3  5  Among 


Fig.  4  Gas-phase  velocity  fluctuations  for  fully  developed  flow  and 
atomization  breakup. 


Table  1  Drop-size  properties* 


Liquid  surface,  12.5 

Cross  section. 

s  x/d  S  100 

x/d 

=  100 

u„,  SMD. 

MMD/ 

SMD. 

MMD/ 

Row 

m/s"  utn 

SMD 

/xm 

SMD 

Slug 

3-9  174 

1.21 

95 

1.18 

Fully  developed 

8-15  1330 

1.20 

220 

1.22 

■SMD  and  MMD  are  based  on  equivalent  spherical  drops  as  described  in  Ruff 
et  al.2 

"Gas  velocity  estimated  from  5-Mm-diam  drops  near  surface. 


others,  Levich24  has  proposed  a  formula  to  estimate  drop  sizes 
due  to  primary  breakup  for  these  conditions  based  on  Taylor’s23 
theory  of  aerodynamically  induced  growth  of  surface  waves, 
assuming  that  the  average  diameter  of  drops  formed  at  the 
surface  is  proportional  to  the  wavelength  of  the  unstable  sur¬ 
face  waves  having  the  highest  growth  rates.  Using  SMD  to 
represent  the  average  drop  diameter  then  yields 

SMD  =  67rCa<7/[pf(tiA  -  ii„)2]  (2) 

where  Cs  is  a  constant  of  proportionality  near  unity. 

For  present  conditions,  df,  is  roughly  equal  to  iq,  =  56.3 
m/s,  based  on  LHF  computations,  as  well  as  measurements 
of  the  velocity  of  protuberances  from  the  liquid  surface.  Then, 
matching  the  primary  breakup  SMD  in  Table  1,  using  an 
average  value  of  ii,,,  yields  C*  =  0.4.  Since  C*  is  expected 
to  be  on  the  order  of  unity,  the  aerodynamic  primary  breakup 
correlation  of  Eq.  (2)  appears  to  be  reasonable  for  nontur¬ 
bulent  conditions  within  dense  sprays.  However,  additional 
evaluation  of  the  expression  for  other  test  conditions  is  clearly 
needed. 
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Fig.  5  Drop-size  distributions  at  liquid  surface  and  across  mixing 
layer  for  slug  and  fully  developed  flow  and  atomization  breakup. 


Turbulent  Primary  Breakup 

The  much  larger  drop  sizes  after  primary  breakup  for  tur¬ 
bulent  jet  exit  conditions  are  best  explained  by  a  turbulent 
breakup  mechanism  with  relatively  small  aerodynamic  effects. 
Lee  and  Spencer24  report  early  observations  of  the  turbulent 
breakup  mechanism,  based  on  tests  at  low  ambient  pressures 
where  aerodynamic  effects  are  negligible.  Reitz27  and  Reitz 
and  Bracco2*  also  identify  breakup  induced  by  liquid-phase 
turbulence  as  one  of  the  possible  mechanisms  contributing  to 
the  atomization  process.  This  mechanism  involves  turbulent 
distortion  of  the  liquid  surface  to  form  drops,  with  surface 
tension  acting  to  stabilize  the  distortion,  so  that  maximum 
drop  sizes  should  be  comparable  to  the  radial  integral  scale 
of  the  liquid  turbulence.  Taking  the  SMD  to  be  representative 
of  maximum  drop  sizes,  this  implies 

SMD/A  =  C,'  (3) 

where  C,'  is  a  constant  on  the  order  of  unity.  Laufer14  29  finds 
Md  =  0.10-0.15  for  turbulent  pipe  flow  while  properties 
within  the  liquid  core  do  not  vary  significantly  from  jet  exit 
conditions  based  on  the  LHF  computations.1  Then,  based  on 
the  SMD  measurements  in  Table  1  for  turbulent  jet  exit  con¬ 
ditions,  Eq.  (3)  yields  C,'  =  0.9- 1.3,  with  the  range  resulting 
from  Laufer’s29  range  of  A Jd.  Additionally,  the  ratio  of 
streamwise  to  cross-stream  integral  scales  for  turbulent  pipe 
flow  is  roughly  4: l,29  which  is  comparable  to  the  mass-aver¬ 
aged  ellipticity  (ratio  of  maximum  to  minimum  liquid  element 
dimensions)  observed  for  turbulent  breakup  conditions  near 
the  liquid  surface.1  Thus,  the  turbulent  breakup  mechanism 
provides  a  plausible  description  of  the  large-scale  features  of 
the  sprays  after  primary  breakup. 

Additional  evidence  supporting  the  turbulent  breakup 
mechanism  can  be  found  by  considering  the  criterion  for  this 
mechanism  to  be  active.  This  can  be  estimated,  similar  to 
transition  from  the  drip  to  the  jetting  regime  of  liquid  jets  in 
gases,30  by  equating  the  momentum  available  from  turbulent 
fluctuations  of  scale  dp  to  the  surface  tension  forces  resisting 
the  formation  of  a  drop  of  similar  size.  Using  the  relationship 
between  velocity  fluctuations  and  eddy  size  within  the  inertial 
region  of  the  turbulence  spectrum,29  the  expression  for  the 
range  of  drop  sizes  that  can  be  formed  by  the  turbulent  pri¬ 
mary  breakup  mechanism  is 

pr\V}(dJ\Y'Va  >  6C,  (4) 


where  C,  is  a  constant  on  the  order  of  unity.  Taking  C,  =  1. 
Eq.  (4)  implies  that  the  minimum  drop  diameter  resulting 
from  turbulent  breakup  for  the  conditions  of  Table  1  is  ca. 
250  ^m,  which  is  in  good  agreement  with  observations.  In 
contrast,  Eq.  (4)  is  not  satisfied  even  when  d„  =  A  for  the 
first  wind-induced  test  condition  for  the  present  sprays. 1 2  This 
also  corresponds  to  observations  because  the  liquid  surface 
only  was  roughened  by  turbulent  distortion  for  this  test  con¬ 
dition  and  no  drop  breakup  occurred.  Thus,  the  turbulent 
breakup  mechanism  is  consistent  with  present  observations 
and  deserves  additional  study  over  a  wider  range  of  test  con¬ 
ditions. 

Equations  (3)  and  (4)  also  help  to  put  observations  of  effects 
of  liquid  turbulence  on  primary  breakup  for  the  present  large- 
scale  sprays  into  perspective.  For  example,  small  injector  di¬ 
ameters  imply  relatively  small  maximum  drop  sizes  from  this 
mechanism,  while  the  aerodynamic  breakup  correlation.  Eq. 
(2),  is  essentially  independent  of  injector  diameter  (as  long 
as  Oh  is  sufficiently  low.  Thus,  aerodynamic  breakup  is  likelv 
to  be  more  dominant  for  conventional  sized  injectors,  because 
it  yields  larger  drops  than  turbulent  breakup,  even  when  jet 
exit  conditions  are  fully  turbulent.  Additionally,  most  prac¬ 
tical  injectors  have  small  length-to-diameter  ratios  so  that 
liquid  turbulence  is  not  highly  developed  even  when  jet  Rey¬ 
nolds  numbers  are  large.  Based  on  these  considerations,  tur¬ 
bulent  primary  breakup  is  probably  not  a  major  feature  of 
most  practical  sprays,  as  suggested  by  Reitz27  and  Reitz  and 
Bracco.28 

Secondary  Breakup 

Criteria  for  secondary  breakup  are  frequently  stated  in  terms 
of  a  drop  Weber  number,  defined  as  follows: 

Wep,  =  PAK  -  u,y/a  (5) 

One  criterion,  discussed  by  Clift  et  al.,30  is  based  on  water 
drops  accelerated  by  shock  waves,  yielding: 

Wep,  >  We*,  =  6.5  (6) 

Another  criterion,  discussed  by  Pruppacker  and  Klett,31  was 
obtained  for  liquid  drops  falling  in  still  air: 

We„  >  We •„  =  8/ Co  (7) 

For  present  test  conditions,  drop  drag  coefficients  are  in  the 
range  of  0.8- 1.2  for  the  largest  drops,  based  on  the  standard 
drag  coefficient  for  spheres5;  therefore,  Eqs.  (6)  and  (7)  yield 
similar  values  of  We’,. 

In  order  to  quantify  the  potential  for  secondary  drop  breakup 
in  the  mixing  layer,  mass-weighted  Weber  numbers  were  com¬ 
puted  as  follows: 

We,  =  2  ( p,d,K  -  “AAA.  /  S  (8) 

,-i  /  ... 

given  N  drop  measurements  at  a  point.  The  resulting  distri¬ 
butions  of  We,  for  the  sprays  having  slug  and  fully  developed 
jet  exit  conditions  are  plotted  in  Fig.  6.  The  values  of  We, 
are  greatest  near  the  liquid  surface  and  decrease  monotoni- 
cally  with  increasing  distance  from  the  surface.  In  general, 
values  of  We  near  the  liquid  surface  exceed  We’,  for  breakup 
from  Eqs.  (6)  and  (7),  supporting  significant  effects  of  sec¬ 
ondary  breakup  in  the  mixing  layer.  The  larger  liquid  elements 
ejected  from  the  surface  when  the  liquid  is  turbulent  imply 
greater  effects  of  secondary  breakup  for  these  conditions. 

Primary  followed  by  secondary  breakup  is  likely  to  be  a 
general  mechanism  for  dense  sprays.  For  example,  the  value 
of  Wep ,  for  a  drop  having  a  diameter  equal  to  the  SMD  is 
roughly  8  for  aerodynamic  primary  breakup  from  Eq.  (2). 
Thus,  drops  having  diameters  of  the  order  of  the  SMD  and 
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Fig.  6  Favre-averaged  drop  Weber  numbers  for  slug  and  fully  de¬ 
veloped  (low  and  atomization  breakup. 

larger,  which  comprise  more  than  half  the  mass  of  the  liquid 
in  the  spray  for  MMD/SMD  =  1.2,  have  Weber  numbers  that 
exceed  We’t  from  Eqs.  (6)  and  (7).  Naturally,  the  turbulent 
breakup  mechanism,  if  active,  yields  even  larger  values  of 
SMD  and  a  greater  propensity  for  secondary  breakup.  These 
observations,  as  well  as  the  relatively  small  liquid  volume 
fractions  within  the  multiphase  mixing  layer,1  suggest  that 
secondary  breakup  rather  than  drop  collisions  is  a  dominant 
feature  of  dense  sprays. 

Separated  Flow  Effects 

The  results  illustrated  in  Figs.  1  and  2  are  somewhat  mis¬ 
leading  because  small  velocities  give  the  impression  that  sep¬ 
arated  flow  effects  are  small  near  the  edge  of  the  flow.  This 
is  not  the  case,  which  can  be  seen  from  the  plots  of  Favre- 
averaged  separated-flow  factor,  defined  as  follows: 

S  =  (“/  "  us)/uf  (9) 

Distributions  of  S  are  illustrated  in  Fig.  7  for  both  slug  and 
fully  developed  jet  exit  conditions.  Values  of  5  approaching 
unity  are  observed  throughout  the  mixing  layer  near  the  jet 
exit,  and  near  the  liquid  surface  for  all  streamwise  positions. 
Values  of  5  near  the  outer  edge  of  the  flow  tend  to  decrease 
with  increasing  distance  from  the  injector,  however,  ap¬ 
proaching  values  in  the  range  0. 1-0.2  for  fully  developed  jet 
exit  conditions  at  x/d  =  100.  In  spite  of  the  presence  of  larger 
drops  near  the  liquid  surface,  values  of  S  are  generally  lower 
for  fully  developed  than  slug  flow  jet  exit  conditions. 

The  trends  seen  in  Fig.  7  result  from  both  the  generation 
of  drops  at  the  liquid  surface  and  momentum  exchange  be¬ 
tween  the  phases  in  the  mixing  layer.  Velocities  within  the 
liquid  core  remain  near  liquid  injection  velocities  so  that  newly 
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Fig.  7  Favre-averaged  separated-flow  factor  for  slug  and  fully  de¬ 
veloped  flow  and  atomization  breakup. 

formed  drops  near  the  injector  exit  and  along  the  liquid  sur¬ 
face  generally  have  velocities  that  are  significantly  greater 
than  gas  velocities,  yielding  relatively  large  values  of  S.  Drops 
near  the  outer  edge  of  the  flow,  however,  must  have  had 
significant  residence  times  in  the  gas  so  that  turbulent  dis¬ 
persion  can  transport  them  to  this  region.  Therefore,  drops 
near  the  edge  of  the  flow  tend  to  relax  toward  gas  velocities, 
particularly  when  residence  times  are  long  (at  large  x/d).  Ad¬ 
ditionally.  fully  developed  jet  exit  conditions  increase  drop 
concentrations  in  the  mixing  layer  in  comparison  to  slug  flow 
exit  conditions,  causing  greater  acceleration  of  the  gas  within 
the  mixing  layer.  This  compensates  for  the  higher  velocities 
associated  with  larger  drops  produced  by  fully  developed  exit 
conditions  and  tends  to  reduce  S.  Nevertheless,  the  values  of 
5  on  the  order  of  unity  throughout  most  of  the  mixing  layer, 
separated-flow  effects  are  clearly  important,  which  explains 
why  predictions  based  on  the  LHF  approximation  are  poor 
for  present  test  conditions. 

Scaling 

Characteristic  Times 

Characteristic  times  are  considered  in  the  following  to  eval¬ 
uate  secondary  breakup  and  effects  of  separated  flow,  and  to 
help  relate  present  observations  for  large  sprays  to  practical 
pressure-atomized  sprays.  Characteristic  times  of  interest  are 
flow  residence,  liquid  core  residence,  drop  response,  and  drop 
breakup  times.  The  largest  drops  are  of  interest,  which  have 
velocities  of  roughly  Ug-,  therefore,  an  appropriate  flow  resi¬ 
dence  time  is 

Tf  =  XlUg  (10) 

where  x  is  to  be  interpreted  as  the  distance  from  the  position 
where  the  drop  entered  the  mixing  layer  after  primary  breakup. 
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It  is  also  of  interest  to  assess  response  over  the  length  of 
he  liquid  core  as  a  measure  of  the  length  of  the  dense-spray 
1  l0n.  Taylor's23  aerodynamic  breakup  analysis  yields  the 
following  expression  for  the  length  of  the  liquid  core: 

LJd  =  CJpJp,)1'2  (11) 

Chehroudi  et  al."  find  Cc  in  the  range  7-16  for  the  contig¬ 
uous  length  of  liquid  core  for  pressure  atomized  round  jets 
in  still  gases  within  the  atomization  breakup  regime.  Taking 
the  characteristic  velocity  of  large  drops  to  be  u„,  the  char¬ 
acteristic  residence  -  me  of  the  liquid  core  becomes: 

tc  =  LJuq  =  Cc(dJu„)(pf/pt)'2  (12) 

The  largest  drops  have  the  longest  characteristic  response 
times.  The  results  illustrated  in  Fig.  6  suggest  that  the  largest 
drops  have  Weber  numbers  comparable  to  WePt  from  Eq.  (6) 
for  both  laminar  and  turbulent  liquid  conditions,  and  this 
condition  will  be  used  as  a  measure  of  maximum  drop  sizes. 
This  yields  a  critical  drop  size  of 

d’p  =  <TWe;tl{p,u-,)  (13) 

assuming  that  relative  velocities  are  comparable  to  u„  based 
on  present  measurements  for  large  drops.'  A  characteristic 
response  time  for  drops  can  be  defined  in  terms  of  their  rate 
of  deceleration,  as  follows: 

rp  =  -(dup/dx)-‘  =  -up/(dupldt)  (14) 

Virtual  mass  and  Basset  history  forces  generally  can  be  ig¬ 
nored  when  computing  drop  motion  in  sprays  because  p/p, 
»  1,  while  pressures  within  the  mixing  layer  are  constant; 
therefore,  conservation  of  momentum  for  a  drop  yields5: 

dup/dt  =  -3 p,CD(up  -  up)-/(4p^dp)  (15) 

Streamwise  velocities  have  been  assumed  to  be  large  in  com¬ 
parison  to  cross-stream  velocities  to  obtain  Eq  (15).  which 
is  justified  due  to  the  boundary-layer  character  of  the  mixing 
layer.  Substituting  Eq.  (15)  into  Eq.  (14),  assuming  ut  «  up 
"•  un  as  before,  and  using  Eq.  (13)  to  identify  the  critical  drop 
size,  then  yields: 

t,  =  4We;np^C'Dp\ul)  (16) 

Because  Ca  does  not  vary  appreciably  for  large  drops  in  sprays,5 
Eq.  (16)  shows  that  increased  gas  densities  and  jet  exit  ve¬ 
locities  reduce  drop  response  times,  both  due  to  effects  of 
breakup  which  reduces  drop  sizes  and  the  relative  increase  of 
drag  in  comparison  to  inertia  of  the  drops  for  a  particular  size 
[see  Eq.  (15)]. 

Finally,  the  characteristic  time  required  for  secondary 
breakup  is  an  important  feature  of  dense  sprays.  As  noted 
earlier,  drop  Weber  numbers  are  not  much  greater  than  the 
breakup  limit,  so  that  deformation  and  stripping  breakup  are 
the  main  mechanisms  of  secondary  breakup.32  Therefore, 
adopting  the  correlation  of  Ranger  and  Nicholls33  for  stripping 
breakup  yields 

=  Cbdp(p/p,)'V(up  -  “,)  (17) 

where  the  empirical  factor  Cb  *  4.  Taking  dp  =  dp  and  u , 
«  up  -  ,  as  before,  Eq.  (17)  becomes: 

r»  =  C*We>(p,/p,)3'2/(p,u3)  (18) 


Characteristic  Time  Ratios 

A  variety  of  characteristic  time  ratios  can  be  formed  to 
highlight  various  properties  of  dense  sprays.  Forming  the  ratio 
tJtp  from  Eqs.  (16)  and  (18)  yields: 

Tb'TP  =  (3CC»'4)(p,  p,)1 2  U9) 

Small  values  of  r*.  tc  imply  that  secondary  breakup  times  are 
small  in  comparison  to  the  time  required  tor  drops  to  relax 
to  the  local  gas  velocity.  Taking  C‘n  =  1.  as  before,  implies 
that  r,;r.  is  ca.  0. 1  for  present  test  conditions  so  that  secondary 
breakup  is  relatively  fast,  which  agrees  with  observations  that 
the  largest  drops  were  close  to  the  liquid  surface  Under  these 
conditions,  the  outcome  of  secondary  breakup  has  the  greatest 
effect  on  dense-spray  properties.  At  high  pressures,  however. 
tJtp  =  1  and  breakup  extends  over  distances  required  for 
drops  to  approach  gas  velocities;  then,  drop  breakup  must  be 
treated  like  a  finite  rate  process,  similar  to  effects  of  drag  and 
vaporization  on  drop  motion. 

Another  perspective  on  secondary  breakup  can  be  found 
from  the  ratio  of  Eqs.  (12)  and  (18) 

rblrc  =  CbWe;,/(CeWe,)  (20) 

which  is  the  ratio  of  secondary  breakup  to  liquid  core  resi¬ 
dence  time.  For  present  test  conditions,  r Jr,  ca.  0.005,  so 
that  breakup  times  are  short  in  comparison  to  residence  times 
in  the  mixing  layer.  This  also  agrees  with  observations  that 
the  large  drops  resulting  from  primary  breakup  were  associ¬ 
ated  with  the  region  near  the  liquid  surface.  However,  this  is 
not  always  the  case.  For  example,  transition  from  wind-in¬ 
duced  to  atomization  breakup  requires  Weg  >  We *  =  13- 
40. 23  Near  the  transition  limit,  t„Itc  is  in  the  range  0. 1-0.3. 
which  implies  that  secondary  breakup  proceeds  rather  slowly 
within  the  mixing  layer  (particularly  because  new  drops  are 
formed  all  along  the  liquid  core).  Under  these  conditions, 
rates  of  secondary  breakup  should  control  processes  in  dense 
sprays,  somewhat  like  drop  drag  and  vaporization  tends  to 
control  processes  in  dilute  sprays. 

Insight  concerning  effects  of  separated  flow  can  be  obtained 
from 

rr’Tp  =  ( 3C  ydWepg)(p,!pf)-(p,xu}yCr )  (21) 

where  large  values  of  t/tp  correspond  to  regions  of  the  flow 
where  effects  of  separated  flow  are  small.  For  present  flows, 
drops  continue  to  be  formed  at  the  liquid  surface  throughout 
the  region  where  measurements  were  made.  Therefore,  char¬ 
acteristic  distances,  x,  and  the  corresponding  values  of 
remain  small,  which  is  consistent  with  the  large  values  of  (ti, 
-  tig)/ufseen  in  Fig.  7  in  this  region.  In  contrast,  drops  near 
the  edge  of  the  flow  originate  near  the  jet  exit.  Then,  esti¬ 
mating  the  empirical  parameters  in  Eq.  (21)  as  before.  rrrp 
varies  in  the  range  0.5-4  as  xld  varies  in  the  range  12.5-100. 
This  implies  significant  effects  of  separated  flow  near  the  jet 
exit,  with  separated-flow  effects  decreasing  near  the  edge  of 
the  flow  at  larger  streamwise  distances,  which  is  also  con¬ 
sistent  with  the  results  illustrated  in  Fig.  7. 

For  a  round  jet.  most  of  the  liquid  enters  the  mixing  layer 
near  the  jet  exit,  and  large  values  of  the  ratio 

tc/t„  =  ( 3  Cc  C  q/4  We'p,){pjp, ) '  -Wef  (22) 

implies  that  separated-flow  effects  are  small  within  the  dense 
spray  region.  The  effect  of  the  density  ratio  p/p,  is  somewhat 
reduced  for  tJ t„  in  comparison  to  t/tp,  because  the  length 
of  the  liquid  core  is  reduced  as  this  density  ratio  increases 
through  Eq.  (11).  For  present  test  conditions,  LJd  is  in  the 
range  200-  460  from  Eq.  (11)  while  tJtp  is  in  the  range  11- 
25  from  Eq.  (22).  Thus,  the  values  of  tJtp  imply  modest 
effects  of  separated  flow  over  the  length  of  the  liquid  core, 
although  separated-flow  effects  are  still  significant  for  the 
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relatively  low  range  of  xld  ( xid  s  100)  where  present  mea¬ 
surements  were  made.  More  generally,  Eq.  (22)  implies  that 
large-diameter  jets,  like  those  of  the  present  experiments, 
have  reduced  effects  of  separated  flow  because  of  their  large 
We,  in  comparison  to  practical  sprays  having  similar  jet  exit 
velocities.  Thus,  the  LHF  limit  is  approached  at  large  We,. 
with  this  limit  approached  more  rapidly  as  pr'p ,  increases,  as 
long  as  Oh  is  sufficiently  small.  Additional  measurements, 
however,  are  needed  to  confirm  the  trends  suggested  bv  Eqs. 
(19-22). 


Conclusions 

The  multiphase  mixing  layer  in  the  near-injector  region  of 
large-scale  pressure-atomized  sprays  was  investigated,  con¬ 
sidering  both  nonturbulent  slug  and  fully  developed  turbulent 
pipe  flow  at  the  jet  exit.  The  major  conclusions  of  the  study 
are  as  follows: 

1 .  Present  measurements,  combined  with  the  earlier  results 
of  Ruff  et  al.*  for  the  same  flows,  tend  to  support  the  classical 
view  of  atomization:  a  surprisingly  dilute  multiphase  mixing 
layer  (liquid  volume  fractions  less  than  1%)  dominated  by 
secondary  breakup  and  separated-flow  processes,  while  ef¬ 
fects  of  collisions  appear  to  be  small. 

2.  Use  of  the  LHF  approximation  was  not  effective  for 
estimating  the  properties  of  the  multiphase  mixing  layers  for 
present  test  conditions,  because  large  drops  are  continuously 
formed  along  the  liquid  surface  and  require  significant  times 
to  relax  to  gas  velocities.  Characteristic  time  considerations 
suggest  that  locally  homogenous  flow  is  approached  at  large 
Wef,  with  this  limit  approached  more  rapidly  as  p,ip,  increases. 

3.  For  nonturbulent  jet  exit  conditions,  the  SMD  for  pri¬ 
mary  breakup  agreed  with  aerodynamic  breakup  theory,24 
with  C,  =  0.4.  For  turbulent  jet  exit  conditions,  SMD  for 
primary  breakup  was  explainable  by  turbulent  breakup  theory 
with  th-.  SMD  roughly  equal  to  the  radial  integral  scale  of  the 
liquid  turbulence. 

4.  Prop-size  distributions  after  aerodynamic  primary 
breakup,  turbulent  primary  breakup,  and  secondary  breakup 
all  exhioited  MMD/SMD  ratios  of  roughly  1.2  and  agreed 
reasonably  well  with  Simmons's  root-normal  drop-size  distri- 
butioi.  -J  This  is  a  useful  simplification  because  the  entire 
drop-size  distribution  is  known  if  a  single  moment,  like  SMD, 
can  be  found. 

5.  An  expression  for  the  range  of  drop  sizes  that  can  be 
formed  from  the  turbulent  primary  breakup  mechanism  was 
deve!o->ed  and  indicated  the  dependence  of  drop  size  on  the 
integral  scale  and  intensity  of  liquid-phase  turbulence.  The 
expre*  ion  was  generally  supported  by  present  measurements. 
Howt  >  er,  the  importance  of  this  breakup  mechanism  for  prac¬ 
tical  i.  .ectors  having  small  diameters  is  questionable. 

6.  Eased  on  present  measurements,  more  than  half  of  the 
mass  of  the  spray  is  contained  in  drops  that  are  unstable  to 
second  try  breakup  after  aerodynamic  breakup,  while  tur¬ 
bulent  primary  breakup  yields  larger  drops  that  have  a  greater 
propei  sity  for  secondary  breakup.  Thus,  secondary  breakup 
appears  to  be  an  intrinsic  process  of  dense  sprays  in  the  atom¬ 
ization  '.pray  breakup  regime. 

The  conclusions  are  based  on  limited  measurement  for  liq¬ 
uid  jets  having  large  diameters  which  have  relatively  slow  flow 
development  rates  in  comparison  to  practical  injectors  and 
are  only  provisional  pending  additional  experimental  evalu¬ 
ation. 
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A  theoretical  and  experimental  study  of  the  near  jet-exit  region  of  nonevaporating  round  liquid  jets  in  still 
gases  is  described,  emphasizing  effects  of  ambient  gas  density  in  the  atomization  breakup  regime  where  liquid 
breakup  begins  right  at  the  jet  exit.  Mean  liquid  volume  fraction  distributions  were  measured  for  9.5-mm-diam 
water  jets  in  still  air  at  pressures  of  1-8  atm.  Mixing  was  strongly  affected  by  the  gas/liquid  density  ratio  and 
the  degree  of  flow  development  at  the  jet  exit,  with  the  largest  gas/iiquid  density  ratio  and  fully  developed 
turbulent  pipe  flow  yielding  the  fastest  mixing  rates.  Flow  properties  were  predicted  using  the  locally  homo¬ 
geneous  flow  approximation,  where  relative  velocities  between  the  phases  are  assumed  to  be  small  in  comparison 
to  mean  flow  velocities.  Predictions  were  in  good  agreement  with  measurements,  including  representation  of 
effects  of  gas/liquid  density  ratio  and  flow  development  at  the  jet  exit,  but  only  at  relatively  high  mixture 
fractions.  In  contrast,  separated  flow  effects  caused  predictions  to  overestimate  rates  of  flow  development  at 
low  mixture  fractions. 


Nomenclature 

d  -  jet-exit  diameter 

/  =  mixture  fraction 

k  =  turbulence  kinetic  energy 

L  =  injector  passage  length 

Oh  =  Ohnesorge  number.  do-)1 : 

p  -  pressure 

r  -  radial  distance 

Re,  -  jet  Reynolds  number.  p,du„  pt 

u  =  streamwise  velocity 

v  =  radial  velocity 

w  =  tangential  velocity 

We,u  ~  Weber  number  based  on  phase  i  and  the  jet 
diameter.  p,du;,  o 
x  =  streamwise  distance 

a  =  volume  fraction 

f  =  rate  of  dissipation  of  turbulence  kinetic  energy 

p.  =  molecular  viscosity 

p  =  density 

<r  =  surface  tension 

<t>  =  generic  property 

Subscripts 

c  =  centerline  value 

/  =  liquid-phase  property 

g  =  gas-phase  property 

o  =  jet-exit  conditions 

*  =  ambient  conditions 

Superscripts 

(  ),(  )'  =  time-averaged  mean  and  rms  fluctuating  quantities 
(").(  )”  =  Favre-averaged  mean  and  rms  fluctuating  quan¬ 
tities 
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Introduction 

N  experimental  and  theoretical  investigation  of  the  near 
jet-exit  region  of  nonevaporating  round  liquid  jets  in  still 
gases  is  described.  This  flow  is  of  interest  for  a  variety  of 
atomization  and  gas/liquid  mixing  processes;  it  also  merits 
study  as  the  multiphase  counterpart  of  the  single-phase  tur¬ 
bulent  jet.  The  objective  of  the  research  was  to  extend  earlier 
studies  of  round  water  jets  in  still  room  air.1  ’  to  consider 
effects  of  ambient  gas  density  on  the  structure  and  mixing 
properties  of  the  flow.  Measurements  of  liquid  volume  frac¬ 
tion  distributions  were  completed  for  water  jets  in  still  air  at 
various  ambient  pressures.  Tests  were  limited  to  atomization 
breakup,  where  liquid  breakup  into  drops  and  the  develop¬ 
ment  of  a  multiphase  miving  layer  begins  right  at  the  jet  exit, 
because  this  regime  is  most  important  for  practical  applica¬ 
tions.1  ■  Jet-exit  conditions  included  both  fully  developed  tur¬ 
bulent  pipe  flow  and  low-turbulence  intensity  slug  flow  due 
to  the  known  importance  of  the  degree  of  flow-  development 
at  the  jet  exit  on  flow  properties. 1  * ' *  The  measurements  were 
used  to  evaluate  predictions  of  effects  of  gas.  liquid  density 
ratio  on  flow  properties,  found  from  an  earlier  model  based 
on  the  locally  homogeneous  flow  ( LHF)  approximation  where 
relative  velocities  between  the  phases  are  assumed  to  be  small 
in  comparison  to  mean  flow  velocities.1 

Only  the  main  features  of  past  work  are  considered .  because 
reviews  treating  multiphase  jets  and  dense  sprays  have  ap¬ 
peared  recently. 7H  Past  measurements  have  established  the 
main  features  of  the  near  injector  region  of  liquid  jets  in  still 
gases  within  the  atomization  breakup  regime.  The  flow  con¬ 
sists  of  a  liquid  core,  much  like  the  potential  core  of  a  single¬ 
phase  jet.  surrounded  by  a  dispersed  drop  gas  mixing  layer. 
Conductivity  probes  have  been  used  to  study  the  le-gth  of 
the  liquid  core  at  various  ambient  gas  densities. " 1 The  liquid 
core  tends  to  become  shorter  at  high  ambient  gas  densities 
due  to  increased  entrainment  rates.  Nevertheless,  the  liquid 
core  and  its  multiphase  mixing  layer  are  prominent  features 
of  the  flow  for  typical  gas/liquid  density  ratios,  extending  200- 
400  jet-exit  diameters  for  liquid  jets  in  still  gases  at  atmos¬ 
pheric  pressure.1' 

Recent  measurements  of  Ruff  et  al.‘~'  have  yielded  some 
information  about  the  structure  of  the  near  jet-exit  region  of 
liquid  jets  in  still  gases.  These  experiments  involved  large  (9  5- 
and  19.1-mm  initial  diameters)  water  jets  injected  in  still  air 
at  atmospheric  pressure.  Measurements  included  liquid  vol¬ 
ume  fraction  distributions  using  gamma-rav  absorption,  drop 
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sizes  and  velocities  in  the  mixing  layer  using  double-pulse 
holography,  and  air  entrainment  rates  using  laser  velocimetry . 
The  measurements  of  mean  liquid  volume  fraction  distribu¬ 
tions  showed  that  the  rate  of  development  of  the  flow  was 
affected  by  both  the  breakup  regime  and  the  degree  of  flow- 
development  at  the  jet  exit:  atomization  breakup  and  fully- 
developed  turbulent  pipe  flow  at  the  jet  exit  yielded  the  fastest 
mixing  rates.  Predictions  based  on  the  LHF  approximation 
were  effective  for  atomization  breakup  when  mean  liquid  vol¬ 
ume  fractions  were  greater  than  0  2.  however,  they  progres¬ 
sively  failed  as  the  flow  became  dilute  and  as  jet-exit  condi¬ 
tions  approached  the  wind-induced  breakup  regime.  Difficulties 
with  the  LHF  approach  were  identified  by  drop  size  and  ve¬ 
locity  measurements  in  the  mixing  layer,  which  showed  sig¬ 
nificant  effects  of  separated  flow  in  regions  where  liquid  vol¬ 
ume  fractions  were  low.  Thus,  the  success  of  the  LHF 
approximation  at  large  liquid  volume  fractions  was  attributed 
to  the  small  proportion  of  the  momentum  associated  with  the 
gas  phase  at  these  conditions  due  to  the  small  gas,  liquid  den¬ 
sity  ratio  of  the  flow.  Similarly,  entrainment  rates  are  gov¬ 
erned  mainly  by  processes  within  the  dilute  portions  of  the 
mixing  layer:  therefore.  LHF  predictions  generally  overesti¬ 
mated  entrainment  rates  due  to  separated-flow  phenomena 
in  the  mixing  layer  However,  the  findings  suggested  im¬ 
proved  performance  of  LHF  methods  at  large  He.,  with  this 
limit  approached  more  rapidly  at  large  gas  liquid  density  ratios. 

The  objective  of  the  present  investigation  was  to  explore 
effects  of  gas  liquid  density  ratio  on  the  mixing  properties  of 
liquid  jets  in  still  gases  and  the  adequacy  of  predictions  of 
flow  properties  based  on  the  LHF  approximation.  Experi¬ 
ments  involved  measurements  of  liquid  volume  fraction  dis¬ 
tributions  within  large-scale  (9.5-mm  initial  jet  diameter)  water 
jets  in  still  air  at  ambient  pressures  of  1-8  atm.  The  mea¬ 
surements  were  compared  with  predictions  based  on  the  ear¬ 
lier  LHF  approach.1-  to  determine  whether  this  methodology 
could  treat  effects  of  varying  gas/  liquid  density  ratios. 

The  paper  begins  with  brief  descriptions  of  experimental 
methods  and  the  LHF  computations.  Results  are  then  de¬ 
scribed.  considering  liquid  volume  and  mass  fraction  distri¬ 
butions  in  turn. 

Experimental  Methods 

Apparatus 

Experimental  methods  were  similar  to  Ruff  et  al. .'  except 
that  the  flow  was  contained  within  a  large  pressure  vessel  so 
that  the  pressure  of  the  ambient  air  could  be  changed.  A 
sketch  of  the  apparatus  appears  in  Fig.  1.  The  arrangement 
involves  a  steady  water  jet  injected  vertically  downward  in 
still  air  within  a  large  windowed  pressure  vessel  (1.5-m  di¬ 
ameter  x  4.5-m  long  with  two  pairs  of  opposite  windows 
having  maximum  opening  diameters  of  250  mm).  Access  to 
the  interior  of  the  pressure  vessel  was  provided  by  a  manhole 
at  the  top  of  the  vessel.  A  second  port  at  the  top  provided 
entry  points  to  supply  water  to  the  jet  and  electrical  cables 
needed  to  control  the  jet  traversing  system. 

City  water  was  fed  to  the  injector  using  a  centrifugal  pump. 
The  water  was  collected  in  the  bottom  of  the  tank  with  the 
liquid  level  maintained  by  outflow  to  a  drain.  A  grate  located 
above  the  liquid  level  reduced  splash  back  into  the  test  area. 
The  rate  of  water  flow  was  adjusted  using  a  bypass  system 
and  was  measured  using  a  turbine  flow  meter  that  was  cali¬ 
brated  by  collecting  water  for  timed  intervals. 

The  water  injectors  had  exit  diameters  of  9.5  mm  and  were 
the  same  as  those  used  by  Ruff  et  al.:1  one  yielding  slug  flow 
with  low  turbulence  intensities,  the  other  yielding  fully  de¬ 
veloped  turbulent  pipe  flow.  The  slug  flow  injector  consisted 
of  a  honeycomb  flow  straightener  ( 1  6-mm  cells.  25-mm  long) 
and  two  screens  to  calm  the  flow  ( 16  x  16  square  mesh.  0. 18- 
mm-diam  wire)  followed  by  a  13.6:1  area  contraction  to  the 
jet  exit.  The  contraction  was  designed  following  Smith  and 
Wang"  to  yield  uniform  velocities  at  the  jet  exit.  The  fully 
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Fig.  t  Sketch  of  the  variable  gas  density  apparatus. 

developed  flow  injector  had  the  same  flow  straightener  and 
contraction  area  ratio,  but  the  contraction  was  followed  by  a 
constant  area  passage  41  jet-exit  diameters  long. 

Instrumentation  was  mounted  rigidly  so  that  flow  structure 
was  measured  by  traversing  the  jet.  Horizontal  traverses  were 
carried  out  using  a  stepping  motor-driven  linear  positioner 
with  a  positioning  accuracy  of  5  jum.  Vertical  traverses  in¬ 
volved  moving  the  injector  assembly  along  linear  bearings 
with  a  manual  positioner  having  a  positioning  accuracy  of 
0.5  mm. 

Gamma-Ray  Absorption  Measurements 

Distributions  of  mean  liquid  volume  fractions  were  mea¬ 
sured  using  gamma-ray  absorption  similar  to  Ruff  et  al.1  An 
iodine- 125  isotope  source  (5  mCi.  emitting  primarily  at  27.20. 
27.47  and  31.00  keV)  provided  a  soft  gamma-ray  source  with 
good  absorption  levels.  The  source  was  placed  in  a  lead  casket 
with  an  outlet  aperture  of  1.6-mm  diameter  and  13-mm  long. 
Gamma  rays  passing  through  the  flow  were  detected  and 
counted  with  a  Bicron  X-ray  probe  (Model  1  x  M. 040'  1.54) 
and  a  EG&G  Ortec  single-channel  analyzer  and  counter  timer 
(models  556.  590A.  and  974).  A  lead  aperture  (1.5-mm  di¬ 
ameter  and  12-mm  long)  was  placed  in  front  of  the  detector 
to  define  the  path  observed  through  the  flow.  The  energy 
window  of  the  detector  was  set  at  22-32  keV  to  minimize 
spurious  counts  due  to  background  radiation  and  Compton 
scattering.  The  source  and  detector  were  placed  in  recessed 
mounts  within  opposing  windows  of  the  apparatus  ( not  show  n 
in  Fig.  1)  so  that  they  were  separated  by  a  distance  of  500 
mm  with  the  radiation  path  crossing  horizontally  through  the 
chamber  axis. 

Absorption  measurements  (based  on  20.CXXI-25.000  counts) 
were  made  for  30-60  parallel  paths  through  the  flow  and 
deconvolved  in  the  same  manner  as  Santoro  et  al.,;  and  Ruff 
et  al.1  The  narrow  absorption  path  minimized  potential  errors 
due  to  the  orientation  of  liquid  elements  to  less  than  5<~r.n 
Experimental  uncertainties  were  largely  due  to  finite  sampling 
times  and  background  from  the  small  drops  dispersed  within 
the  pressure  vessel.  They  are  estimated  (95^  confidence)  to 
be  less  than  30 ac  for  af  >  0.01. 

Test  Conditions 

Test  conditions  are  summarized  in  Table  1.  Operating  con¬ 
ditions  were  selected  to  yield  the  same  mean  jet  exit  velocity 
for  fully  developed  and  slug  flow  at  pressures  of  l,  2.  4.  and 
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Tabic  1  Summary  of  test  conditions' 


Jet-exit  diameter  (mm) 

9  5 

Ambient  pressure  (atm) 

1.  2.  4.  and  8 

Jet  flow  rate  (kgs) 

3  47 

Injector  pressure  drop  (kPa): 

Fullv  developed  flow 

:ro 

Slue  flow 

2I1«) 

Averaee  jet  exit  velocity  (m  si 

44  1 

Re. 

4A2.00») 

wv.. 

-12.000 

wv<a 

3 so.  'no.  (520.  3040 

Oh 

0  00(21 

JPres$ure-atomized  jet  injected  verticalK  Uo^nu^rd  in  still  jir  jt  var¬ 
ious  pressures  and  29$  r  2K;  in  atomization  breakup  retime  tor  both  due 
tlo*  and  lulls  developed  lurbulent  pipe  flow  \L  J  =  41 1  ict-exit  conditions 


S  atm  within  the  atomization  breakup  regime.  Due  to  the 
limitations  of  the  pump,  this  required  a  water  flow  rate  roughly 
13*7  lower  than  the  atomization  breakup  condition  used  by 
Ruff  et  al. '  1  for  the  9.5-mm-diam  injector  However,  all  test 
conditions  are  well  into  the  atomization  breakup  regime  de¬ 
fined  by  Ranz4  and  Miesse.'  and  exhibited  initial  liquid  breakup 
right  at  the  jet  exit,  which  is  characteristic  of  this  breakup 
regime. 

Flow  properties  at  the  jet  exit  were  measured  earlier  using 
taser  velocimetryv  For  fully  developed  flow  ,  mean  streamwise 
velocity  distributions  were  in  good  agreement  with  values  in 
the  literature  for  the  same  Reynolds  number  range.'4  How¬ 
ever.  rms  streamwise  and  radial  velocity  fluctuations  were 
more  uniform  across  the  central  region,  yielding  streamwise 
and  radial  turbulence  intensities  of  7  and  4*7  near  the  axis, 
which  are  somewhat  larger  than  literature  values.14  For  slug 
flow,  mean  streamwise  velocities  were  uniform  over  the  cen¬ 
tral  region  of  the  flow  and  then  declined  near  the  wall  (within 
3-5*7-  of  the  injector  radius)  due  to  boundary -layer  growth 
in  the  nozzle  passage.  Streamwise  and  radial  rms  veloctty 
fluctuations  were  roughly  1*7  of  the  mean  streamwise  velocity- 
over  the  central  region  for  slug  flow.' 


streamwise  mean  momentum,  mean  mixture  fraction,  tur¬ 
bulence  kinetic  energy,  and  the  rate  of  dissipation  of  turbu¬ 
lence  kinetic  energy”  The  specific  formulation,  all  empirical 
constants,  and  a  discussion  of  calibration  of  the  approach  for 
a  variety  of  constant  and  variable  density  single-phase  jets, 
appears  elsewhere. 

The  specification  of  initial  and  boundary  conditions,  and 
the  details  of  the  numerical  computations,  can  be  found  in 
Ruff  et  ai.' :  For  fully  developed  flow  ,  initial  profiles  of  u.  k. 
and  e  were  taken  from  Schlichting'4  and  Hinze  '  because  thev 
are  good  approximations  of  present  measured  jet-exit  con¬ 
ditions  and  are  readily  available  to  others.'  For  slug  flow, 
properties  were  assumed  to  be  uniform  except  for  bounding 
estimates  of  properties  in  the  boundary  layer  along  the  wall 
for  Ld  =  0  and  5.  with  the  latter  conditions  found  assuming 
clean  entry  and  no  vena  contracta  along  the  nozzle  passage 
from  Schlichting.14 


Results  and  Discussion 
Mean  Liquid  Volume  Fractions 

Experimental  Sensitivity 

The  mean  liquid  volume  fraction  distnbutions  provide  a 
quantitative  indication  of  effects  of  gas.  liquid  density  ratio 
and  jet-exit  conditions  on  flow  properties.  However,  it  is  im¬ 
portant  to  recognize  the  relationship  between  liquid  volume 
fractions  and  mixing  levels  when  interpreting  these  results 
The  present  flows  have  large  liquid  gas  density  ratios,  which 
implies  that  liquid  volume  fractions  vary  rapidly  with  mixture 
fraction.  This  can  be  seen  from  the  state  relationship  for  liquid 
volume  fraction: 


Of  =  /'■</  *  (PrPfM  l  -  f»  (-') 

Using  Eq.  (2).  an  analagous  expression  can  be  obtained  for 
time-averaged  mixture  fraction  57  as  a  function  of  Favre- 
averaged  mixture  fraction  /as  follows: 

a.  =  /</  -  (p,p,)(l  -  /))  (-M 


Theoretical  Methods 

Predictions  of  flow  properties  were  limited  to  use  of  the 
LHF  approximation  similar  to  past  work.1 ;  ’  In  addition  to 
the  LHF  approximation,  the  major  assumptions  of  the  model 
are  as  follows:  steady  (in  the  mean)  axi-vmmetric  flow  with 
no  swirl,  boundary-layer  approximations  apply,  negligible  ki¬ 
netic  energy  and  viscous  dissipation  of  the  mean  flow,  buoy¬ 
ancy  only  affects  the  mean  flow,  equal  exchange  coefficients 
of  all  species  and  phases,  and  negligible  mass  transport  be¬ 
tween  the  phases  (no  evaporation).  These  assumptions  are 
either  conditions  of  the  experiments  or  are  justified  by  past 
practice,  except  for  the  LHF  approximation,  which  will  be 
evaluated  by  the  measurements.  In  particular,  operation  of 
well-atomized  sprays  within  a  closed  container  saturated  the 
air  with  water  vapor  so  that  there  was  no  potential  for  liquid 
evaporation. 

Under  these  assumptions,  the  instantaneous  mixture  frac¬ 
tion  (defined  as  the  fraction  of  mass  that  originated  from  the 
injector). is  either  0  (in  the  gas)  or  1  (in  the  liquid).  Then, 
time-  and  Fav re-averages  of  any  scalar  property  <t>  can  be 
found  in  terms  of  the  Fav  re-averaged  mean  mixture  fraction, 
as  follows: 

d>  =  d>,(l  -  /)  *  <b,f  (1) 

Z  =  (<t>,p,{  1  -  /)  -  6..pJ)  (p,{  1  -  />  -  pJ)  (2) 

Given  Eqs  ( 1)  and  (2).  the  flowfield  can  be  found  from  a  sim¬ 
plified  version  of  the  conserved-scalar  formation  of  Lockwood 
and  Naguib.1'  but  using  Favre  averages  following  Bilger  ' 
Governing  equations  are  solved  for  conservation  of  mass. 


Table  2  is  a  summary  of  f  for  a.  =  0. 1  and  0.01  (or  equivalently 
/for  the  two  values  of  a.)  over  the  present  ambient  pressure 
range.  It  is  evident  from  the  table  that  low  levels  of  mixing 
cause  large  reductions  in  liquid  volume  fractions,  even  at  the 
highest  ambient  pressures  of  the  present  test  range.  Thus.  a. 
is  an  unusually  sensitive  indicator  of  mixing  levels  near  the 
jet  exit. 

Axial  Distributions 

Measured  and  predicted  time-averaged  mean  liquid  volume 
fractions  along  the  jet  centerline  are  illustrated  in  Figs  2  and 
3  for  fully  developed  and  slug  flow-jet  exit  conditions,  re¬ 
spectively.  Results  are  plotted  as  a  function  of  distance  from 
the  injector,  normalized  by  the  injector  diameter,  with  am¬ 
bient  pressure  as  a  parameter.  Predictions  for  slug  flow  jet- 
exit  conditions  are  shown  for  L  d  =  0  and  5  which  bounds 
the  potential  degrees  of  flow  development  within  the  nozzle 
passage,  as  noted  earlier  Measurements  of  Ruff  et  at  ;  at  1 
atm  for  atomization  breakup  also  are  shown  on  the  plots  for 
both  slug  and  fully  developed  flow,  they  agree  very  well  with 
present  results,  although  jet-exit  velocities  are  slightly  differ¬ 
ent.  This  is  expected  based  upon  the  LHF  predictions,  which 
exhibit  little  variation  of  liquid  volume  fraction  distributions 
with  jet-exit  velocity  for  the  high  Rey  nolds  numbers  of  present 
flows.1 

For  fully  developed  flow.  Fig.  2.  the  region  near  the  let 
exit  Ixd  <  3-S)  exhibits  mean  liquid  volume  fractions  near 
unity  Just  beyond  this  region,  however,  mean  liquid  volume 
fractions  decrease  rapidly  .  The  initial  reduction  of  a.  occur- 
at  progressively  smaller  v  alues  of  t  d  as  the  pressure  increase- 
with  values  of  a.,  at  a  given  value  of  r  d  generally  being  lower 
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Table  2  /via,  for  air/water  mixtures'* 


0  10  0.990  0.979  0.960  0  923 

001  0.897  0.813  0.684  0.520 


‘Air  waier  mixtures  ut  300  K  and  various  pressures 


at  higher  pressures  as  well.  These  trends  indicate  faster  mixing 
rates  at  higher  ambient  gas  densities,  analogous  to  effects  of 
flow  density  ratio  for  single-phase  turbulent  jets."  There  is 
good  agreement  between  measurements  and  predictions,  in¬ 
dicating  that  the  LHF  approach  correctly  treats  effects  of  the 
density  ratio  of  the  flow  on  mixing  properties.  However,  con¬ 
ditions  illustrated  in  Fig.  2  represent  relatively  low  levels  of 
mixing:  for  example,  results  in  Table  2  suggest  that  Favre- 
averaged  mixture  fractions  are  generally  greater  than  0.85. 
For  such  low  levels  of  mixing,  predictions  based  on  the  LHF 
approximation  have  been  reasonably  good  in  the  past.1  be¬ 
cause  separated  flow  effects  due  to  relative  velocity  differ¬ 
ences  between  the  gas  and  liquid  are  not  very  significant  when 
the  mass  of  the  flow  is  predominantly  liquid  Based  on  past 
evaluations  of  the  methodology.'  performance  of  the  LHF 
approach  is  likely  to  be  poorer  as  the  dilute  dispersed  flow 
regime  (where /«  1)  is  approached:  behavior  in  this  regime 
is  considered  later.  Finally,  although  the  variation  of  a,c  sug¬ 
gests  a  relatively  short  liquid  core,  this  is  not  the  case  when 
viewed  in  terms  of  mixture  fraction.  Favre-averaged  mixture 
fractions  are  greater  than  0.85  for  all  the  results  illustrated  in 
Fig.  2.  so  that  even  low  levels  of  flapping  of  the  liquid  core 
can  explain  the  reductions  of  a,c. 

The  slug  flow  results  illustrated  in  Fig.  3  exhibit  slower  rates 
of  mixing  than  the  fully  developed  flows.  First.  a,c  remains 
at  unity  until  x id  is  in  the  range  20-50.  Similar  to  the  results 
for  fully  developed  flow,  however,  the  value  of  x/d  where  afc 
first  begins  to  decrease  from  unity  progressively  decreases  as 
the  pressure  increases,  implying  faster  rates  of  mixing  at  higher 
ambient  densities.  The  strong  effect  of  the  degree  of  flow 
development  at  the  jet  exit  (cf.  Figs.  2  and  3)  is  similar  to 
earlier  observations  at  atmospheric  pressure.1  This  occurs  be¬ 
cause  the  liquid  density  is  large  in  comparison  to  the  gas; 
therefore,  the  fully  developed  flow  carries  significant  levels 
of  turbulence  energy  into  the  mixing  layer,  which  enhances 
mixing  rates.  Predictions  of  properties  for  slug  flow  condi¬ 
tions  are  very  sensitive  to  the  degree  of  flow  development  at 
the  jet  exit.  Thus,  there  are  significant  differences  between 
predictions  for  pure  slug  flow  (Lid  =  0)  and  allowance  for 
boundary-layer  growth  within  the  injector  passage  (Lid  =  5). 
These  conditions  bound  the  range  of  possibilities  for  present 
tests,  and  it  is  encouraging  that  the  two  predictions  tend  to 
bound  the  measurements  except  at  x/d  =  100  and  pressures 
of  4  and  8  atm.  Tiie  discrepancies  at  x/d  =  100  occur  in  a 
region  where  the  streamwise  variation  of  flow  properties  is 
rapid,  and  tends  toward  dilute  conditions.  Thus,  because  both 
predictions  tend  to  overestimate  the  rate  of  development  of 
the  flow,  separated  flow  effects  are  probably  responsible  for 
the  difficulty. 

Radial  Distribution  ( Fully  Developed  Flow ) 

Predicted  and  measured  radial  profiles  of  mean  liquid  vol¬ 
ume  fractions  for  fully  developed  flow  are  illustrated  in  Figs. 
4  and  5  for  ambient  pressures  of  1  and  8  atm  (results  at  2  and 
4  atm  are  similar).  These  results  involve  a,!a,c  plotted  as  a 
function  of  radial  distance  normalized  by  the  injector  radius, 
so  that  the  actual  width  of  the  flow  can  be  seen.  Results  are 
shown  for  various  x/d  £  100  because  larger  distances  risked 
disturbances  of  the  flow  from  the  chamber  walls.  The  mea¬ 
surements  of  Ruff  et  al.1  at  1  atm  are  illustrated  in  Fig.  4 
along  with  the  present  results.  The  two  sets  of  experiments 


x/d 

Fig.  2  Time-averaged  liquid  volume  fractions  along  the  axis  at  var¬ 
ious  ambient  pressures  for  fully  developed  flow. 
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Fig.  3  Time-averaged  liquid  volume  fractions  along  the  axis  at  var¬ 
ious  pressures  for  slug  flow. 

agree  within  experimental  uncertainties  except  at  x  d  =  100 
where  the  confinement  of  the  present  jets  may  be  responsible 
for  somewhat  reduced  flow  widths.  As  discussed  earlier,  pre¬ 
dictions  are  essentially  the  same  for  the  Reynolds  number 
differences  of  the  two  tests.  The  measurements  in  Figs.  4  and 
5  show  a  progressive  increase  of  flow  width  with  increasing 
distance  from  the  jet  exit,  with  flow  widths  increasing  at  a 
faster  rate  at  the  higher  pressure.  Apparent  flow  radii  based 
on  a/a,r.  however,  are  much  smaller  than  for  single-phase 
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Fig.  4  Radial  profiles  of  mean  liquid  volume  fractions  for  fully  de¬ 
veloped  flow  at  I  atm. 


Fig.  5  Radial  profiles  of  mean  liquid  volume  fractions  for  fully  de¬ 
veloped  flow  at  2  atm. 


jets  because  of  the  sensitivity  of  c7  to  the  extent  of  mixing, 
discussed  earlier  Notably,  flow  widths  based  on  mean  void 
fraction  distributions  for  gas  jets  in  liquids  are  unusually  large 
for  similar  reasons.1'  For  both  turbulent  liquid  jets  in  gases 


and  gas  jets  in  liquids,  however,  predictions  using  the  LHF 
approach  indicate  relatively  normal  flow  widths  far  from  the 
jet  exit,  in  terms  of  /:  this  is  discussed  more  fully  later. 

The  comparison  between  predicted  and  measured  liquid 
volume  fraction  distributions  in  Figs.  4  and  5  is  generally  quite 
good,  consistent  with  the  good  predictions  of  properties  along 
the  axis  for  fully  developed  flow  observed  in  Fig.  2.  The  main 
exception  is  the  farthest  downstream  position,  x  J  =  100.  at 
1  atm.  where  measurements  exhibit  a  wider  flow  than  pre¬ 
dictions.  Ruff  et  al.‘ :  show  that  this  difficulty  is  largely  due 
to  effects  of  separated  flow  as  the  multiphase  flow  becomes 
more  dilute.  Predictions  at  x  d  =  100  are  improved  at  S  atm. 
which  is  expected,  because  the  larger  gas  density  should  yield 
smaller  drops  after  secondary  breakup,  tending  to  reduce  er¬ 
rors  due  to  effects  of  separated  flow.:  1 

Radial  Distributions  {Slug  Flow) 

Predicted  and  measured  radial  profiles  of  mean  liquid  vol¬ 
ume  fractions  for  slug  flow  are  illustrated  in  Figs.  6  and  7  for 
ambient  pressures  of  1  and  8  atm  (results  at  2  and  4  atm  are 
similar).  Measurements  of  Ruff  et  al. 1  at  1  atm  for  slug  flow 
are  illustrated  in  Fig.  b  along  with  the  present  measurements, 
fn  this  case,  both  sets  of  measurements  agree  everywhere 
within  experimental  uncertainties,  in  accord  with  expectations 
from  the  predictions.  Better  agreement  between  the  two  sets 
of  measurements  for  slug  flow  than  fully  developed  flow  at 
xd  =  100  is  reasonable  because  the  slug  flow  is  narrower, 
which  reduces  effects  of  confinement  for  present  measure¬ 
ments.  Predictions  for  L  d  =  0  and  5  are  shown,  similar  to 
Fig.  3.  Results  illustrated  in  Figs.  6  and  7  exhibit  a  relatively 
sharp  transition  between  the  liquid  core  of  the  flow  and  the 
region  where  a,  decreases,  at  least  for  x  d  £  50.  Additionally, 
the  extent  of  radial  spread  of  the  slug  flows  is  less  than  the 
fully  developed  flows.  Both  these  observations  are  consistent 
with  slower  mixing  rates  for  slug  flow  than  fully  developed 
flow.  As  before,  however,  flow  widths  at  a  particular  x  d  are 
significantly  larger  at  8  atm  than  1  atm.  suggesting  faster 
mixing  rates  at  higher  pressures. 


Fig.  6  Radial  profiles  of  mean  liquid  volume  fractions  for  slug  flow 
at  I  atm. 
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Fig.  7  Radial  profiles  of  mean  liquid  volume  fractions  for  slug  flow 
at  2  atm. 

Predictions  of  slug  flow  properties  in  Figs.  6  and  7  are 
reasonably  good  for  x/d  £  25.  within  the  bounds  of  the  limiting 
estimates  of  the  degree  of  flow  development  at  the  jet  exit. 
Farther  from  the  injector,  however,  predictions  are  less  sat¬ 
isfactory.  In  particular,  errors  are  large  in  the  region  where 
arc  first  begins  to  decrease  from  unity  (see  Fig.  3).  This  be¬ 
havior  is  caused  by  poor  estimates  of  afe,  due  to  uncertainties 
in  the  initial  degree  of  flow  development,  because  afc  is  used 
to  normalize  both  predictions  and  measurements  in  Figs.  6 
and  7.  Beyond  this  region,  predictions  are  in  better  agreement 
with  measurements  but  this  is  largely  fortuitous  because  afc 
is  still  not  predicted  very  well  by  either  limiting  condition  (see 
F'g.  3).  These  difficulties  are  probably  due  to  the  effects  of 
separated  flow  in  the  rapidly  developing  region  near  the  tip 
of  the  liquid  core  for  slug  flow,  as  discussed  earlier. 

Sensitivity  Study 

The  sensitivity  of  present  computations  was  examined  sim¬ 
ilar  to  past  work.1 2  Predictions  were  very  sensitive  to  initial 
mean  velocity  distributions,  as  can  be  seen  from  the  results 
illustrated  in  Figs.  3.  6,  and  7  for  slug  flow  at  the  limits  Lid  = 
0  and  5.  Predictions  were  also  sensitive  to  initial  values  of  k 
and  e.  with  10%  changes  in  these  properties  causing  5-10% 
changes  in  a,r  for  x/d  in  the  range  20-60.  However,  these 
properties  were  reasonably  well  known  for  the  fully  developed 
flows,  whereas  effects  of  low-turbulence  levels  in  the  case  of 
the  slug  flow  were  not  very  significant.  Thus,  the  uncertainties 
of  the  predictions  were  largely  governed  by  difficulties  in 
specifying  the  mean  velocity  distribution  for  slug  flows,  and 
the  less  quantifiable  limitations  of  k-e  turbulence  models  using 
the  LHF  approximation  for  multiphase  boundary-layer  flows. 

Favre-averaged  Mixture  Fractions 

Axial  Distributions 

Computations  indicated  that  distributions  of  Favre-averaged 
mixture  fractions  were  much  less  influenced  by  variations  of 


Fig.  8  Favre-averaged  mixture  fractions  along  the  axis  for  fully  de¬ 
veloped  and  slug  flow. 


70 


ambient  gas  density  than  time-averaged  liquid  volume  frac¬ 
tions;  therefore,  both  present  measurements  and  those  of  Ruff 
et  al.1-3  were  examined  in  terms  of  this  fundamental  mixing 
property.  This  was  done  by  computing  /  from  the  measured 
values  of  a  using  Eq.  (4). 

It  is  well  known  that  properties  along  the  axis  of  single¬ 
phase  vanable-density  jets  should  scale  in  terms  of  a  density- 
weighted  streamwise  distance,  which  becomes  (p,'p,)' :  x  d  for 
present  measurements.'*  Chehroudi  et  al.1*  also  find  similar 
scaling  for  their  correlation  of  the  length  of  the  liquid  core. 
Thus,  predicted  and  measured  fc  are  plotted  in  terms  of  this 
variable  for  both  fully  developed  and  slug  flows  in  Fig.  8. 
Measurements  ate  identified  by  ambient  pressure  and  source. 
When  plotted  in  the  manner  of  Fig.  8.  predicted  effects  of 
ambient  pressure,  jet-exit  Reynolds  number,  and  Lid  (for  slug 
flow)  were  small;  therefore,  only  single  predicted  lines  are 
shown  for  fully  developed  and  slug  flow. 

In  contrast  to  afc.  the  sensitivity  of  present  experiments  to 
/c  is  relatively  low  due  to  the  large  density  ratios  of  the  flows. 
Thus,  measured  values  of/c  are  generally  greater  than  0.85 
as  noted  in  connection  with  the  plots  of  afc  (Figs.  2  and  3). 
Over  this  narrow  range  of /.,  however,  the  measurements  at 
various  ambient  pressures  correlated  reasonably  well  in  terms 
of  density-weighted  streamwise  distance. 

Predicted  and  measured  fc  shown  in  Fig.  8  are  in  reasonably 
good  agreement  for  (p/p/)1 :  x/d  <  5.  However,  this  obser¬ 
vation  is  much  less  definitive  than  comparisons  between  pre¬ 
dicted  an  teasured  values  of  afc  due  to  the  reduced  sensi¬ 
tivity  of/c.  similarly,  the  striking  effect  of  jet-exit  turbulence 
on  the  variation  of  afc  is  much  less  evident  for  fc.  The  meas¬ 
urements,  however,  yield  larger  values  of  /c  than  predicted  at 
density-weighted  distances  greater  than  five,  particularly  for 
slug  flow  where  predicted  properties  vary  most  rapidly  in  the 
streamwise  direction.  This  implies  slower  rates  of  mixing  along 
the  axis  than  LHF  predictions,  suggesting  significant  effects 
of  separated  flow  in  the  region  just  downstream  ot  the  end 
of  the  liquid  core.  This  is  plausible,  because  breakup  of  the 
end  of  the  liquid  core,  the  same  as  primary  breakup  along  its 
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surface,  should" yield  larger  drops  having  significant  relative 
velocities.  Thus,  the  predicted  universal  behavior  of  fc  was 
not  approached  tor  present  test  conditions  as  the  flow  become 
more  dilute  for  density-weighted  streamwise  distances  greater 
than  five.  This  is  much  more  evident  from  the  plots  of  /.  in 
Fig.  8  than  those  for  afc  in  Figs.  2  and  3  due  to  effects  of  large 
density  ratio  on  sensitivity  to  these  properties. 

Radial  Distributions 

Due  'o  the  length  of  the  liquid  core,  present  flows  are 
transitional  between  a  mixing  layer  and  a  jet  when  viewed  in 
the  radial  direction.  Thus,  radial  profiles  of/ were  limited  to 
mixing  layer  conditions  (/,  =  1 )  to  avoid  the  complications  of 
this  transition.  Within  the  mixing  layer,  predictions  showed 
that//  was  relatively  independent  of  the  density  ratio  when 
plotted  in  terms  of  ( r  -  r,Y:,t)<x  so  that  measurements  will 
be  considered  in  terms  of  these  variables  In  order  to  extend 
the  range  of  /  that  could  be  exanmed.  the  present  gamma- 
ray  absorption  measurements  were  supplemented  bv  holog¬ 
raphy  measurements  in  the  multiphase  mixing  laver  from  Ruff 
et  al.: 

Predicted  and  measured  radial  profiles  of  //  are  plotted  in 
Fig.  9  for  both  fully  developed  and  slug  flows  As  mentioned 
earlier,  predictions  exhibit  small  changes  with  densitv  ratio 
and  position  for  present  test  conditions.  The  present  gamma- 
ray  absorption  measurements  are  in  reasonably  good  agree¬ 
ment  with  predictions  when  plotted  in  the  manner  of  Fia  9. 
however,  they  are  limited  to  //  >  0  85. 

The  agreement  between  predictions  and  the  holography 
measurements  illustrated  in  Fig.  9  is  much  less  satisfactory. 
The  holography  measurements  were  obtained  in  the  dispersed 
flow  region  where  the  liqi.  ;d  is  present  as  irregular  liquid 
elements  and  drops,  which  is  typical  of  dilute  dispersed  flows. 
Thus,  poorer  performance  of  the  LHF  approximation  is  ex¬ 
pected.  based  on  observations  in  other  dilute  dispersed  flows.’' 


Fig.  9  Radial  profiles  of  Fas  re-averaged  mixture  fractions  in  the 
multiphase  mixing  layer  for  fully  developed  and  slug  flo’v. 


Nevertheless,  the  holographic  measurements  for  slug  flow, 
conditions  do  approach  predictions,  with  predicted  overesti- 
mation  of  the  mixing  rates  near  the  outer  edge  of  the  flow 
being  attributable  to  separated  flow  effects.  Discrepancies 
between  predictions  and  measurements,  however,  are  much 
larger  for  fully  developed  flow.  Nevertheless,  this  behavior  is 
reasonable  because  drop  sizes  after  primary  breakup  were 
much  larger  for  the  fully  developed  flow  than  the  slug  flow 
(SMD  of  1330  p. m  as  opposed  to  174  pm>.  which  implies 
greater  effects  of  separated  flow.  -  Thus,  improved  atomiza¬ 
tion  conditions  might  yield  the  universal  results  suggested  by 
the  predictions:  however,  additional  measurements  are  needed 
to  evaluate  this  possibility. 

A  final  matter  to  be  considered  in  connection  with  Fig.  9 
is  why  the  holography  measurements  show  faster  mixing  than 
predicted  for  fully  developed  flow  but  slower  mixing  than 
predicted  for  slug  flow.  This  behavior  is  caused  by  interactions 
between  primary  breakup  and  separated-flow  phenomena.' 
Primary  breakup  of  a  turbulent  liquid  generates  large  liquid 
elements  with  radial  velocities  comparable  to  radial  velocity 
fluctuations  in  the  liquid.  Thus,  these  large  liquid  elements 
are  projected  across  the  dispersed  flow  before  they  break  up 
or  relax  to  local  flow  velocities,  which  enhances  mixing.  In 
contrast,  primary  breakup  of  nonturbulent  liquids  yields  small 
radial  velocities,  while  large  drops  have  poor  response  to 
turbulent  dispersion  due  to  their  inertia,  which  retards  mixing. 


Conclusions 

The  structure  and  mixing  properties  of  nonevaporating  liq¬ 
uid  jets  in  still  gases  were  studied  for  various  ambient  gas 
densities  at  atomization  breakup  conditions,  emphasizing  the 
region  near  the  jet  exit.  The  major  conclusions  of  the  study 
are  as  follows: 

1)  Increasing  gas/liquid  density  ratios  reduces  the  length  of 
the  liquid  core  and  increases  the  flow  width,  implying  in¬ 
creased  rates  of  mixing  analogous  to  effects  of  density  rai:o 
for  single-phase  turbulent  jets. 

2)  Turbulence  levels  and  the  degree  of  flow  development 
at  the  jet  exit  have  a  strong  effect  on  mixing  rates,  with  tur¬ 
bulent  flows  mixing  much  faster  than  nonturbulent  slug  flows 

3)  Use  of  the  locally  homogeneous  flow-  approximation,  in 
conjunction  with  a  Favre-averaeed  turbulence  model,  yielded 
good  estimates  of  effects  of  gas  liquid  density  ratio  and  initial 
liquid  vorticity  on  time-averaged  liquid  volume  fraction  dis¬ 
tributions  and  Favre-averaged  mixture  fraction  distributions 
in  the  region  near  the  liquid  core  where  /  >  0.85.  The  main 
reason  for  this  is  that  separated  flow  effects,  like  small  ve¬ 
locities  of  the  gas  and  small  drops,  are  not  very  important  t'oi 
these  low  degrees  of  mixing.  Nevertheless,  it  is  helpful  that 
predictions  provide  reasonable  estimates  of  the  initial  stages 
of  mixing. 

4)  Distributions  of  /  were  relatively  independent  of  the 
density  ratio  of  the  flow  in  terms  of  (p,  p,)' 1  x  d  along  the 
axis  and  (r  -  r,,05).x  across  the  mixing  It"  er  for  /  >  0.85. 
therefore,  many  of  the  unusual  mixing  properties  of  these 
flows  in  terms  of  a,  are  due  to  the  strong  nonlineanties  of  the 
state  relationship  between  a,  and  /  rather  than  unusual  be¬ 
havior  of  the  fundamental  mixing  variable  /  itself  Pi  edictions 
suggest  that  this  scaling  of  /  should  extend  to  the  full  range 
of  /if  effects  of  separated  flow  are  small:  however,  this  regime 
was  not  reached  during  the  present  experiments 

5)  For  /  <  0.85.  significant  effects  of  separated  flow  were 
observed,  causing  poor  performance  of  predictions  using  the 
locally  homogeneous  flow  approximation  and  effects  of  the 
density  ratio  on  the  scaling  of  distributions  of  /  This  agrees 
with  other  determinations  of  important  effects  of  separated 
flow  in  dilute  dispersed  flows,  however,  quantitative  estimates 
of  conditions  where  separated-flow  must  be  considered  re¬ 
quire  more  information  on  primary  breakup  properties  than 
■  currentlv  available. 
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Present  conclusions  are  based  on  large-scale  sprays  (9.5- 
mm  injector  diameter)  that  have  much  lower  rates  of  decelera¬ 
tion  than  practical  injectors,  and  regions  of  the  flow  having 
relatively  high  mixture  fractions  (generally  greater  than  0.9 
along  the  axis)  where  the  momentum  of  the  gas  does  not  have 
a  strong  influence  on  flow  dynamics.  These  factors  favor  use 
of  the  locally  homogeneous  flow  approximation,  so  that  pres¬ 
ent  observations  are  not  necessarily  in  conflict  with  earlier 
work  showing  significant  separated-flow  effects  within  dense 
sprays  for  smaller  injector  diameters  and  mixture  trac¬ 
tions.1  " 1 '  Additional  information  concerning  liquid  breakup 
properties  in  dense  sprays  clearly  is  needed  in  order  to  provide 
a  rational  means  of  evaluating  separated-flow  effects  and  the 
adequacy  of  the  locally  homogeneous  flow  approximation  tor 
particular  conditions. 
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Dispersed-Phase  Structure  of  Pressure-Atomized  predict  observed  trends  of  decreasing  drop  sizes  with  increasing  gas  density  due  to 

Sprays  at  Various  (Jas  Densities  aerodynamic  effects;  in  contrast,  the  laminar  primary  breakup  predictions  successfully 

treated  the  relatively  small  effects  of  gas  density  for  this  breakup  mechanism,  liffects  of 
by  liquid  turbulence  at  the  jet  exit  were  qualitatively  similar  to  single  phase  flows,  yielding 

faster  mixing  rates  with  increased  turbulence  levels  even  though  drop  sizes  tended  to 
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INTRODUCTION  This  has  included  measurements  of  liquid  volume  fraction  distributions  using  gamma  ray 

absorption  for  water  jets  in  air  at  various  ambient  pressures.1"*  These  results  showed  that 
The  present  study  extends  past  work  in  this  laboratory  concerning  (he  properties  of  mixing  rates  were  influenced  by  the  breakup  regime,  turbulence  properties  at  the  jet  exit  and 

the  near  injector  dense  spray  region  of  pressure-atomized  sprays.1  6  The  first  phase  of  the  the  ambient  gas  density.  Predictions  based  on  the  locally-homogeneous  flow  (LHF) 

investigation  involved  measurements  of  liquid  volume  fractions  and  dispersed-phase  approximation,  where  relative  velocities  between  the  phases  are  neglected,  were  in  good 
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Past  work  in  this  laboratory  has  largely  concentrated  on  the  properties  of  the  near-  breakup  by  observing  water  jets  in  still  air  at  pressures  of  1,  2  and  4  atm.  with  both 

injector  region  during  atomization  breakup,  using  relatively  large  diameter  liquid  jets  (3.5  -  turbulent  and  nonturbulent  jet  exit  conditions.  Holography  was  used  to  find  drop  sizes  and 

19  I  mm)  to  provide  known  jet  exit  flow  conditions  and  reasonable  spatial  resolution  velocities,  and  liquid  volume  fractions  and  fluxes,  within  the  dispersed-flow  region  beyond 

because  the  flow  diameter  is  not  thought  to  affect  dense  spray  structure  appreciably.1'6  the  liquid  surface.  Following  earlier  work,1-4  predictions  of  flow  properties  using  the  LIIF 


approximation  were  used  !o  help  assess  separated-flow  effects  on  the  structure  of  the  flow.  measured  and  the  maximum  and  minimum  diameters  of  the  ellipsoid  having  the  same 

Ihe  correlations  for  primary  breakup  properties  also  were  evaluated  for  effects  of  varying  crossectional  area  and  perimeter  were  computed  to  find  the  effective  sphere  diameter  as 

ambient  gas  density.  The  present  discussion  is  brief;  additional  details  and  a  complete  before.  The  ellipticity,  defined  as  the  ratio  of  the  maximum  and  minimum  ellipsoid 

tabulation  of  data  can  be  found  in  Tseng  diameters,  also  was  computed  for  each  object. 

Experimental  Methods  Drop  velocity  measurements  were  based  on  the  motion  of  the  centroid  of  the  image 


the  injecior  radius),  due  io  boundary  layer  growth  in  the  nozzle  passage,  while  rms  velocity  Tseng23  for  similar  plots  at  all  test  conditions.  Figures  1-3  illustrate  flow  properties  at 

fluctuations  were  roughly  1%  of  the  mean  streamwise  velocity  over  the  central  portion  of  pressures  of  1,  2  and  4  atm.  at  the  fined  streamwise  position  of  x/d  =  6.  Figures  3  and  4 

the  flow  For  fully-developed  flow,  mean  velocity  distributions  were  in  good  agreement  illustrate  effects  of  increasing  x/d  from  6  to  25  at  a  fixed  pressure  of  4  atm.  The  plots 

with  literature  values  for  the  same  Reynolds  number  range,  while  rms  velocity  fluctuations  include  volume-averaged  ellipticity,  ep,  SMD  and  drop  velocities  for  dp  =  10,  50,  100  and 

near  the  axis  were  somewhat  larger  than  literature  values.25  26  200  pm.  The  range  of  positions  where  the  surface  of  the  liquid  core  was  observed  and 
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median  diameter  (MMD)  to  the  5MU  typical  ot  past  oDservations,-'-'1''^'  ivnviu/aiviu  -  distance  rrom  tne  jet  exit  mat  is  characteristic  ot  turbulent  primary  breakup.0  Smaller 

I  2  for  present  turbulent  and  nonturbulent  flows  so  that  the  entire  drop  size  distribution  will  relative  velocities,  caused  by  faster  mixing  rates  at  elevated  pressures,  also  are  probably  a 

be  represented  by  the  SMD  in  the  following.  factor  in  this  behavior. 

Typical  distributions  of  drop  properties  across  the  mixing  layer  for  fully-developed  Distributions  of  drop  velocities  in  Figs.  1-4  have  been  normalized  by  the  mass- 

flow,  illustrating  flow  properties  for  turbulent  primary  breakup,  appear  in  Figs.  1-4  See  averaged  velocity  at  the  exit  of  the  jet.  In  general,  there  are  substantial  differences  between 
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bound  for  gas  velocities  because  their  velocity  relaxation  limes  are  too  large  for  them  to  be  breakup. 

accurate  seeding  particles  to  find  gas  velocities.  However,  this  effect  is  not  very  significant  The  LHF  predictions  of  u  are  generally  not  very  satisfactory  in  Figs.  5-8,  because 

because  the  large  drops  contain  most  of  the  momentum  of  the  flow.  Favre-averaged  present  measurements  emphasize  near-injector  conditions  where  effects  of  separated  flow 

remain  important.  Nevertheless,  they  tend  to  improve  with  increasing  distance  from  the  jet 


exit  and  ambient  pressure,  in  accord  with  earlier  scaling  arguments.5  The  tendency  toward  'he  entire  drop  size  distribution  similar  to  Ref.  3,  however,  Weg  was  found  from  the 

better  agreement  of  LHF  predictions  comes  about  because  large  drops  generally  dominate  measured  values  of  phase  velocities  and  SMD,  as  follows: 

the  momentum  content  of  the  tlow,  similar  to  behavior  at  larger  liquid  volume  fractions  ^eg  =pg(Uf .  ug)2SMD/a  (2) 

where  the  liquid  core  is  present  2  4  The  small  drops  and  the  gas  all  have  velocities  that  are  Th's  definition  has  the  advantage  that  maximum  values  of  individual  drop  Weber  numbers 


supported  by  the  relatively  small  effect  of  ambient  pressure  on  liquid  surface  properties.  behavior  anticipated  from  the  deformation  criterion:  deformed  drops  were  quite  prevalent 

Secondary  Breakup  near  <he  liquid  surface  while  drops  near  the  edge  of  the  (low  were  generally  spherical.  This 

Similar  to  Ruff  at  al.,}  the  propensity  for  secondary  breakup  in  the  mixing  layer  behavior  contributes  to  decreasing  values  of  ep  with  radial  distance,  approaching  unity  near 

was  assessed  by  computing  mass-averaged  Weber  numbers.  Rather  than  summing  over  lhe  edge  of  the  flow,  although  the  presence  of  irregular  ligament-like  drops  contributes  to 

large  values  of  ep  near  the  surface  as  well. 


icreased  gas  densities  than  those  controlling  nonturbulent  primary  breakup  (see  Eq.  (9)).  Ignore  the  presence  of  the  liquid  surface  and  attached  ligaments,  as  noted  earlier;  and  in 

Thus,  higher  pressures  than  considered  here  may  involve  merging  of  primary  and  part,  because  of  the  significant  effects  of  separated  flow  in  the  dispersed  flow  region, 

secondary  breakup,  similar  to  present  observations  of  turbulent  primary  breakup:  this  The  liquid  flux  measurements  and  predictions  illustrated  in  Figs.  11  and  12  are 

possibility  should  be  explored.  qualitatively  similar  to  the  liquid  volume  fractions.  Fluxes  for  fully-developed  flow 
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outer  edge  of  the  region  where  the  liquid  surface  is  present,  which  involves  liquid  volume  turbulent  primary  breakup  was  in  excellent  agreement  with  present  measurements  at 

fractions  greater  than  0  1  Thereafter,  the  LUF  predictions  only  bear  a  crude  resemblance  atmospheric  pressure  —  the  same  ambient  pressure  used  to  develop  the  correlation, 

to  measured  properties  in  the  dispersed  flow  region:  in  part  because  the  measurements  However,  the  correlation  did  not  represent  measured  trends  of  reduced  SMD  after 


aerodynamic  effects,  or  merging  of  primary  and  secondary  breakup,  that  were  not  suggest  that  turbulence  generation  by  drops, 7'S  is  a  major  feature  of  the  dispersed 

considered  when  the  comdauon  was  developed.  The  difficulty  is  largely  associated  flow  region. 

with  an  effect  of  density  ratio;  therefore,  the  correlation  of  Ref.  6  should  only  be  7.  Predictions  based  on  the  LflF  approximation  were  only  qualitatively  useful  due  to 

used  for  pf/Pg  >  580,  which  was  the  range  originally  used  to  define  it.  significant  separated-flow  effects  in  the  dispersed  flow  region.  Performance  of  the 

Measured  trends  of  effects  of  ambient  pressure  and  distance  from  the  injector  on  LHF  approach,  however,  tended  to  improve  with  increasing  liquid  volume 
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regime  were  comparable  in  ihe  initial  jcl  iliamcler  anil  wcic  probably  helical  waves  as 
oR  0Jnwc.  6^  "lAcil  eailier  Ihen  assiiining  llial  Ihe  SMI)  alter  primary  breakup  is  pmporlional  In  the 

lJ"'3  thickness  nl  Ihe  siiitaec  hnumlary  layer,  lhal  the  distance  ihe  boundary  layer  grows  is 

pioportional  to  ihe  |el  cut  diameter.  lhal  Ihe  relative  veins  Hy  nl  Ihe  layer  scales  as  u,„ 
ami  ill. 11  Ihe  surface  hnumlary  layer  is  laminar.  Iliere  results 
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APPENDIX  F:  WU  ET  AL.  (1992) 


PRIMARY  BREAKUP  IN  GAS/LIQUID  MIXING  LAYERS 
FOR  TURBULENT  LIQUIDS 


105 


t  Author  to  whom  correspondence  should  be  sent. 


The  paper  begins  with  a  description  of  experimental  methods.  Results  are  then  sample  size.  These  relatively  short  time  periods  were  not  problematical,  however,  because 

considered  treating  flow  visualization,  drop  size  distributions,  drop  velocutes.  onset  of  ^ow  development  times  are  generally  less  than  1%  of  the  flow  time  for  the  near-injector 

breakup  and  drop  sizes  and  shapes,  in  turn.  The  paper  concludes  with  discussion  of  region  where  measurements  were  made  while  pulsed  photography  and  holography  only 

potential  aerodynamic  effects  on  turbulent  primary  breakup,  based  on  recent  measurements  required  times  in  the  range  I  -32  ps  for  triggering  and  data  accumulation*  Jet  exit  velocities 

in  air  at  elevated  ambient  pressures  due  to  Tseng  et  al.(  I2|.  time  of  the  measurements  were  calibrated  using  an  impact  plate  as  described  in  Ref. 
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dominated  by  finite  sampling  limitations,  with  the  accuracy  of  finding  the  centroid  of  large  This  provided  the  following  ranges  of  fluid  properties:  liquid/gas  density  ratios  of  590- 

incgular  liquid  elements  ah  being  a  factor  for  the  relatively  small  radial  velocities  that 
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photographs  are  shown  one  right  at  the  jet  exit  (which  can  be  seen  at  the  top  of  this 

wave  number  range  of  the  turbulence  spectrum  disappearing  first,  tending  to  reduce  fine- 

photograph)  in  the  region  of  initiation  of  breakup,  and  the  other  two  centered  at  x/d  =  10  ... 

grained  distortion  of  the  liquid  surface  accordingly.  It  seems  likely  that  similar  processes 
continue  to  act  even  in  the  presence  of  primary  breakup  yielding  the  properties  of  turbulent 
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pressure  with  d  =  6.4  mm,  showing  the  variation  of  drop  sizes  and  the  distance  from  the  jet  ihf  srrMmu/it»  k„  i  ,  . .  .  .  ,  ,  - 

v  lne  streamwise  direction  by  the  local  mean  velocity,  which  is  taken  to  be  u„  based  on  the 

exit  where  breakup  begins,  SMD,  and  Xj.  as  a  function  of  the  mean  jet  exit  velocity. 

Correlations  of  these  properties  found  from  Eqs.  (5)  and  (1 1)  for  A  =  d/8  from  Hinze  (14), 
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where  only  proportionally  is  implied  to  account  for  effects  of  ellipticity.  nonuniform  vicw  of  the  approximations  used  to  develop  the  correlating  express, on  and  expenmental 

velocities  within  the  eddy  and  the  efficiency  of  conversion  of  kinetic  energy  into  surface  uncertainties.  The  coefficient  of  Eq  (5)  is  relatively  large  but  this  can  be  anticipated  from 

energy.  In  order  for  drops  to  be  formed  at  all.  0,  must  be  less  than  the  largest  eddies  Eq.  (4)  if  CS1  is  on  the  order  of  unity  because  (i^o')6'5  is  large  for  fully-developed 

present,  which  is  comparable  to  A  1 10),  while  Ik  <  1,  by  definition  for  this  primary 


turbulent  pipe  flow  1 14|.  Finally.  SMDj  estimates  from  Eq.  (5)  are  seen  to  be  in  excellent  where  onset  of  breakup  occurs,  and  is  supported  by  the  appearance  of  the  (low;  therefore, 

agreement  with  the  specific  measurements  for  water  illustrated  in  Fig.  4.  this  breakup  dine  will  be  used  in  the  following. 
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present  test  range  (note  that  A  =  d/8  has  been  used  to  plot  Eq.  ( 16)  in  Fig.  8,  as  before) 
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correlations  reported  here  heyond  the  present  test  range. 
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Fig.  2  Hologram  reconstructions  of  flow  near  the  liquid  surface  at  various  relative 
velocities  for  water  injection  into  still  air. 
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Pulsed  shadowgraph  photography  also  was  used  to  find  the  streamwise  location  of  because  primary  breakup,  particularly  near  the  onset  of  breakup,  yields  relatively  few 

the  onset  of  turbulent  primary  breakup.  The  experimental  uncertainties  (95%  confidence)  ol  drops.  Within  the  limitations  of  the  definitions  of  object  sizes  and  ellipticiues.  experimental 
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icsi  conditions  Thus,  present  experiments  involve  primary  breakup  at  the  liquid  surface  complete  their  growth  faster  than  large  disturbances  and  should  be  observed  sooner 

rather  than  conditions  after  breakup  of  the  enure  liquid  volume  itself  Secondly,  reduced  levels  of  turbulence  production  as  liquid  velocities  become  uniform. 

once  the  retarding  effect  of  the  injector  passage  walls  is  removed,  cause  the  liquid 
lei  exit  velocities  were  in  the  range  22-67  m/s,  yielding  the  following  ranges  of  jei  turbulence  10  decay  with  the  small-scale  high  wave  number  range  of  the  turbulence,  and  ihe 
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l»  unary  breakup  roughly  correspond  to  ihe  sixes  of  ligaments  and  other  protrusions  from  the  orientation  of  the  experiment  although  effects  of  gravitauonaJ  forces  were  negligible  due 

ihe  liquid  surface  Additionally,  the  sixes  of  surface  protrusions,  and  corresponding  drops.  to  the  high  liquid  velocity  The  liquid  core  of  the  jet  is  toward  ihe  left  of  each  photograph 

progressively  increase  with  increasing  distance  from  the  jet  exit  for  a  given  injection  Finally,  a  900  pm  diarneier  pm  is  visible  in  all  dure  p holographs  in  order  to  prov.de  a  s./e 

voikIiiioii  I  his  behavior  was  attributed  to  two  effects  |2|  First  of  all.  small  disturbances  reference 


Other  features  of  the  photograph  in  Fig.  2  at  the  lowest  liquid/gas  density  ratio 

The  photographs  o.  Fig  2  illustrate  some  of  the  mam  features  of  aerodynamic  SUPP°"  P«*"«  of  significant  aerodynamic  effects  at  this  coodiuon  as  well  fo, 

phenomena  on  turbulent  primary  breakup  F.rst  of  all.  it  is  evident  fwm  the  left  and  center  «■*  »**"  S'neni'*  «  denec,ed  u>**nl  ,he  *  «"•  “  »"UC,P»*«1  whe"  *»* 

. . uni.,  . . .  . . .  out  rn  an?  h.«  vifiualtv  nn effect  on  nronenier  nnr  forces  are  more  significant,  rather  than  the  random  orientations  seen  at  the  two  higher 


known  ,f  the  SMD  is  known  Thus,  drop  sues  w.ll  be  descnbed  .n  .eons  of  (he  SMD  alone  breakup  of  20-40*  |20|  are  consistent  with  stmtlar  reductions  of  sirean.w.se  drop 

in  (he  following,  similar  to  pas.  work  1 1.  2.  4.  6.  7|  velocities  dunng  present  observations  for  conditions  where  aerodynamic  effects  are 

significant. 
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IH  -  ip,  SMDVo)'/2  Similarly,  by  analogy  10  ihe  results  of  Ref  19  and  27.  ihe  10  play  a  role  in  reducing  drop  sires  near  onset  conditions  of  turbulent  primary  breakup 

characteristic  secondary  breakup  time  was  taken  to  be.  fR  -  (p,/pg)'«  SMD/uj  Then  when  Itquid/gas  density  ratios  are  less  than  500.  Furthermore,  the  rather  close 

using  Cq  <  1 8|  to  eliminate  SMD  from  the  ratio,  the  characteristic  time  rauo  was  taken  to  be  correspondence  between  effects  of  secondary  breakup  on  the  primary  breakup  process 

observed  here,  and  secondary  breakup  of  drops  exposed  to  shock  wave  disturbances 


ticguis  can  be  explained  by  interactions  among  liquid  turbulence,  Rayleigh  breakup  ol  neaf  region  0f  large  scale  pressure-atomired  sprays  were  studied,  considering  liquid 

ligaments  and  fast  (merged)  secondary  breakup.  Additionally,  the  classical  aerodynamic  m  it|„  gases  w|(h  fully  developed  turbulent  pipe  flow  at  ihe  jct  exit  The  major 

dice  i  that  causes  a  reduced  pressure  over  protuberances  from  ihe  surface  (21.23),  appears  conclusions  of  ihe  siudy  are  as  follows 


The  presence  of  aerodynamic  phenomena  for  turbulent  primary  breakup  largely  is  averaged  drop  velocities  after  aerodynamic  turbulent  primary  breakup  approximate 

controlled  by  the  liquid/gas  density  ratio,  and  affect  both  conditions  at  the  onset  of  mean  and  rms  velocity  fluctuations  of  the  liquid  in  the  streamwise  and  crosstream 

breakup,  and  drop  sizes  and  velocities  (to  a  lesser  extent)  after  breakup,  when  this  directions,  respectively,  although  there  was  a  tendency  for  streamwise  velocities  to 


=  tune-averaged  r  m  s.  fluctuating  property 
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